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INTRODUCTION 


Measurements of fluid flow speed may be made by utilizing the 
Doppler shift of laser light scattered from small particles suspended 
in the flowing fluid. The principle of the Doppler shift is of course 
well known, but only recently was a technique introduced by Yeh and 
Cummins (1964) to utilize the Doppler shift of a laser radiation to 
successfully measure fluid flow speeds. Since that time there have 
been a number of separate investigations reported in the literature 
(see references). The instrument utilized in this investigation was 
developed by a team of scientists at NASA/MSFC (Huntsville, Alabama), 
Raytheon Company (Sudbury, Massachusetts) and Lockheed Missiles and 
Space Company (Huntsville, Alabama). Much of the technology used was 
originally developed in assembling a system to be used in subsonic and 
supersonic gas flows with large quantities of particle entrainment 
[Rolfe et al. (1968)]. The system used in this study involved only 
aerosols and particulate matter suspended naturally in the atmosphere. 

Interest in application of the instrument has broadened currently 
(1972) to a variety of practical situations where a remote-sensing 
instrument has particular advantages over conventional velocimeters. 

Two applications currently under research is for use as an airport 
warning system for wake vortex detection and as an air-borne system for 
clear-air turbulence detection. A potentially important use of the 
instrument is in meteorological investigations of the atmospheric 
boundary layer. Further uses of the instrument could be for remote 
air-pollution detection and for measurement of mass and momentum fluxes 
in a variety of fluid flow fields. 
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In principle it is possible to measure "point" velocities in the 
flow field with complete vector directional resolution. A laboratory 
three-dimensional instrument is presently being investigated at 
NASA/MSFC (Huntsville, Alabama), where also an atmospheric three- 
dimensional arrangement is under research and development. The instru- 
ment used in this investigation was a one-dimensional co-axial system, 
using a 25-watt C0 2 laser and back-scattered radiation. The direction 
of wind velocity was resolved by utilizing an ordinary wind-vane 
direction sensor. 

The purpose of this research project was to obtain measurements 
of atmospheric velocities and turbulence with the laser Doppler system 
and to compare the results with cup anemometer and hot-wire measurements 
in the same wind field. 


BASIC PRINCIPLES 

The frequency of laser light scattered by moving particles in a 
flow field is shifted by the Doppler effect. The Doppler shift is 
detected by optical mixing of the emitted or incident and scattered 
beams. A variety of optical configurations is possible to accomplish 
the optical mixing. In the present arrangement the back-scattered 
radiation along the axis of the incident beam was redirected into the 
laser to combine with the original laser beam. The resultant hetero- 
dyne or "beat" frequency is equal to the difference in frequencies of 
the emitted and scattered frequencies, and is directly proportional to 
the particle speed. If the scatterers are small, and no relative velo- 
city exists between the particle and the fluid, then fluid velocity is 
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measured. An infrared detector was used to convert the Doppler-shifted 
frequency to a measurable electrical signal. The arrangement of the 
system is shown schematically in Figure 1. 

The laser Doppler velocity measurement system (hereinafter referred 
to as the laser Doppler velocimeter and mnemonically denoted LDV) is 
almost instantaneous and has the advantage that no prior calibration is 
required as with other velocity instruments. The range of detectable 
velocities is very large. There is minimal perturbation of the fluid 
flow field by the laser radiation. The spatial resolution which is 
fixed ultimately by diffraction limitations can be controlled to a large 
degree by size and optical quality of the lenses and mirrors. 

A nonrelativistic derivation of velocity determination from the 
Doppler shift frequency follows. A definition diagram relative to the 
derivation is shown in Figure 2. For purpose of clarity, the scattered 
beam is shown at an arbitrary angle e from the direction of particle 
motion. In the case of a coaxial system, e = a. 

The emitted monochromatic laser radiation of wave length a q and 
speed c illuminates a particle having a velocity V. The direction of 
the incident beam is defined by the unit vector r Q . If the particle 
is motionless, the number of waves incident on the particle per unit of 
time is f^ = c/A where c is the speed of the laser radiation and 
is the wave length. 

If the particle is in motion at an angle a with respect to the 
incident beam, the frequency of the waves per unit of time relative to 
the moving particle is 

, c + VcOSd 
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Figure 2. Definition diagram for Doppler shift frequency. 


which is also the frequency of the scattered waves relative to the 
particle. The scattered radiation is directed toward a fixed point 
along a direction r s from point P. The frequency of the scattered 
radiation relative to the particle is fp, but to a fixed observer 
along r s> the wave length appears to be 

c - Vcose c - Vcose 

s T C + VC0Sa A 0 

P 


and the frequency of the scattered radiation appears to be 

f - C _ C /C + Vcosou 

T s " x "a x - Vcose' 
s o 

which is rearranged to give 


f 


s 


1 + 


'o 1 


VcOSa 
C \ 

Vcose ' ' 


c 
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The apparent shift in frequency, the Doppler shift, is 


or, 


f n = y — [ V (cosa + cose)] 
u o 

using the approximation that « 1. 

For backscatter along the incident laser beam, e = a, thus 

£ _ 2VC0Sa 

° ». 


and 

x f n _ f n 

u - o D - c D 

V 2cOSa 2coSa f Q 

In particular the component of the particle velocity along the laser 
beam axis V Q is always determinable from 

x f n cf n 

= VcOSa = ^ • 

0 2 2f 0 

The wavelength of the C0 2 laser was 10.6 microns, thus the velocity is 
given by 

V = 5.3 x 10~ 6 f n m/sec 

o u 


or. 


V Q = .53 cm/sec/KH2 Doppler shift. 


DESCRIPTION OF THE LASER DOPPLER VELOCIMETER 

The optical configuration of the LDV is shown schematically in 
Figure 1. It consists of a 25-watt, 10. 6u, C0 2 laser, beam splitters, 
mirrors and attenuators, an f8 12-inch Newtonian telescope and a 
liquid-helium cooled Ge-Hg infrared detector. 
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Based on relative power of 100 percent of the laser output 
(nominal 25 watts), the power at the focal region F was about 60 per- 
cent. The focal region is the sample space or volume from where the 
scattered signal is effectively heterodyned. The relative power at 
the detector was about 1 percent. 

The laser radiation is focused at the desired range by a 2-in. 
focusing lense L. A diagonal , 1-7/8 by 2-21/32 inches mounted on a 
spider within the 15-in. diameter tube of the telescope, directs the 
beam to a 12-in. diameter schlieren mirror mounted at the end. The 
mirror is adjustable on a 3-point mount. Physical limitation of the 
focusing lense movement limited the near range of the telescope focus 
to about 60 feet from the mirror. The other limit of the telescope 
focusing range is limited to about 250 feet by the size of the diagonal. 
Of course if the power loss from beam "spill over" at the diagonal is 
not of concern the range can be extended. A curve of focal distance as 
a function of lense movement is shown in Figure 3. The reference posi- 
tion of the lense is arbitrary and made relative to 60 feet in the 
figure. The range of the telescope relative to "performance" is also 
diffraction limited [cf. Lockheed Missiles and Space Company (LMSC) 
progress report D162417, July 23, 1970]. 

Spatial Resolution 

The spatial resolution of the system is specified in terms of the 
signal -to-noise, S/N, ratio. A calculation of S/N was made by LMSC 
(cf. Appendix A, Interim Report D225028, June 1971). It has been shown 
[Thomson and Dorian (1967)] that only radiation scattered from the 
region near the focus of the telescope contributes most significantly 



Displacement of Lens (in.) 



Figure 3. Range positioning as a function of lens position 
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to the Doppler signal. Nevertheless, there is some amount of heterodyned 
signal attributable to scattered returns in the whole space of illumina- 
tion. The ratio of S/N from the focal region in comparison to the total 
S/N, then, is a method of defining the spatial resolution. A curve of 
spatial resolution (axial dimension) as a function of focal range is 
reproduced from the LMSC calculations as Figure 4. 

Signal Processing 

There are several options for discriminating the Doppler shift in 
frequency from the detector. These are: 

1. Spectrum analyzer 

2. Wi de-band frequency discriminator 

3. Filter bank 

4. Doppler frequency tracker 

5. Phase-locked receiver. 

The merits, advantages and disadvantages are discussed by Rolfe et al . 
(1968). In this system principal use was made of a spectrum analyzer 
and to a limited extent of a frequency tracker. 

Sgeetrm Analyzer - The Hewlett-Packard Model 8553B/8552A spectrum 
analyzer used in this investigation is a swept superhetrodyne receiver. 

A simplified block diagram is shown in Figure 5. Essentially the signal 
frequency is compared with a harmonic of the local oscillator frequency 
and the analyzer displays the signal directly in the frequency domain 
as a carrier with its side bands. The center frequency is tuneable, 
and a scan of the total band is selectable. The spectrum analyzer 
resolution is determined by a selectable IF bandwidth. The scan time 
can vary from 1 millisecond to 100 seconds for the selected scan width. 



Resolution (ft) 


LO 
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Figure 5. Spectrum analyzer block diagram. 


If time intervals are too small, power output of the signal may be 
too small to measure. On the other hand, for large time intervals the 
output reflects the spectrum of particle speeds passing through the reso- 
lution focal volume of the beam, and can give therefore only a spectrum 
of velocities (Doppler frequencies) and not an "instantaneous" velocity 
as a function of time. Clearly, for "instantaneous" velocities 
the time interval should be consistent with the focal resolution volume, 
convected particle speed and S/N ratio of the spectrum analyzer. 

In order to convert spectral information in frequency space to 
velocity, use is made of the linear variation of velocity with Doppler 
frequency shift. The frequency bandwidth of the spectrum analyzer is 
"swept" at a rate consistent with resolution of the analyzer and the 
power contained in the bandwidth is recorded on a conventional FM re- 
corder in time space. Conversion from time to frequency hence to velocity 
in principle is simple, requiring only a reference zero frequency and 
known bandwidth or alternatively a calibrated external frequency. The 
rate at which the spectral bandwidth is swept is controlled externally 
to the spectrum analyzer. A schematic arrangement of the process inclu- 
ding preconditioning of the detector signal is shown in Figure 6. A 
typical time, frequency trace of the power output is shown in Figure 7. 






Bias Supply Recorder 


*— * 

ro 


Figure 6. Block diagram of signal detection circuitry. 
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Figure 7. Typical spectrum analyzer output for 
calibration and Doppler frequencies. 


floppier frequency, tracker - A device which provides an output 
voltage proportional to a given Doppler frequency is termed a Doppler 
frequency tracker, or simply frequency tracker. The technique is also 
known as "frequency compressive feedback" or "frequency-locked loop" 

[cf. Rolfe et al . (1968)]. The Doppler frequency, f D ( t ) , is heterodyned 
with a local oscillator frequency. The local oscillator frequency, f LQ , 
is varied so that the difference f^Q - fg is constant and equal to the 
center frequency of a discriminator. The driving voltage of the local 
oscillator is then proportional to f D> hence to the velocity. A sche- 
matic representation of the tracker is shown in Figure 8. 
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E Out 

Figure 8. Block diagram of frequency tracker. 


TEST FACILITIES 

The field site for the experiments was selected at the Colorado 
State University airport (Christman field) located approximately three 
miles west of the city of Fort Collins, Colorado (see Figure 9). The 
test site has a clear field from northwest to northeast, and from south 
to southwest. There are buildings and trees in the range from south to 
east, although the nearest building is some 1100 feet away. To the west 
is the foothills of the Rocky Mountains about a mile distant. The site 
was selected on the basis of land and power availability and proximity to 
the research center about 1/2 mile away. The dominant wind directions in 
the area are north-south, as evident from the alignment of the runway, 
although strong winds also blow over the foothills directly from the west. 




Runway 


Boundary of Emergency 
Runway — . 


rSite of 
\Towers 



o 



v r 

Hangers-"^ 


Figure 9. Field site at Christman Field. 
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The site facilities included two towers and two trailer vans to 
house the instruments and the LDV system. The arrangement shown in 
Figure 10 was to provide clear wind fields to the north and south. As 
winds seldom blow from the east, the instrument vans were located so 
as to cause as little interference as possible to the wind field. 

The 60-ft high tower was used to mount the wind profiling anemo- 
meters. The 40-ft tower was used to mount mirrors to direct the laser 
beam and also to mount the comparison instruments, a climet anemometer 
and wind vane, and a hot wire for turbulence measurements. Photographs 
of the established arrangement are shown in Figures 11 and 12. 


INSTRUMENTATION 

The arrangement of the various instruments in the laser instrument 
van is shown in the photograph of Figure 13. The total instrumentation 
for data taking and recording included the following: 

Spectrum analyzer - The function and description of the spectrum 
analyzer was given in a previous section. 

Frequency tracker - This instrument was also discussed in the 
earlier section. 

Wide hand freq-uency generator - A frequency generator of MHz range 
was used to establish a calibration point for the spectrum analyzer. 
Depending upon the prevailing wind speed, a calibration frequency was 
selected near the extreme of the wind speed range and the scan width of 
the spectrum analyzer was selected to contain this calibration frequency. 

Frequency counter - A frequency counter was used to determine the 


calibration frequency. 








18 



Figure 11. Instrument vans at 
test facility. 


Figure 12. Towers at test facility. 

Profile tower is at left. 




Figure 13. Instrument arrangement inside 
laser van. 
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Fi motion genemtor - A stable function generator was used to drive 
the sweep of the spectrum analyzer IF at a rate consistent with the 
spectrum analyzer scan time. A finite sweep time and ’'flyback" is 
involved. A given combination of sweep duration and scan width has its 
optimum IF filter bandwidth. A table of sample rates for various scan 
time settings is given in Table 1, and bandwidths as a function of scan 
width and scan time is given in Table 2. These tables were reproduced 
from the LMSC report No. D162840 describing the operating procedures of 
the LDV system. 

Miypov "position indicator and drive - The upper mirror on the 
40-ft tower had a motor drive to rotate the mirror about its vertical 
axis. This permitted orientation of the laser beam into nominal align- 
ment with the wind direction. The position of the mirror was indicated 
by a 357 degree potentiometer. There were 3 degrees of ambiguity from 
357 degrees to 360 (zero) degrees. The position pot of the mirror was 
oriented so that zero was due east. 

Climet wind translator - The translator presented wind direction 
and speed as sensed by the cup anemometer and wind direction sensor 
into recordable analog signals. The wind direction sensor was oriented 
so that zero output coincided with due east. The analog signals were 
then monitored on a dual channel strip chart recorder. 

FM tape recorders - Two 14 channel FM tape recorders were used to 
record the analog signals, one a CP-100 Ampex unit and the second an 
FR-1300 Ampex recorder. 

Temperature sensor - A standard bridge and amplifier circuity was 
constructed for this study to measure the deviations in temperature of 
the various thermistors from a reference unit. 
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TABLE 1. MAXIMUM SAMPLE RATES FOR SELECTED SCAN TIMES 


Spectrum Analyzer 
Scan Time 
(Millisec/cm) 

Maximum Sample 
Rate 
(Hz) 

External 
SYNC Period 
(sec) 

0.5 

165 

0.006 

1.0 

69 

0.0145 

2.0 

40 

0.025 

5.0 

18.2 

0.055 

10.0 

5.0 

0.200 

20.0 

3.3 

0.303 


TABLE 2. MINIMUM BANDWIDTHS IN kHz FOR COMBINATIONS 
OF SCANWIDTH AND SCAN TIME 


Scan 

Width/ cm 


Scan 

Time, Millisec/Division 



1.0 

2.0 

5.0 

10.0 

20.0 

50.0 

0.02 

kHz 

0.3 

0.3 

0.1 

0.1 

0.1 

0.1 

0.05 

kHz 



0.3 




0.1 

kHz 

1.0 



0.3 



0.2 

kHz 


1.0 



0.3 


0.5 

kHz 

3.0 


1.0 



0.3 

1.0 

kHz 


3.0 


1.0 



2.0 

kHz 



3.0 


1.0 


5.0 

kHz 

10.0 



3.0 


1.0 

10.0 

kHz 


10.0 



3.0 


20.0 

kHz 

30.0 


10.0 



3.0 

0.05 

MHz 


30.0 


10.0 



0.1 

MHz 

100.0 


30.0 


10.0 


0.2 

MHz 


100.0 


30.0 


10.0 

0.5 

MHz 

300.0 


100.0 


30.0 


1.0 

MHz 


300.0 


100.0 


30.0 

2.0 

MHz 

— 


300.0 


100.0 


5.0 

MHz 

--- 

— 


300.0 


100.0 

10.0 

MHz 

— 

— 



300.0 
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Hot-wire anemome tev - A constant temperature hot-wire anemometer 
was used to measure the atmospheric turbulence. A 100-ft long cable 
was used for the probe and a cable capacitance compensator was used 
for the long-length cable. The hot wires were calibrated with the 
extra cable and compensator. 

Time code generator - A time code generator in I RIG B format was 
used to synchronize the two tape recorders. Usually the times were 
synchronized with the National Bureau of Standards time broadcasts. 

RECORDING OF TEST DATA 

There were in all 26 separate pieces of continuous information 
desired for each test. Two analog 14 channel FM recorders were needed. 
However, two recorders were not available for all tests and some infor- 
mation was therefore sacrificed. The sample data recording sheet shown 
in Figure 14 indicates the data recorded on each channel of each 
recorder. They were arranged in such a way that temperature and 
humidity data were sacrificed when the second recorder was unavailable. 

The data can be grouped into the following sets. On the 60-ft 
tower, six levels of wind speed were obtained to establish the vertical 
profile of the wind field in which comparison data were taken. These 
were grouped in the CP-100 Ampex recorder. Also, on the same tower, 
there were six levels of temperature measurements to determine the 
temperature profile and four levels of wet bulb temperatures to establish 
the humidity profile. These were grouped on the FR-1300 Ampex recorder. 
On the 40- ft tower the comparison instruments, the cup anemometer, the 
wind vane, and the hot wire were mounted. These data together with the 
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Q NASA"MSFC Field Test Site, Huntsville, Alabama 

□ Airport Field Test Site, Foothills Campus, Colorado State Univ», Ft. Collins 

□ Other 

Test Conducted Between a,m,/p,m. and a. m. /p. m. on 


METEOROLOGICAL DATA 

Air Pollution Index: . Visibility: QGood; GFair; Q poor 

Sky Condition: □ Clear; □ Light Clouds; □ Medium Clouds; QHeavy Overcast 

Temperature °F ; Relative Humidity % or Dew Point °F; 

Barometric Pressure mb; Anemomete r(s) Locale 


Time into Test (min) 

— 

0 

15 

30 

45 

60 


Mean Wind Speed (knots, mph, ft/sec) 







Mean Wind Direction (deg wrt north) 







Laser Coolant Temperature (°F) 








General Weather Conditions {frontal presence, rain in past 12 hours, etc.); 


OPTICAL CONFIGURATION 

Mirror Orientation deg (wrt north) 

Telescope mirror to lower tower mirror distance: ft in. 

Distance between lop and bottom mirrors on tower* ft in* 

Total distance from telescope mirror to focus vol: ft 

Homodyne configuration: □ Mach Zehnder; □ internal cavity 

Laser power into telescope: watts; Power at focal volume; watts; He dewar check □ . 

Telescope mirror size; in. diam.; Lens focal length: in; Detector type: _____ - 

Comments: 


SPECTRUM ANAEYZER/AVERACER DATA 

Sweep Rate: ms/cm; Sample Rate: samplea/aec. 

Number of sweeps averaged per sample; _____ . 

Frequency dispersion: MHz/cm. Filter bandwidth: kHz; Bandwidth: kHz* 

Othe r: 


F* RECORDER DATA - MODEL CP-100 AMP EX 
Label Tape Reel with Test No. 

Tape No. ; Tape Speed lps; Response Hz. 


Cha. ^o. 

Contents ! 

1 

Voice: Test ident. etc. 

2 

Spcct, analy. sync, pulse: 

3 

Spect. analy. Y out. FreQ. disp. 

4 

Kind dir. Climet volts — > ^wrt 

5 

Mind speed Cl met volts — > fps 

6 

Mirror azimuth volts — > wrt 

7 

Hot wire guvesoeeter 

8 

Fixed tower dat3 - Nincf sp. level 1 

9 

Wind sp. level 2 

10 

wind sp. level 3 

11 

Kind sp. level 4 

12 

Wind sp. level 5 


is i nd sp. level o 

1 14 

Tine code ident. 1R1C b j 


FM RECORDER DATA - MODLL FR-1300 AMP EX 

Label Tape Reel with lest so. 

Tape No. ; Tape Speedy ips; Response Hz. 


Chj. No. 

Content. 

F i 

xcd lower Data) 

1 

I herrt) 

j. ' ct\ 

d. ib . 

tc.:.p. level J 

2 'Therm 

. d : f , . 

d iv 

•iir level 3 

3 

Therm 

, Jiff. 

Jr> 

air level 1 

4 

Then:! 


lit : 

Du Jo level 1 

S 

Therm 

, Jiff. 

Jr^ 

air level 4 1 

6 

Therm 

, Jiff, 

“C L 

bulb level 4 

7 

Therm 

, Jiff 

Jr> 

a j r 1 v ve 1 5 

8 

TWmi 

, J i ff. 

M Ot 

Ufl b > 

9 

Therm 

, Jiff. 

dr> 

a:r level o 

10 

Therm 

, dift. 

wet 

l 'ulh level b 

11 

hind i 

direct. on L 

u'! > • > 

12 



Li . 

14 

Time 

:ode ide 

nt . 

IK-Iti B : 

Aux 

Voice 

Test 

.etc. 1 


AUDIO RECORD 


OTeat Identification Number; Q Spectrum Analyzer settings (sweep rate, number samples averaged, 
etc.); OMean Wind; Q Distance from telescope mirror to focal volume; Q Visual quality of signal, 
and O Problems and other comments. 


(Signed) Test Engineer 


Figure 14. Sample data sheet. 
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spectrum analyzer signal and appurtenant data were grouped into the 
CP-100 recorder. On one channel of each recorder there was an IRIG B 
time code for referencing the two sets of data to corresponding times. 

A voice channel (direct record) was reserved for verbal description of 
conditions and problems which occurred during a test. 

Data with a frequency tracker were taken during a period when the 
second tape recorder was unavailable. Since two additional channels 
were required to record the signals from the tracker, two levels of wind 
speed data were sacrificed on the CP-100. These were levels 2 and 4. 

TEST PROCEDURE 


Pre-Test Preparation 

Preparations for recording one-hour of continuous wind data and 
associated documentation was elaborate and time-consuming. For any 
given test, or attempted test, the following routine was necessary. 

Cooling the ge-flg deteator - The Ge-Hg detector was pre-cool ed 
with liquid nitrogen for a period of about one hour before filling 
with liquid helium. This procedure was followed primarily to conserve 
liquid helium, which is comparatively many times more expensive than 
liquid nitrogen. Just prior to data-taking, after all preparations 
were completed, the dewar of the detector was filled with liquid helium. 

Ogtios alignment - Before each test, alignment of the optics was 
necessary. A specific alignment procedure progressing outward from the 
laser to the tower was necessary. Although the beam splitters and 
mirrors did not require frequent adjustment, the optical beam on which 
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the focusing mirror was mounted required frequent adjustment. As the 
scattered radiation was redirected back into the laser, slight mis- 
alignment of the optical axis caused poor to no heterodyning, hence weak 
or no Doppler detection. Since alignment of the focusing lense mount is 
coupled with the diagonal and the schlieren mirror, a sequence of trial 
and readjustment was usually necessary. 

After the optical beam was adjusted, the diagonal required minute 
adjustment to center the diverging radiation on the schlieren mirror. 

The schlieren mirror in turn required adjustment to direct the converging 
beam through the end of the 9-ft long tube. Thereafter, the entire 
mounting table required movement to center the beam on the lower external 
mirror near the base of the 40-ft tower. If the optics were bumped out 
of alignment during this process, then the entire procedure was restarted. 

Once the laser beam was centered on the lower mirror, then the lower 
mirror was adjusted to center the beam on the upper one, and finally the 
upper mirror was rotated to direct the beam as closely as possible either 
directly into the prevailing wind direction or downwind along the wind 
direction, checking also to see that the beam was parallel with the 
ground. To establish the latter, an identification mark on the adjacent 
60-ft tower was used to place the line of sight parallel to the ground, 
hence the axis of the laser beam was in the horizontal plane of the mean 
wind. 

Profile tower - The thermistors on the 60-ft tower were arranged in 
a radiation shield, with a suction pump arranged to draw 2 ft/sec air 
velocity over the "dry bulb" thermistor and 30 ft/sec over the "wet bulb" 
thermistors. Distilled water was forced up the tower by air pressure 
into water wells with wicks leading to the "wet bulb" thermistors. 
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These thermistors were checked before each test and wicks were prewetted 
to insure that the distilled water would be drawn up from the wells. 

Hot-wire anemometer - The hot-wire anemometer which was dismounted 
during a non-test period was remounted. The wire was placed in a verti- 
cal axis and the probe was oriented toward the wind and in a location 
such that there was no interference from the mirror, the cup anemometer 
or the tower itself. 

Pre-Test Calibration 

Ta£e recorder - The FM record amplifiers of the tape recorder are 
subject to slight deviations from calibrated conditions from day to day. 
To account for these deviations, a five-level DC signal was provided as 
a calibration of tape-recorded (and playback) voltage against a "true" 
voltage registered by a calibrated digital voltmeter (DVM). Since in 
the data set, a continuous square-wave signal was recorded, the calibra- 
tion set did not include a sinusoidal signal of known rms value. 

Climet anemometers - Both climet anemometer translators were 
calibrated for zero and full scale 1 volt outputs and recorded on the 
tape recorder. Prior to mounting the anemometers in the towers, all 
cup anemometers as well as the hot wires were calibrated in the Colorado 
State University wind tunnel against a pitot probe of known performance. 
Calibrations were performed twice, in February and June 1971. 

Mirror position - The mirror position, with zero oriented directly 
east for convenience, was calibrated for zero and full scale ouput, with 
the assumption of linearity with angular position. Since the position 
was indicated by a potentiometer, the assumption seems justified. 

Spectrum analyzer - Proper settings of the spectrum analyzer 
controls were established consistent with the prevailing wind speeds. 
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The sweep of the spectrum analyzer was triggered by an external square 
wave from a stable function generator by a change from negative to posi- 
tive voltage. A known deviation frequency was input to the spectrum 
analyzer and the resultant signal from the IF output was recorded as the 
frequency band was swept. This calibration thus provided the reference 
for determining velocity from the Doppler shifted frequency of the back- 
scattered radiation. 

Noise gaMbvation - The final pre-test calibration was made of the 
background noise emitted from the detector. With the detector dewar 
charged with liquid helium, and the main laser beam to the telescope 
blocked, the output signal from the detector which consisted only of 
noise was recorded. 

Data Recording 

After completing the pre-test preparations and calibrations, data 
were recorded on the tape recorders for nominal periods of one hour 
duration. Constant monitoring of the data was provided, and instrument 
adjustments when necessary were properly recorded as to time and nature. 

The turbulence range of interest extended only to a maximum of 5 Hz, 
thus the CP-100 tape recorder was operated at inches per second (ips) 
and the FR-1300 recorder at 1-7/8 ips. The higher speed of the CP-100 
recorder was necessary to record the Doppler signals from the spectrum 
analyzer. At 7h i ps the recorder amplifiers are responsive to 2.5 KHz. 

Anomalies in the data noticed were recorded on a voice channel 
(direct record) of the tape recorder, as well as on the data record 
sheet (Figure 14). 
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DATA REDUCTION PROCEDURE 

All data for this investigation were analyzed digitally, the 
digitizing being done in prescribed sets in simultaneous sample and 
hold mode at the NASA-MSFC computer center. The digitized data were 
analyzed at the Colorado State University computer center. 

Selection of Digitizing Rates 

The turbulent frequencies of interest in this study are less than 
5 Hz, thus the digitizing frequency should be at 10 samples per second, 
and also, because in general the recorded information should be related 
to the same instants of time, a simultaneous sample and hold mode was 
used in digitizing. The analog signals were filtered at 5 Hz (real- 
time base). 

The scan rate of the spectrum analyzer for Doppler frequencies was 
16 Hz. Since the Nyquist frequency is equal to one-half sampling 
frequency, 

f - f ° 
f N ' T~ 

the highest frequency information contained in the recorded signal is 
8 Hz. However, the usual criterion of digitizing rate to obtain this 
frequency information does not apply. The objective in data reduction 
was to determine the location (time base) of the Doppler signal with 
reference to zero frequency, hence of Doppler frequency and of wind 
velocity. The bandwidth and resolution of the spectrum analyzer deter- 
mines the nature of the Doppler signal. If we view the peak signal in 
the bandwidth as depicting the mean velocity in the prescribed resolu- 
tion interval, then the digitizing rate of the Doppler signal is 
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independent of the spectrum analyzer settings. Thus with a view to 
maximizing the frequency resolution (of the peak) in a given sweep, 
a choice of 250 points per sweep was made. The choice of this digiti- 
zing rate does however affect the total quantity of digitized data. Two 
channels of information, the external function generator and the IF 
output of the spectrum analyzer, were digitized at this higher rate, 
multiplexed on digital magnetic tape in binary format. The sampling 
rate for these channels was thus 4 KHz/channel and the data were filtered 
at 2 KHz. 

Multiplexed Data Groups 

The 26 channels of analog information were digitized in three 
separate groups. 

Group 1 - The sweep signal (square wave) and the spectrum analyzer 
IF(y) output were multiplexed and digitized at a rate of 4 KHz/channel. 

Group 2 - The climet anemometer and wind direction sensor, the 
mirror position, the hot wire output and six levels of wind speeds on 
the profile tower were multiplexed and digitized at a rate of 
10 Hz/channel . 

Group 3 - The ten channels of thermistor data were multiplexed and 
digitized at a rate of 10 Hz/channel. 

The remaining four channels of voice, time codes and wind direction 
on the profile tower were not digitized and were retained for reference. 
The time code information was of course used to identify the regions of 
the analog tape which were digitized. 
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Data Format 

The A/D converter used at NASA/MSFC generated data words of 10 bits 
plus sign. The packed format of the multiplexed data therefore were 
written in groups of 11 bits. The CDC 6400 at Colorado State University 
is a 60-bit word machine, thus some tape reading problems were presented 
with the original format of the generated data tape. In order to reduce 
the reading problem, the original data tapes were reformated to give data 
words which were 11 bits plus sign, or 12 bit words where a zero was 
inserted into the most significant bit. The packed 12-bit data words 
were thus conveniently separated and sorted from the 60-bit computer 
word. 

The data included a record of header information at the beginning 
of each data set, and a 24-bit time word at the beginning of each data 
record. This time word is a reference digitizing time, and relates to 
real time in accordance with the ratio of record to playback tape speed. 
However, for records of the order of 60 minutes real time duration, the 
time word (expressed in milliseconds) becomes excessively large. Thus 
the digitizing clock which recycles after 100 seconds requires accounting 
of the cycles to convert digitizing time to real time as well as the 
ratio of record to playback speeds. 

Data Reduction 

Laaev_ Doggier signals - The bulk of data reduction involved the 
Doppler signals. The view adopted in computer program formulation was 
to devise a general, automatic program. This was successful to a 
degree, however sufficient problems with data anomalies were encountered 
that some initialization is necessary. Considerable time was spent in 
developing this feature of a data reduction program. In retrospect, 
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perhaps less automatic, sequential programs would be more economical in 
terms of total effort. The flow chart for the program is shown in 
Figure 15 and a listing is given in Appendix A. 

The essential technique is as follows: Data from Group 1 (identi- 

fied above), and the first channel of the multiplexed data of Group 2, 
are necessary to convert the spectrum analyzer data to wind speed. If 
the mirror direction varies in the data period, that information is also 
required. 

The cup anemometer wind speed, the hot wire data and the profile 
information can be processed separately, but because the two groups of 
data were arranged on different tapes and had to be read in "simulta- 
neously" to analyze the Doppler signal, the program included processing 
of these data at the same time. It should be noted here that several 
alternative methods were recognized from the outset, and a one-pass 
automatic program seemed feasible and most desirable. Ultimately the 
profile data program was separated from the others and analyzed in a 
separate pass. The flow chart in Figure 15 reflects this variation to 
the original technique. 

The program first determines the input-output calibration of DC 
voltage. This calibration enables conversion of such data as wind and 
mirror directions, cup anemometer speeds and hot wire turbulence velo- 
cities from tape output voltage to true voltage hence to the physical 
quantities. The next step in the analysis is to determine the calibra- 
tion Doppler frequency. That is, the known frequency input is identi- 
fied in the time space (number of points) from zero frequency, and 
since velocity is linear with Doppler frequency, then calibration is 
obtained for the velocity component along the laser axis. In order to 
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Figure 15. Simplified Flow Chart of Doppler Data Reduction. 
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distinguish the Doppler "peak" from the background noise, the noise 
calibration established the noise level across the entire frequency 
band of the spectrum analyzer. In the program the S/N ratio is a 
variable and may be set at any level compatible with the recorded 
Doppler signal. 

The first step in the data analysis is to read in one record from 
the multiplexed Group 2 data. Each digital value is converted to velo- 
city, and reference times for each value are calculated. The velocities 
and reference times are stored. The cup and hot wire data are digitized 
at identical times, thus one reference time serves both channels of 
information. Means and variances are calculated. Wind direction voltages 
are averaged for 10 seconds (one record) and converted to angle with 
respect to the laser beam. The value is temporarily stored. The mirror 
azimuth (direction) is averaged and checked. If no change occurred 
(i.e. the mirror was not rotated) the information is redundant and 
discarded. 

The first record of Group 1 is then read in. Each spectrum analyzer 
scan, approximately 250 points, is searched for zero frequency (the change 
in voltage of the square wave from negative to positive identifies the 
beginning of the sweep) and the Doppler signal. The reference time for 
Doppler-converted velocity is referenced to the beginning of the sweep. 
Successive sweeps and time words at the beginning of each record references 
the true time of the calculated Doppler-measured velocity. The first 
identifiable Doppler peak is accepted as the measured velocity. To deter- 
mine the peak value, comparison is made to successive points, and if the 
signal level (voltage) drops, the previous point is accepted as the 
Doppler frequency. It is possible that in a given sweep there is no 
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Doppler signal {signal dropout), in that event, the velocity determined 
in the previous sweep is recorded. The Doppler-indicated velocity is 
then converted to wind velocity by the 10-second average angle of the 
wind direction with respect to the laser beam axis (mirror direction). 

There are 3000 data points (2 channels) in each record of the 
Group 1 data. This corresponds to 6 sweeps of the spectrum analyzer 
and 0.375 second in terms of real time. Successive records of Group 
1 are read in and analyzed until the real time reference period 
exceeds the real time period of the data read in from the Group 2 
data. Additional Group 2 data are then read and reduced, and the 
process repeated. 

The stored values of velocities and reference times are periodically 
purged from storage and written on a magnetic tape. Thus the entire test 
record is converted to velocity-time history with the same reference 
times for the cup anemometer and hot-wire data, but a different reference 
time for the Doppler-indicated velocities. 

The generated velocity-time history tape is then reprocessed to 
obtain the statistical characteristics of the turbulent wind data. These 
characteristics are the mean, variance (standard deviation), probability 
density and spectral densities (power spectrum). 

Veloaitu profiles - The velocity profiles are calculated in a 
straight forward manner, using the other six channels of data in Group 
2. Only the mean values are of concern, and ten-minute average velo- 
cities are calculated for each anemometer. The calibration data for 
voltages, and the prior wind-tunnel calibrations, are all that 
are required. A program listing is given in Appendix A. 
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Temperature profiles - Temperature and humidity profiles likewise, 
are relatively straightforward requiring manufacturer's calibration data 
for the thermistors and conversion of average tape voltage to true 
voltage. The resistances are calculated from a standard bridge equation, 
hence temperatures are determined. The program listing is given in 
Appendix A. 


EXPERIMENTAL RESULTS AND DISCUSSION 


Calibrations 

Climet anemometers - Calibration curves of the climet anemometer. 
Series No. 828, are shown in Figure 16. The calibration was performed 
in a wind tunnel with the translator set for 1 volt output at 1896 Hz 
input (signal frequency generated by the cup) for the 60 scale setting 
on the translator. Ordinarily, the translator is adjusted to output 1 
volt for specific input frequencies on each scale. However, for purpose 
of this calibration, adjustment was made for 1 volt output on the 60 scale 
only (any frequency would have served as well) and outputs read from both 
30 and 60 scales. In setting the translator during an experiment, there- 
fore, adjustment was always made only for the 60 scale. The output is 
linear with velocity as seen in the figure. 

The CP-100 tape recorder has a low input impedance, causing a 
loading problem with virtually all the instruments connected to it. 

Thus the cup anemometers and hot wires were calibrated with the outputs 
connected to the tape recorder. 

Hot-wire anemometer - A typical calibration curve for the hot-wire 
anemometer is shown in Figure 17. For purpose of this investigation, 
the King's law relationship is shown, and it is seen that in the region 
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Output in Volts 


Figure 16. Calibration curves for climet anemometer. 



36 



E 2 (Volt ) 2 


Figure 17. Hot-wire calibration curve 
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of interest, the curve was linear. A linearizer was not used with the 
anemometer. Instead, each digitized data point was converted to actual 
voltage and velocity calculated from the calibration. 

Measurements of Run 5Q801 (May 8, 1971) 

The data for this test were taken from 1:48 pm to 2:45 pm, covering 

a period of approximately one hour. At the beginning of the test the 
wind was blowing from the south-southeast (30 degrees east from south) 
which gradually changed to south-southwest (15 degrees west from south) 
by the end of the test period. The wind speed was reasonably constant 
at about 4 m/sec (9 mph) throughout the test period. Particle counts in 
the atmosphere were not available for this test; however, with the pre- 
vailing south wind, the pollution from Denver was evident as a blue haze 
along the horizon. This was also reflected in the strength of the Doppler 
signals on the spectrum analyzer. 

Velocity - The velocity profiles for successive 10-minute 

periods throughout the test are shown on Figure 18. The velocity profiles 
were logarithmic as expected; however, the slope of the profiles differ, 

indicating that the effects of accelerating and decelerating winds 
(gusts) are reflected in the profiles. It will be seen in the time 
traces of velocities that the fluctuations are of the same order of 
magnitude as the means, and the mean values change with time. The 
analysis to establish the profiles assumes piece-wise stationarity. 

Spectrum analyser settings - The following settings were made on 
the spectrum analyzer: 

Sweep rate: 5 ms/cm 

Sample rate: 16 sweeps/sec 
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Figure 18. Velocity profiles for test period 50801. 
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Frequency Dispersion: 0.2 MHz/cm 

Filter Bandwidth: 10 KHz 

Bandwidth: 30 KHz 

The calibration frequency was 1.007 MHz (5.34 m/sec) which is 
pictured in Figure 19. The noise level from the detector is shown in 
Figure 20. The photograph is the oscilloscope trace from playback (at 
record time) of the recorded signal on the CP-100. The signal is inver- 
ted to avoid confusion with the square wave shown at the top part of the 
picture. The vertical scale is 200 mv/cm. 

Typical Doppler signals are shown in Figures 21 and 22. As noted, 
the S/N ratio is large, but the spectral bandwidth is also large. Peaks 
in the signal of the kind shown in Figure 21 are relatively easy to 
determine; however, multiple peaks are evident in Figure 22. In these 
instances, the first largest peak is detected, and the others ignored. 
There were undoubtedly particles of different sizes in the focal region 
with different angularity with respect to the laser beam axis which cause 
the multiple peaks in a given sweep. 

Velocity time traces - Time traces of velocity from the cup anemo- 
meter, hot wire and the LDV, for two consecutive 4.26-minute periods are 
shown in Figures 23 and 24. Mean velocities for each 4.26-minute 
interval have been subtracted; the fluctuations thus are referenced to 
zero for each plot. 

As seen in these traces, there is reasonable conformance between 
the cup anemometer, hot wire and LDV outputs. It should be noted here 
that the cup anemometer was at a level 11.3 meters above ground, the 
hot wire was 0.3 meters below the cup level and the laser beam axis was 
at the same level as the hot wire although the focal region was 3 meters 
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Figure 19. Calibration frequency 1.007 MHz. Figure 20. Detector noise calibration. 

Test 50801 Vertical scale is 200 mv/ctn. 

Test 50801 



Figure 21. Sample Doppler signal. 
Test 50801 


Figure 22. Sample Doppler signal. 
Test 50801 
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Figure 23. Time traces of wind velocity. 

Test 50801, Interval 1 

(For means and variances see Table 3) 
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Figure 24. Time traces of wind velocity 
Test 50801, Interval 2 
(For means and variances see Table 3) 
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farther upwind. It should be noticed in making visual comparisons that 
the vertical scales are different for the traces. 

Means and variances - The means and variances from a 34-minute 
interval of the total record were analyzed and are shown in Table 3. 

The choice of a 34-minute period was based largely on the limitations 
of the spectral analysis program. This was also a sufficiently large 
period to reflect a reasonable confidence interval for the spectral 
densities. 


TABLE 3. MEANS AND VARIANCES FOR TEST 50801 


4.26-Minute 

Intervals 

Mean 

Velocities m/sec 

Variances (m/sec) 2 


Cup 

Hot Wire 

LDV 

Cup 

Hot Wire 

LDV 

1 

4.203 

4.232 

4.044 

.612 

.604 

.689 

2 

4.486 

4.488 

4.253 

.539 

.524 

.672 

3 

3.762 

3.799 

3.585 

.340 

.258 

.348 

4 

4.245 

4.270 

4.247 

.458 

.355 

.596 

5 

3.976 

4.000 

3.953 

.444 

.340 

.526 

6 

3.823 

3.847 

3.693 

.342 

.342 

.503 

7 

3.618 

3.642 

3.489 

.623 

.598 

.573 

8 

4.212 

4.235 

4.073 

.461 

.361 

.674 

Averages 

4.041 

4.064 

3.917 

.472 

.413 

.567 


The mean wind speeds detected by the LDV is in overall 3 percent 
agreement with the cup anemometer, and within 5 percent for any given 
4.26-minute interval. The greater spread for smaller time intervals is 
to be expected because of the spatial spread of sampling points for the 
two instruments. 
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The variances for LDV are larger than those detected by either 
the hot-wire or the cup anemometer. It is surprising to note also 
that the variances for the hot wire are less than that for both the 
cup anemometer and LDV measurements. The greater variances for the 
LDV results are due in part to the fact that only mean horizontal 
angularity of the particle motion with respect to the laser axis is 
included in the correction. Thus there are greater variations of 
velocities from the mean. This is observed also in comparing the 
mean speeds for the three data sets. The mean is lower for the LDV 
as compared to cup speeds. 

ProbabiMtp distrn^uM^ns - The distributions of velocities about 
the means for the three instruments are shown in Figure 25. These data 
are in terms of standard deviations, and are not normalized so that 
straight lines are drawn from one data point to another. The distribu- 
tions are skewed to the right. This skewness is governed by the nature 
of the turbulence in the atmosphere rather than by instrument response, 
as it can be seen that all three instruments respond similarly. The 
percentage of data near the mean is greater for the cup anemometer than 
for the other instruments, as was suggested in the preceding paragraph, 
the percent of low velocities appear to be greater for the LDV than for 
either cup or hot wire measurements. 

Spectral densities - The spectral densities for measured turbulence 
in the atmosphere are shown in Figures 26, 27 and 28 for the cup, hot 
wire and LDV instruments, respectively, and a comparison of the three 
are shown on Figure 29. 

There are apparent energy concentrations in the spectra for the 
cup anemometer and hot wires at 5 Hz which are also noted at 2.5 and 
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Standard Deviation m/sec 


Figure 25. Distributions of velocities about the mean. 
Test 50801 
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Figure 26. Spectral density distributions for cup anemometer. 
Test 50801 
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Figure 27. Spectral density distributions for hot-wire anemometer. 
Test 50801 
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Figure 28. Spectral density distributions for LDV data. 

Test 50801 



PRECEDING PA^ BLANK NOT FILMED 

1.25 Hz. These must be due to mechanical aliased frequencies from the 
tape recorder, for they appear in the hot-wire and cup anemometer data 
but not in the LDV data. Mechanical aliasing does not appear in the 
LDV data because of the manner in which the velocity-time history is 
generated (see section on data reduction). 

If the aliased spectral densities are ignored, it can be seen that 
the hot-wire and cup anemometer have identical spectra up to 0.4 Hz. 
Beyond that frequency, the spectrum decreases because of the limited 
frequency response of the cup anemometer. The cup anemometer data may 
in principle be corrected by a frequency response function (see Camp 
1965), but in this study the correction was not made, as the comparison 
spectrum for higher frequency is given by the hot-wire anemometer data. 
The response of the constant temperature hot wire used here is up to at 
least 1 KHz and the data were filtered at 5 Hz before digitizing. 

As it is seen on Figure 29, the spectral densities for the LDV- 
measured turbulence is slightly greater for frequencies less than 1 Hz, 
but essentially parallel to the hot-wire data. For higher frequencies, 
there appears to be more energy contained in the LDV-measured turbulence. 
This must be aliased information because the hot-wire data do not show 
this trend. 

The aliasing must arise from the technique used in data reduction. 
While the spectrum analyzer is being swept (sampled) at a rate of 16 Hz, 
thereby effectively establishing the Nyquist frequency, the velocity 
time data cannot be filtered at 8 Hz before the sampling is done. That 
is, turbulence of higher frequency transporting aerosol and solid parti- 
cles in the atmosphere are sensed in the resolution volume of the LDV. 
Thus in calculating the velocity from the sampled spectrum, the aliasing 
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from higher frequency cannot be avoided. What is surprising, however, 
is to note the magnitude of the aliased spectrum in the LDV-measured 
turbulence indicated by the deviation beyond 1 Hz. 

Measurements of Run 32701 (March 27, 1971) 

The data for this test were taken from 3:30 pm to 4:18 pm, a period 
of 48 minutes. The wind was essentially steady from the north-east (60 
degrees east from north) at around 12 m/sec (27 mph). Particle counts 
in the atmosphere were not available for this test. There was an arctic 
front moving in from the north and the air was "clean." Visibility was 
virtually unlimited. The laser beam axis was directed downwind in this 
test because the direction of the wind was such that the laser beam axis 
would have been close to a vertical leg of the tower. 

Velocity profiles - The velocity profiles for successive 10-minute 
intervals are shown on Figure 30. The profiles are logarithmic and the 
mean velocities increased in the first 20 minutes of the 50-minute period 
and decreased thereafter. The spread of mean velocities for the total 
period varied from about 10.7 to 13 m/sec at the level of the focal re- 
gion of the laser beam. 

Spectrum analyser settings - The settings of the spectrum analyzer 
were as follows: 

Sweep rate: 5 ms/cm 

Sample rate: 16 sweeps/sec 

Frequency dispersion: 0.5 MHz/ cm 
Filter Bandwidth: off 


Bandwidth: 


30 KHz 



Height above Ground in meters 
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Figure 30. Velocity profiles. Test 32701 
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The calibration frequency was 4.009 MHz (21.2 m/sec) which is shown 
in Figure 31. The noise level from the detector is shown in Figure 32. 
The vertical scale in the oscilloscope trace is 100 mv/cm. 

Typical Doppler signals are shown in Figures 33 and 34. As noted, 
the S/N ratio is small and the spectral dispersion is also small. There 
were larger periods of signal dropout, that is sweeps when there were 
no detectable signals. In these instances the analysis was made assuming 
that the velocity indicated in the current sweep was equal to that of the 
previously detected velocity. 

yetooit.y time ti\zee_s - Time traces of velocity from the three 
instruments are shown in Figures 35 and 36 for two representative 
4-minute time intervals. 

There is reasonable agreement between the cup anemometer and hot- 
wire traces in general trend of mean velocities. However, the turbulent 
fluctuations in the hot-wire signals are greater than that indicated by 
the cup anemometer traces. The LDV signals have several peculiarities. 
The fluctuations are clipped at both the upper and lower limits. These 
clipped signals are results of the low S/N ratio and the computer pro- 
gram. As indicated previously, the low particle concentration in the 
atmosphere often caused no detectable signal in a given sweep of the 
spectrum analyzer. In such instances the velocity was set equal to the 
immediately- previous calculated velocity. At the lower end, the signal 
was lost in the noise (see the noise calibration trace of the oscillo- 
scope) and a previously higher value was then identified as the velocity 
for that sweep. There are noticeable high peaks in the LDV trace. It 
is believed that these signals are spurious, resulting from identifica- 
tion of high noise peaks as Doppler signals. The trend of mean 
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Figure 31. Calibration frequency 4.009 MHz. Figure 32. Noise calibration. 

Test 32701 Vertical scale is 100 mv/cm 

Test 32701 



Figure 33. Typical Doppler signal. 
Test 32701 


Figure 34. Typical Doppler signal. 
Test 32701 
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Cup Anemometer 


Hot Wire 


LDV 


Figure 36. Time traces of wind velocity 
Test 32701, Interval 5 
(For means and variances see Table 4) 
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velocities is generally identifiable, but the comparison is not as 
favorable as for test 50801. 

Means and variances - The means and variances from a 34-minute 
interval of the total record are given in Table 4. 


TABLE 4. MEANS AND VARIANCES FOR TEST 32701 


4.26-Minute 

Intervals 

Mean Velocities m/sec 

Variances (m/sec) 2 

Cup 

Hot Wire 

LDV 

Cup 

Hot Wire 

LDV 

1 

12.041 

12.152 

11.697 

2.686 

5.067 

4.326 

2 

13.659 

13.835 

13.460 

4.951 

4.281 

6.111 

3 

13.164 

13.203 

13,990 

6.497 

7.258 

9.897 

4 

13.973 

14.094 

14.226 

4.117 

5.382 

5.415 

5 

13.486 

13.575 

14.557 

5.167 

5.429 

6.833 

6 

12.658 

12.697 

12.441 

2.812 

3.349 

6.620 

7 

12.417 

12.578 

11.570 

4.000 

6.290 

3.193 

8 

11.453 

11.551 

10.071 

2.934 

3.826 

2.802 

Averages 

12.856 

12.961 

12.751 

4.093 

5.040 

5.448 


The average wind speed detected by the LDV in the 34-minute period 
is within 1 percent of the cup and hot wire averages. There are larger 
variations however for the shorter 4.26-minute intervals, and as the 
time traces would suggest, variations become greater for even shorter 
periods. As noted in the preceding section, these are undoubtedly caused 
by the spurious signals in the velocity calculations. The mean velocities 
measured by the hot wire were generally larger than the cup anemometer, 
and the variances as expected are definitely greater because the frequency 
response of the cup anemometer is limited. 



Over a 34-minute period, the fluctuations (variances) detected by 
the LDV are larger than those of the hot wire. This was also true for 
Test 50801 which had considerably lower mean wind speeds. Again, the 
spurious signals in the LDV velocities contribute significantly to 
variances. 

Probabititii i iietrffcMttons - The distributions of velocities about 
the means for the three instruments are shown in Figure 37. Turbulence 
velocities are skewed to the left for all three instruments. The LDV 
data indicated difficulty in tracing the larger velocities. As explained 
previously, this could be due in part to the three dimensional nature of 
turbulence and only the horizontal angularity was corrected (in the mean) 
in these measurements. This feature of the LDV traces was noted also for 
test 50801. 

SgeatvaX densities - The spectra for the cup anemometer, hot wire 
and LDV data are shown in Figures 38, 39 and 40, respectively. For 
comparison, the three are replotted in Figure 41. Spikes of high fre- 
quency are again noted at 2.5 and 5 Hz in the cup anemometer spectra. 

It was noted that the time traces of the LDV data included spurious 
spikes of high velocity. These spikes are transformed into the spectra 
and are noted particularly as spikes of power near 1 and 3 Hz. These 
spikes in the spectra were ignored in replotting on Figure 41. 

The spectra of turbulence measured by the LDV and hot wire compare 
favorably. This is also indicated by the comparison of variances in 
Table 4. The cup spectra however drops off at around 0.2 Hz because 
of the limited frequency response. Response corrections for the cup 
anemometer were not made. 
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Figure 37. Distributions of velocities about the mean. 
Test 32701 



Spectral Density {m/sec) 2 /H 



Frequency in Hz 

Figure 38. Spectral density distributions for cup anemometer data. 
Test 32701 
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Figure 39. Spectral density distributions for hot-wire anemometer. 
Test 32701 
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Figure 40. Spectral density distributions for LDV data. 
Test 32701 



Spectral Density (m/sec) 2 /Hl 




64 


Measurements of Run 101401 (October 14, 1971) 

The data for this test were taken from 9:16 pm to 9:55 pm, a 
period of 39 minutes. The wind was from the north -northwest across the 
clear grassland. The mean wind speed varied from about 4 m/sec at the 
start of the test to about 5.7 m/sec at the end. The wind direction 
remained constant. With a northern weather front moving in, the air 
was clear, (little pollution), and visibility was good. 

S^ectrwv analyzer setting - The settings of the spectrum analyzer 
were as follows: 

Sweep rate: 5 ms /cm 

Sample rate: 16 sweeps/sec 

Frequency Dispersion: 0.2 MHz/cm 

Filter Bandwidth: 10 KHz 

Bandwidth: 30 KHz 

The calibration frequency was 1.691 MHz, which is shown in Figure 42. 
The noise level is shown in Figure 43. It will be noted that reference 
zero frequency is shifted slightly from the pulse rise of the square 
wave, resulting from a horizontal axis shift of the spectrum analyzer. 

An accounting of this shift was made in data analysis. 

A sample trace of one sweep of the spectrum analyzer is depicted 
in Figure 44. The S/N of the Doppler trace is small but was sufficient 
to discriminate from noise. There were drop outs in Doppler signature 
as indicated by the time traces of wind speeds. 

Velocity time trapes - Time traces of wind speeds from the cup 
and hot wire anemometers and the LDV are shown for representative 
4-minute intervals in Figures 45 and 46. As with the two previous tests, 
the mean trends correspond with apparent differences in turbulence 
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Figure 42. Calibration frequency 1.691 MHz. 
Test 101401 


‘Figure 43. Noise Calibration. 
Test 101401 



Figure 44. Sample Doppler signal. 
Test 101401 







66 




Figure 45. Time traces of wind velocity. 

Test 101401, Interval 1 

(For means and variances see Table 5) 



Velocity m/sec Velocity m/sec Velocity m/sec 
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Hot Wire 



Time -Minutes 


Figure 46. Time traces of wind velocity. 

Test 101401, Interval 5 

(For means and variances see Table 5) 
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fluctuations. The large number of low points in the LDV signature 
resulted from the low S/N ratio; particularly by having to set a low 
level trigger in the computer program. The spurrious high peaks are 
believed to be caused by extraneous signal in the Doppler sweep. There 
are not enough of these to cause difficulty with the statistical anslysis. 

Means and variances - Means and variances for the entire 34-minute 
test period are given in Table 5 for each 4. 26 -minute segment. 


TABLE 5. MEANS AND VARIANCES FOR TEST 101401 


4.26-Minute 

Intervals 

Mean 

Velocities m/sec 

Variances Jm/ sec) 

2 

Cup 

Hot Wire 

LDV 

Cup 

Hot Wire 

LDV 

1 

5.150 

5.154 

5.451 

.760 

.654 

.770 

2 

5.535 

5.543 

5.677 

.736 

.600 

.847 

3 

5.425 

5.479 

5.722 

.940 

.760 

.977 

4 

6.052 

6.092 

6.463 

.813 

.744 

.856 

5 

5.381 

5.406 

5.742 

.714 

.586 

.692 

6 

6.417 

6.426 

6.698 

.822 

.809 

.879 

7 

6.417 

6.457 

6.821 

.702 

.659 

.707 

8 

5.958 

5.996 

6.218 

.745 

.614 

.675 

Averages 

5.792 

5.819 

6.099 

.799 

.678 

.800 


The average wind speed indicated by the LDV measurements is about 
5 percent greater than that indicated by the cup anemometer. This is 
comparably about the same as for Test 50801. The variance for the LDV 
is greater than for the anemometers. Also, the variance for the hot 
wire is less than that for the cup anemometer as was the case also for 
Test 50801. 
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Probabilitii distributions - The distributions of velocities about 
the means for the three instruments are shown in Figure 47. The 
turbulent fluctuations are more normally distributed about the mean 
than was the case for the previous two tests. As before, the probability 
distributions compare favorably one instrument to another. 

Sgeotral densities - A comparison of the spectral density distributions 
with frequency for the three instruments is shown in Figure 48. The 
spectral distribution for the cup anemometer drops off slightly at about 
0.5 Hz, the hot wire spectrum decreases on a constant slope and the LDV 
spectrum tends to level off for higher frequencies. The 2.5 and 5 
hertz spikes were not included in drawing these spectra. The comparisons 
are reasonable to about 1 Hz frequency. 

Frequency, tvaeker - Considerable difficulty was experienced in 
tracking the LDV output with the frequency tracker. The tracker 
required frequent adjustments during the test, and tracking was often 
lost. Consequently the tape recorded output was too intermittent and 
analysis was difficult. 

From observations during the test, it was noted that when tracking 
was achieved, the D.C. output (although slightly nonlinear) corresponded 
with the mean Doppler frequency, hence with the indicated wind speed. 

The A.C. output however did not correspond very well with the turbulent 
fluctuations. For example, in Figure 49, is shown a simultaneous trace 
of the hot wire and the A.C, output from the tracker for Test 101401. 

The hot wire leads the laser focal volume by about 3 meters and the 
average wind speed was about 6 meters per second. The horizontal sweep 
on the oscilloscope was 0.2 sec/cm. 
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Figure 47. 


Distribution of velocities about the mean. 
Test 101401 



Spectral Density (m/sec) /Hz 



0.001 0.01 0.1 I 10 

Frequency in Hz 

Figure 48. Comparison of spectral density distributions. 

Test 101401 
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Figure 49. A.C. Tracker and Hot Wire Traces . 
Test 101401 


The A.C. output (top trace) resembles noise rather than turbulence, 
while the hot wire output is clearly that which traces the turbulence. 
The intermittency of the tracker signal created considerable difficulty 
with digital data analysis. After considerable effort, this part of the 
data analysis was abandoned. The particular frequency tracker used in 
these tests (1971) should be modified to provide long-term uninterrupted 
velocity-time histories. This of course is related to Doppler S/N 
ratio and to the concentration of aerosols which provide the Doppler 
shifted signals. With no Doppler signature (signal drop out) there can 
be no tracking regardless of the quality and design of the frequency 
tracker. 

Measurements of Run 102501 (October 25, 1971) 

Test time was from 2:04 pm to 2:45 pm. The wind was from the 
south-southeast at about 5 m/sec. There were no active weather fronts 
in the vicinity and the sky had been clear for the day. Some pollution 
was evident in the air, but visibility was good. 
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Spectrm analuzer setting s - The settings were as follows: 

Sweep rate: 5 ms/cm 

Sample rate: 16 sweeps/ sec 

Frequency Dispersion: 0.2 MHz/cm 

Filter Bandwidth 10 KHz 

Bandwidth: 30 KHz 

The calibration frequency was 1.678 MHz as shown in Figure 50. 

The noise level from the detector is shown in Figure 51. The vertical 
scale is 200 mv/cm. A sample Doppler trace of one sweep is shown in 
Figure 52. As is observable, the S/N ratio is small which made data 
analysis difficult. 

Velocity time traces - Time traces of velocity from the cup and hot 
wire anemometers and the LDV are shown in Figures 53 and 54. There was 
much more variability of wind speeds during this test than in previous 
tests. The smaller scale turbulence is superimposed on larger scale 
variations. Thus, it should be expected, as will be seen later, that 
the power spectra would indicate greater power at the lower frequencies. 
Some amount of dropout in signals is indicated for the LDV. In general 
comparisons of the time traces appear satisfactory. 

Means and variances - The means and variances for 8 segments of a 
34-minute time period are given in Table 6. 



Figure 50* Calibration frequency 1.678 MHz. 
Test 102501 


Figure 51. Noise Calibration. 

scale is 200 ms/cm. 
Test 102501 



Vertical 


Figure 52. Sample Doppler signal. 
Test 102501 
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Figure 53. Time traces of wind velocity. 

Test 102501, Interval 1 

{For means and variances see Table 6) 
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Figure 54. Time traces of wind velocity. 

Test 102501 , Interval 2 

(For means and variances see Table 6) 
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TABLE 6. MEANS AND VARIANCES FOR TEST 102501 


4.26-Minute 

Intervals 

Mean 

Velocities m/sec 

Variances (m/sec) 

2 

Cup 

Hot Wire 

LDV 

Cup 

Hot Wire 

LDV 

1 

4.397 

4.444 

4.298 

1.077 

.900 

.984 

2 

4.154 

4.169 

3.946 

1.372 

1.273 

1.418 

3 

6.025 

6.010 

5.805 

.762 

.482 

.762 

4 

4.943 

5.000 

4.683 

1.450 

1.162 

1.436 

5 

5.307 

5.315 

4.989 

.992 

.717 

.921 

6 

4.713 

4.748 

4.252 

.873 

.710 

.953 

7 

5.082 

5.102 

4.878 

.933 

.702 

.968 

8 

5.278 

5.284 

5.004 

.628 

.385 

.666 

Averages 

4.987 

5.009 

4.732 

1.011 

.792 

1.014 


The average wind speed indicated by the LDV is within 5 percent of 
the cup and hot wire averages- The comparison is reasonably good. 

Probabili ty distributions - The distributions of velocities about 
the means for the three instruments are shown in Figure 55. Turbulence 
velocities are skewed to the right. The distributions are about the 
same as for the other tests. 

^eotral densities - The spectral distributions of turbulence are 
shown in Figure 56. As was noted earlier the lower frequency variations 
of velocities produced greater power spectral densities at the lower 
frequencies. The cup anemometer response drops off at about 0.5 Hz, 
and the LDV tends to level off for frequencies greater than about 2 Hz. 
The comparison of spectral distributions is reasonably good. 
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Standard Deviation m/sec 

Figure 55. Distribution of velocities about the mean. 
Test 102501 



Spectral Density ( m /sec )/ Hz 



Figure 56. Comparison of spectral density distributions. 
Test 102501 
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OBSERVATIONS AND CONCLUSIONS 

As a consequence of the comparisons presented, the following 
observations can be made regarding the one-dimensional LDV system. 

1. The gross features of atmospheric phenomena in the boundary 
layer are measured by the LDV system. The time traces show 
reproduction of these gross features and comparison with 
other anemometers are favorable. 

2. Mean values determined from the LDV data are in general within 
5% of other anemometer data for long (34-minute) time periods. 
The variations are larger for shorter time periods, chiefly 
because of larger variations in measured velocities. That the 
LDV measures larger velocities is also indicated by the proba- 
bility (percent) distributions of the data and by the spectral 
distributions with frequency. 

3. The confidence of measuring high frequency turbulence (greater 
than 2 Hz in atmosphere) is not yet established. 

4. The technique for data reduction of the LDV data is cumbersome 
in its present form. Immediate improvements can be made by 
including on-line analog to digital equipment including a 
special purpose minicomputer to calculate the velocities from 
the digitized data. Alternatively an analog system to detect 
Doppler signals such as an improved frequency tracker could be 
used. The frequency tracker used in this study required 

very fine tuning, and dependable frequency lock was not achieved. 
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APPENDIX A 


A-l Computer Program for Analysis of Doppler Signals 


A-2 Computer Program for Determination of Velocity Profiles 


A-3 Computer Program for Determination of Temperature and 
Humidity Profiles 



APPENDIX A-l 


Computer Program for Analysis of Doppler Signal 



PROGRAM 


LASOOP 


TRACE 


PAGE 


I 


c 


1 * 


IS 


2 * 


2 S 


3 fl 


3 S 


A* 


AS 


5 " 


5 S 


CQC 6*00 FTn VJ.0-P261 OPT»0 02/10/72 13.0K02* 

PROGRAM LASOOP < INPUT .OUTPUT •TaPE5*1NPUT •TAP£6*0UTPUTiTAPE2* TAPE1 » 

«TAPE3> 

COM MON/ BLOCK 1 /LENARR1 ♦ W [NO I HE I 1001 * NRECOR t tNF ILF 1 . 

# ZEROTH 1 tOIRMlRR ( LOQ ) * VOLT (2.100) .WKlDAfl 
COMMON/flLOCK2/LENARRi>,SYNCUSoO> • YLASER <1500 ) .NREC0H2. NF 1LE2. 

# ZLROTM 2 *WRlOAT 2 «NrAPL 2 
COMMON/BCALIRR/SLOPE 12) •ZERO! aP 1 2) .SLOPEAN. AN INTER tSLOPEHWt 

SLGPEwO. WDINT£ p.SLOP£MD. 0M1NTER 
*COMMON/BT APEC A/NCHANTP • NCAL V AL * V AR 1 TP * ACT VOL t IS) 

COMMON/81 NSTCA/N INSCAL • VAR l lN*f ULScAW .ZEROWD.FULSCAM .2C-R0M0. 

„ FWSTInT 

C0MM0N/8L ASER/NKL Y BAC *NPT 5WR 

COMMON/BLASCAL/NCALREC.NSrfPHEC. I BEGChK , C ALE vEL .CALvELO. VAVLEN* 

# OtVFREQ 
COMMON/BNOrSCA/FLYBACK .NO ISRLC • KNLEVEL < 275 > 
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. TIMA0j3tT IMAUJ4 

COMMON/BSOPT2/ IHEOSK2. IS* IP2 

COMMON/BSrtPEOl/LPACOA) * IOEimT 1 .NFLSKP 1 ♦nRCSkPI 

C0MMQN/BSKP£02/LPAC0a2. NTOTF 12. IDENT2»NFLSKP2.NRCSkP2* 

# MTOTAPEf ILA1 f 2.NT0 Th£C*T IMADJ 
COMMON/HCONSTM/NAYEM | R.DIRECMb .CHG rt iR. t IMEMIR 

COMMON/BVOLT AO/ [SCALE 

COMMON/BSPFEC>/SU M VELO ♦ I SAMPLE. 1DATAHW. SUM VOL T *T IMRAT 1 * CM ANNE 1 • 

„ D I GHA f 1 . T I M ECHO • VOLT CHG.MULT IME • F I MEHW.OCSUPRE 

„ ♦F«STSPOf*PlTAPe,PClNTOK 

COMMON/8 AVF W IN/ JSA^PLE.SUM* I NO* MULT IM1 » TIM AV*0 • AVEWD 1 3000 > 

. tLASTIME 

COMMON/BBUMPUP/T IMEI (702) . VELnC2t?02> 

COMMON TIME2I200) . VtLOL AS 1 200 > *IP0)NT 
COMMON/UNPK 1 / I TIME 1 . I COMWRU I 2ftl ) . I JATWPD 1 1 000 1 
COMMON/UNPk2/ITIME2*LCOm#«D< 601) .LuAl WRD I 3000 > 

RE AD (5 . 1 1 VQLTChG*WRIDAT2*CALTAPE.LALIN$T » CALL AS • 

„ CALNOI S. WR IDAT 1 t IUE lT 1 • luENT2.LENAHR 1 ♦LENARR2* 1BEG5K1 . , 

, ISK IPI 1 1 HEGSK2. IS* Io2«NAvEMlR*NFLYdAC* 1SCALE .NTOTF 12 • 

^ NTOTAP£.LPACOa2*NCHaNTR.nCALVAL.LPACOA1 »N FOTF II tNCALREC* 

# NSWPREC. 1 HEGCHK *NO I SHEC. NlNSCAL *NL A5REC.0 IGRAT 1 * 

T iMECnGt TIMEHW* SLOP^M^ ,T I MR AT l ♦ CM ANNE 1 1 VAR ] TP • I ACT VOLT 
I III *1=1*5) ,T [M A V2U*C ALE VtL* FLYBACK. Of LYBACf VAR IlN* 

# FULSCAw.ZERGwD.FuLSCAM tZcPOMO. DCSUPRt *DE Vf WEQ * W AVLtN. 

. CHANNE2.DIGRAT2.TIMPAT2 

1 FORMAK 9A3.U l4t/t 11I**.5F6.O«/.2F6.G.SIFS.0I . 7F 6.0 * / . 3F6.0 . 2E9.3 p 

# 3F6.0) 

READ (5*5) TIMADJl. TlMAQ J2.T IM Ai)j3. T IMAOJA.T IMAUA.NREC2.NREC3* 

# NREC^.NCAlFIl.NCALT AP.NKEC 

5 FORMAT l5F10.it 315. 212. 15) „ 

REAO(5*B> FHST APE. FRST INT. F PSTL AS. FRSTNOS, NT APFI1.NL ASF I2.NINSFI1. 

„ NOISFI2.wRITAPEt iRUNNO.MRCONST 

8 FORMAT I4A3.^!3*A3* IG. A3) 

WRITE (6*1 0) IWUNNO 

10 FORMAT (1H1 *5X*|NPUT DATA FUK PUN NUMBER*! 7) 

«R1TE(6.2) VOLT CriG.WRI DA |2 tCALT APE. CAL INST. CALL AS» 

„ CALNO IS*WR iDATi . IOE nT 1 » IuENT?.LENARR l « LENARR2* I BEGS* 1 * 



PROGRAM 
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6^ 


7n 
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9S 
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ir 


UASDOP 


TRACE 


COC 6400 FT* V3.0-PZ6I 0PT»0 02/10/72 13.01.02. 


ISK 1P1 t 1BEGSK2# ISK|P 2 *NAvEHlR«NFLY 8 ACt 1 SCALE *NTOTF 12 « 

NTOTAPEtLPACOA2*NCH4NTPiNCALVAL*LPACOAliNrOTFllfNCAL«ECt 

NSwPRECt ! 0 EGCHKtNOXSREC*NlNSCAL*NLASr?ECfOlGRAn* 
TlHECHG.TlHEHW,SLOPFH^,TlMRATUCHANNElfVARITP. UCTVOLT 
(I) *1 = 1 *5 1 « T lMAVwD« CALL VEL * FLYBACK » DFLY 8 AC * VAR I IN* 
f ULSC AW t ZEkOWD* FULSCAHi Z tROHD.DCSUPRE • DEVFREQ. WAVLfcNt 
CHANNE2.0IGRAT2.TIMP4T2 

FORMAT < 1 HO •* VOLtCHG = *A4* WR TO AT2 = *A4» CALTAPE -*A<*« CALINST * # A 
<♦* CALLAS = »A4« CALNOlS = *«<♦• WR 1 DftT 1 = »A4* 10E.NT1 =*A4» 1 Cj£NT2 «* 
A 4 /* LENARR1 = *| 4 * LENAPR2 =*f4* IbEGSKl = *J3* ISK1P1 =*I3* 10EGSK 
2 =*I 3 <> ISK1P2 = • 13* NAVtMIR = *I3* NlLYdAC = *1<** ISCAlE =*I2/* NTO 
TFl? =*I2* NTOTaPl = *12* LP^CrAZ = *Ih* NCmANTP = M3* NCALVAL »*1Z 

* LPACDA1 = <M4* NfOTFIl =*12* NCALrFC =*13* NSwPHEC =*I ?/• IfclEGCHK 
=•13* NOISWtC =*13° NlNSCAL =*I2* NLASHEC =*IJZ* DiGwAU -*F7.l 

* TfMECnG =*F3.1« TlMtHw = **F9.4* SlOPLHw =*F9.^* T I HR AT 1 =°FA.l/ 

* CmANNEI =*FS.l* VARITP =*Ki/j/ # aCTVOLTO THRU 5) =*SFS.l/* TIMA 
VWO =*FS.2* CALt'VtL =°F 9*2* Ft. Yfi ACk =*fb.l* UFLYGAC =*15. I* VAHIIN 

= »F 6 « 3/* FUl SCAv* =*F9.3 * ZEPOWU =*F^.3* FULSCAM =*F9.3* ZEROMO * 
•F9.3* OCSUPRE = *F 7 • 3 * DEVFhEn =*E10.3/* WAVLEN =*E13.6* CHANNE2 - 
•F5.1* OIGkATZ = *F9.1* TIMRAT2 =*FG*i» 

WP I TE <6*6 > FRST APE ^FPSt IMT t FRStL AS«F MSTNOS tNTAPF 1 1 tNLAST I2*NlN$F Ilf 
NOISF I2t WRI T APEtHRCnNST 

FORMAT 41 H ,* FRSTAPE = *A 4 * FksTINT = *Aa* FRSTLAS = *AA* FRSTNOS =*A 
4 * NT A^F 1 1 =*I3 G NL ASF 1 2 **U° NlNSFIl =*IJ« NOISFI2 = *I3* WRITAPE 
= *AJ» WRCONST =• A3 ) 

i»RITE<6*9)TIMADjl tT JMADJ2 * T IMfiD J3* T IHAOJ4* T IMA0A*NREC2«N«EC3*NREC4 
»rjCAl_F lL ♦NCAlT AP * NR fC 

FORMATHH .* T1MAGJ1 =»FA.i« TIHADjZ =*Fi0.1* T IH ADJ3 =»F10.1* TIM 
A0J4 =*F 10,1* T I M ADA s«F10.1/» NkEC2 =*15* NR6C3 =*I6* NPEC4 »*I6 

* NCALFIL =• 1 2* NCALTAP = *I2* NREC = *IA) 

REWIND 1 


JHYES ) REWIND 3 


REWIND Z 
IF (WRITAPE *EQ « 

FRSTPT = 3HYES 
IEKTIME = 0 
MUlTlMt = 1 
SUMVELO = 0.0 
ISAMPLE - 0 
SUMVOLT = 0.0 
JSAPPCE = 0 
SUM WIND =0.0 
CAIVELO = 0.0 
NSWPS - 0 
NWECOR|= 0 
N»ECOR2=0 
NTRIG = 1 
NFIlEl = l 
NFILE2 = 1 
2EPOTM1 = 0.0 
ZEROTM?=0.0 
IE A I T = 3H NO 
IEXIT2 = 3 m NO 
SLASTPT = 10.0 

I Change = o 

FACTORl = SORT 12 # )/(2.**9- 


1 . 0 ) 
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PROGRAM 


UASOOP 


THAC(E 


COC 6400 FTn V3.0-P261 OPT*0 02/10/72 13.01*02 



\\* 


120 


12 * 


13 * 


13*5 


140 


1 <*S 


IS* 


ISS 


160 


16 * 


FRSTSPO = 3HVES 
NTAPE2 * i 
NEXTPTS - 0 
NAVEWD = 0 
NOATAPT = 1 
JCL0CK1 = 0 
JCL0CK2 = 0 
. IWlND = 1 
IPO! NT = 1 
JCOUNT = 0 

IF M0ENT1 .FO. 3MYES) CALL HEADER 1 
IF (IDLNT2 .EQ. 3NYES) CALL HFA0ER2 
IF IFRSTAPE .EQ « 3h OK) GO TO 20 
NFLSKP1 - 0 
NPCSKP1 - 1 
CALL SKPE0F1 

20 IF ICALTAPP .EQ. 3HYES) CALL TAPECAL 

IF 1CALTAPE .EQ. 3HNE0) REA0<5*11) ( SLOPE I I) *2ER0TAP ( 1 ) t I»1 *2) 

11 FORMAT <AF10. 3) 

IFINFILE1 .GT ♦ NTAPFll) GO TU 21 
NFLSKP 1 * i 
NPCSKP) = 0 
CALL SKPE0F1 

21 IF K Al INST .EQ, 3HY£S I CALL INSTCAL 

IF (CALINST .EQ. 3HNE0I READ IS * 11 ) SLOPEWOfSLQPEMD* POINTER • OH INTER 

IF (NFILEI *6T ♦ N1NSF11) GO TO 23 

NFLSKPI = 1 

NRCSKP) x 0 

CALL SKPE0F1 

23 IF IFRSTLAS .EQ. 3H 0K> GO TO ?4 
NFLSKP? a 0 
NRCSKP? = 1 
CALL SKPE0F2 

24 IF KALLAS .EO. 3HYES) CALL LASCAL 

IF KALLAS .EQ. 3MNE0) READIS.12) CAL VELO.NPTSPP 

12 FORMAT (F10. 3. Ul 

FLTPACK = FlT d ACK - DFLYBAC 

IF (NFILE2 »GT , NL ASF 1 2 ) GO TO 25 

NFLSKP? x i 

NRCSKP? = 0 

CALL SKPEOF? 

2S IF (FRSTNGS .EQ. 3H OKI GO TO 26 
NFLSKP? x 0 
NRCSKP2 - 1 
CALL SkPEOF 2 

26 IF KALNOIS .EQ. 3MYES) CALL nOISCAL 

IF KALNOIS *E0. 3HNE0) REALMS, 41 < XNLEVEL II I ♦ 1 *1 ♦ 202) 

4 FORMAT < ( 13F6.0) I 

IF (NF1LE2 .6T. N0ISFI2I GO TO 27 
NFLSKP2 = I 
NRCSKP? - 0 
CALL SKPEOF? 

27 CALL VOLT AO J 

IF < MRCONST .EQ. 3HYtS> CALL COnSTHR 
IF t MRCONST ,NE. 3HNE0) GO TO 75 



PROGRAM 


17 ft 


175 


180 


1*5 


19 ft 

'•<] ,« 


200 


205 


21 ft 


215 


22 ft 


lasoop trace 


COC 6400 FTN V3.0-P261 OPT»0 02/10/72 13*01.02* 


READ 45*14) CHGM1«*T1MEMIR 
READ (5* 1 3 ) OIRtCMR 

13 FORMAT (F10. 3) 

14 F0RMAT(A3»F6.2J 

WRITE I6» IS) DIRECMR»CMGM IRt TI mEH IR 
15 FORMAT UH0*5X*DIRECmR =*F10.1* CHGMft -*A3* TlMCMlR **f 10.31 
75 JCLOCK ! * 0 
ZER0TM1 = 0.0 
JCL0CK2= 0 


16 

100 


7 


ZEP0TM2 - 0.0 
IDATAhW - i 
MULTJH1 = 1 
IEXT1ME = 0 
LASTIME * 1 
MULTIME = 1 
PRINT 16 
FORMAT ( 1H1 > 

CALL SPEED 
CALL AvEWlNO 

JFIPRIimTOK *EQ. 3HYES1 WR ITE < 4 * 7 ) NF 
FORWAT ( 1 HQ *5X*LASER VELOCIT I tS*5X*F 


ILE2*NT APE2 

lLfc*12«SX*T APEM 2/ 10X*T IMEtSEC* 


,10X*VEL0CITY*M/SEC* 10X*RECQftD* 1 
150 CALL BUFLAS2 
ISEGIN = I 

175 00 ?QG M=IBEGIN*LENARR2 

IF <YIASER<M> *GE. FLYBACK) Go TO 300 
200 CONTINUE 
GO TO 150 

300 J COUNT = JCOUNT ♦ 1 
M a M ♦ 1 

IF <M .LE* LENARR2) GO TO 500 
400 CALL BUFLAS2 
M= I 


500 IF < YLASER (M) .GE. FlTBACKI GO TO 300 
IF fJCOUNl «6T • 15) GO TO 600 
JCOUNT = 0 
1 BEGIN = M 

IF (M .LE* LEN ARR2 ) GO TO 175 
GO TO ISO 

600 M = M * IBEGCMK - 1 

IF (M *Lt* LENARR2IGO TO 650 
CALL BUFLAS2 
M = M - LENARR2 
650 XOUNT =• 0 

LAST = NPTSWP W NFL Y8AC 
DO 800 1= IrtEGCHKt LAST 

IF tYLASEH(M) *GE. XNLEVtL < 1 ) # 2.)G0 TO 900 
M= M ♦ l 

IF tM *LE. LENARR2) GO TO 800 
700 CALL BUFLAS2 
M s 1 

800 CCNTINUE 
I8EGIN s M 

IF <FRSTPT #EQ ■ 3HYE5) GO TO 1 75 
VELOLASUHOlNTI s yelolasupotnt - H 
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PROGRAM 


22 ^ 


230 


235 


240 


24=; 


?S0 


2SS 


260 


255 


LAS0OP 


TRACE 


COC 6400 FTn V3.0-P261 OPT*0 02/10/72 13*01.02. 


910 

920 

825 

675 

900 

925 


TIHE2 < IPOlNT ) s TINE2 < IPOlNT-1 > ♦ |NPTS*P*CHANNE2*T IMRAT2I /0JGRAT2 

IF (T |ME2< IPOlNT) .GT • IWIN0*TIMAVWD) 1WIN0 *1WIND ♦ l 

IF <AVE*D<IWINC>) .NE .0 .0 1 60 TO 62b 

IF < 1 EX I T .EG. 3HYES) 60 TO tt?0 

IF (WRIT APE .EQ* JHYES) CALL I.ASWRIT 

IPOlNT =1 

CALL SPEED 

CALL AVCWlND 

IF (PRINTOK .EG. 3HYES) WRI TE <6. 7MFILE2.NTAPE2 

60 TO 825 

IWIND - 1 W I nD - 1 

IF IAVEWDUWINO) .EQ. 0 . 0 1 GO TO B20 

IF (PRINTOK .NE. 3MYES) 60 TO 975 

WPITE (6*3) T IME2 ( IPOlNT > , VfcLOnS UPUlNT 1 tNR£C0R2 


IPOlNT = I PO I NT 
60 TO 175 
IF (YLASER(M*1) 
I = I ♦ 1 
H * M * 1 
60 TO 900 
TIME2UP0lNT> - 


l 


• LT. YLASER {(A) 1 GO TO 925 


TIHRAT2*( ( I T ImE2~Z£RGTH2> /I 0000* < CM 
* DIGRAT*) 

FRSTPT = 3H NO 

IF (T 1ME2 I IPOlNT ) .GE. TINEHIP .A. CHGMfR .EQ 
IF (T JM£2 ( IPO INT) .61. I WIND * TIMAVWD) IWlNO 
IF (AVE40UWIND1 .NE. 0.0160 TO 93b 
IF C I EX I T .EQ. 3HYES* 60 TO 930 
IF (WRIT APE *E0. 3HYES) CALL [ ASWRIT 
IPOlNT =1 
926 CALL SPEED 
CALL AVEWINU 

IF (PRINTOK .EQ. 3HYESI WRITE <6* 7 )nFILE2>NTAP£2 
GO TO 935 

930 I WIND = I WIND - I 

IF (AvtWDUWIND) .EQ. 0.01 GO TO 930 
935 WDIREC = AvEwOUWlNO) - DlMEC^R 

WOIPEC = (wQIREC * 2. * 3.141/ 3b0. 

VELOLAS ( IPOlNT) = HI ♦ 41 *CAl YEL0)/C0S(WOIREC» 

IF (PRINTOK .NE. 3HYESI GO TO 1000 
WRITE (6*3) TIME2 (IPOlNT) .VELO laS* IP oINT) .NREC0R2 
3 FORMATdH , 8XF6.3»1SXF6.3»14<I4J 
1000 IPOlNT = IPOlNT ♦ 1 
I BEGIN = M 

IF (M .LE. LENARR21GO TO 175 
GO TO 150 
END 


1 *CHANNE21 / 


3HYESI CALL CONSTHR 
1WIND * 1 
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SUBROUTINE TAPfCAL TRACE 


CDC 6400 FTN V3. 0-P261 0PT»0 02/10/72 13.01.02 


PAGE 


1 




I* 


1* 






2* 


3n 


3S 


4" 


4«; 


so 


ss 


SUBROUTINE T APE CAL 

COHMON/BT APECA/NCHANTP»NCALVAL* VAR ITP* ACT VOLT (S) 
COMMON/BLOCK1/LENARR1.*IND1RL( 100) »NREC0R1 tNF ILE 1 t 
• ZEROTMltDrKMlRR(100)«VOLT<2»100> tWRIOATl 

COHHON/BCAL IBP/ SLOPE (2) *ZER0TftP(2) .SLOPE AN* AN INTER • SLOPE HV. 

. SLOPEWD. VOInTEo.SlOpEMD.DM INTER 

DIMENSION SUM CAL ( 2)‘.SUMT AP (2 ) «SQ VAlUE 421 .SUM ACT (2) .ACT X TAP(2>. 

. SUM5Q (2*5) .WECMt AN (/) . TOT ME AN < 2 *5) • TEPMEAN (2) .SUMEAN(2) ■ 

TEMPSUM ( 2 ) .ST ANDEV ( ? . 5) 

1CMECK = 0 
NSAMPLE = 0 

lastcal = 0 
icalval = i 

00 100 la 1 . NCHANTP 
SUMEANU) = 0.0 
TEMPSUMU) S 0*0 
SUMCALU) = 0.0 
SUMACTU) = 0.0 
SUMTAPm = 0.0 
SOVALUEm - 0.0 
AC I X T AP ( I > a 0.0 
TEPMEANII) = 0.0 
RECMEANUI = 0.0 
DO 100 j=i,ncalval 
TOTMEAnU»JI - 0.0 

100 SUMSQCI.J* = 0.0 

110 call BUFLAS1 

GOTO0UF = 3H NO 

IF (ICALVAL .tQ. NCALVALI LASTCAL = LASTCAL ♦ 1 
DO 120 1=1. NCHANTP 
00 120 K=l,LENARRl 

!20 SUMCALU) = SUMCALU) ♦ VOLTU*K> 

IF (ICHECK .CT . 01 GO TO 1S1 
NSAMPLE a NSAMPLE ♦ 1 
00 140 1=1. NCHANTP 
RECMEANU) = SUMCAL(I) /LENARpI 
IF (NREC0N1 .EG. II GO TO 125 

IF (RECMEAn(I) .GT. TOTME AN U . ICALVAL ) ♦ VAR1TP .0. RECMEANU) 

■ .LT. TOT ME AN (I ♦ ICALVAL) - VARITP) GO TO 150 
12S DO 130 K=1.LENARR1 

130 SUMSOUUCALVAU = SUMSOlIf ICALVAL) ♦ VOLTUtK)*** 

SUMEAN(I) a SUMEANU) * RECMEANU) 

SUMCALU) = 0.0 

TOTMEANtlt ICALVAL) * SUMEAN (I ) /NSAMPLE 
IF U *EO • 1) 

. WRI TEI6tl)N«EC0ft I * ICAL V AL* ACT VOLT ( ICALVAL) 

1 FORMAT UHOtSX«RECORL) ME ANS*4x*R£C0R0 NUMBER*I4* 7X*CALIBRAT I0N*I2* 
.4XMNPUT VALUE *F5*l/ll X*CMA nnFL*1 0a*HEAN*1 3X*CUMJL AT I V£ MEAN* 6X 
.•NUMBER RECORDS FOR CUMULATIVE MEAN*) 

140 WRITE lo»2) 1 .RECMEANU) .rOTMEANUUCALVALi *NS AMPLE 
2 FORMAT (1H tl2A*J2'10X*F8.4.l«4X.F6.4*25X.I3) 

GO TO no 

J50 NSAMPLE = NSAMPLE - 1 
ICALVAL - ICALVAL ♦ l 

1S1 IF (ICALVAL .GT. NCALVAl .A. i.ASTCAL .GT. 3) GO TO 180 


l 





SUBROUTINE tapecal TRACE 


COC 6400 FTN V3#0~P26I OPT*0 02/10/72 13.01.02* 


60 


6S 


70 


75 


HO 


85 


9n 


95 


100 


105 


no 


ICHECK * ICHECK ♦ 1 

*RITE<6*5>NR£C0R1»!CALVAL»ACTv 0LTUCALVALI 
5 FORMAT <1 HO *5X* TEMPORARY MEANS*dX*RECORD NUMBER* 14*1 0X*CALI8RAT ION 

« * 1 2 * 4x* I nput value •F5.i/nx*cMANNtL*iox*MEAN«i 
DO 170 1 = 1 #NCrtANTP 
RECMEAN<!> - SUMCALIU/LENARRI 
MR I TE ( 6* 6 ) I * RECmEAn i I > 

6 FORMATUH , 12* I 2 • 1 0XF8.41 
SUMCALlU =0*0 
IF l I CHECK .€0. II GO TO 160 
1)0 155 K=1.LENARR1 

155 SUMSQ ( I f ICAL VAL I - SUMSQ ( I ♦ ICaLVAL ) ♦ V0LT<IiK)**2 

TEMPSUM ( 1 1 = T EMPSUM ( 1 1 ♦ RECmEANUI , . k 

160 IF (RECMEAN til *GT . TOTMEAN U . 1CALV AL- 1 » ♦ VARITP .0. RECMEANlII 
„ ,LT. T0TMEAH(I*ICALVAL-1>-VAR1TP1 GOTO0UF = 3MYES 

170 CONTINUE 

JFUCHECK *GT . 3) GO TO 180 
IF (GOTOBUF .tO. JrirESI GO TO 110 
00 ] 75 1=1 *NCH ANTH 
TEMOSUMU) = 0.0 


175 SUMSOUtICALVAL) = 0.0 
ICHECK = 0 

ICALVAL = ICALVAL - l 
GO TO HO 

]60 lEND = ICALVAL - 1 

*RlTEt6*&MEND.ACTV0LTUENU> ua 

8 FORMAT ( 1H0. /»5x*5TAN0AR0 DEV I *T ION5 # 10X*CALtBRAT 10N*I2t5X*INPUT VA 
# LUE*F5.1/1 1X*CHANNEL*10X*KM$*> 

00 190 1=1 *NCHANTP 

STAKOEV (I t lENOi * SORT < 5UMSO I T • l£N0) / INS AMPLE*LENARR1 1 - 
, TOTmEAN (I * IEnDI **2 1 

190 IfRlTE 16*91 I *STANOEV(It I END I 

9 FORM AT < 1 H . 12X12WXF9.3I 
NS AMPLE = ICHECK “ I 

DO 195 lel.NCHANTP 
SUMEAN(l) = TtMPSUMU) 

TOTMEAN I I 4 ICALVAL) = TEMPSUM I f I /NSAMPLE 
195 TEMPSUM U> = 0.0 
ICHECK = 0 

IF UCALVAl .LE. NCALVAL) GO TO 110 

10 FORMAT n HO !«!"* ACTUAL VS TAPE vOL T AtiE* 10X*LE AST SOUAftE H£TH00*5K*NU 
*M8ER RECORDS USED FOR CALCUL AT IONS* 131 

OO 210 I=1*NCHANTP 
DO 200 J=1«NCALVAL 

SUMTAP(I) = SUMTAPIl) ♦ TOTMEAN ( I f J> 

SOVALUEU) = SQvALUE(I) ♦ TOTmEAN l I ♦ J) **2 
ACT X TAPU) = ACT X TAP(I> ♦ TOTHE AN < 1 1 J) • ACTVOLT ( J1 
200 SUMACT(I) = SUM ACT ( 1 1 ♦ ACTVOlT(J) 

SLOPE < I 1 = (SUM ACT (II *SUMTAP(U -NCALVAL* ACT X TAPUM/ 

. (SUMACT (11**2- NCAl VAL*SOVALUE (1)1 

ZEROTAPin = I SUMACT < I I *ACT X TAP(]I - SUMTAP t l > *SQV ALUE < I II / 

( SUMACT (11**2- NCAl VAL # SUVALUE( I) 1 
’•RITE ( 64 ID I* (ACTVOLT i J) 4 TO IMF An ( I . J > .J=l .NCALVAL) 

11 FORMAT (lHO . 10X*C H A N N E L * l 3/ 1 5x *V ALOES USED FOR LEAST SQUARE C 
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SUBROUTINE TAPECAt TRACE 


COC 6400 FTN V3.0-P261 0PT*0 02/10/72 13.01.02 


.ALCULATION$*lOX*INPUT VALUE *Sx* TAPE VALUE*/ t 69 AF 4.1*11 XF 6.3)1 
210 WRITE (6 ♦ 12) SLOPE 1 1 > *Z£ROT AP < 1 1 

12 FORMAT (1H0.15X*VALUES OBTAINED FROM LEAST SQUARE CALCULATIONS* 7X* 
*slope*2x* intercept */68Xfs*3. i i xfs. ji 
print 13 

13 F0RMATUH1) 

RETURN 

End 



SUBROUTINE 


10 


1* 


2f> 


25 


3n 


35 


4* 


45 


5* 


55 


INSTCAL TRACE COC 6400 FTN V3.0-P261 OPT»0 02/10/72 13*01 *02* 

SUBROUTINE INSTCAL 

COMMON/emSTCA/NINSCALtVARlIN,FULSCA^f2EROrfD#FULSCAM*2ERONO* 

FRST I NT 

*COMMON/8lOCk 1/LENAHR1*WINOIRE| I00> .NRECORl *NF ILU • 

. ZEROTMl'DIRMlRRf 100) * VOLT 1 2* 100 1 1 WR IOAT 1 

COMMON/BCAL IBR/SLOPE (2) *ZEROTaP( 2I » SLOPE AN* ANlNTER*SLOPEH¥« 

# SLOPE^O« WD INTER tSLOPENO* DM INTER 

DIMENSION SUMSQWOO) .SUMSQMO 1 3) *TH£AN*DI3> t THE ANNO C3l •$TDEVWD(3J * 

STDEVMDC3) 

LEVEL = 10HZERO INPUT 
NRECORl = 0 

IF (FRSTINT .EQ. 3HTES) NRECORl = l 

SUMWO = 0.0 

SUMMD = 0.0 

SUMAVE* = 0.0 

SUMAVEM - 0.0 

N5AMPLE = 0 

1NSCAL = 1 

TEMPWD = 0.0 

TEMPMO * 0.0 

ICHfCK = 0 

DO 100 I * I.NINSCAL 
SUMSOVD < I > = 0.0 

ioo suMsanom = o.o 

IF (NRECORl . EO. n GO TO 175 
150 CALL BUFLAS1 

GOTOHUF = 3H NO 
IF CNF 1 LEI .GT. 1) GO TO 850 
175 DO 200 K=1.LENARR1 

SUMWO = SUM WO * WINDIRE(K) 

200 SUMMD = SUMMU ♦ DIRMIRR(K) 

AVEWO * SUMWU/LENARR 1 
AVEMD = SUMMD/LENARR1 
SUMWD = 0.0 
SUMHD = 0.0 

IF (NRECORl .EO. II GO TO 300 
IF ( I CHECK .GT. 0) GO TO 600 

IF (AVE*0 .GT* TMEANWDUNSCAL) ♦ VARIINI 00 TO 500 
300 SUMAtfEw - SUPAvEW ♦ AVEWU 
SUMAVEM = SUMAi/EM ♦ AVEmO 
NS AMPLE = NSAMPLE ♦ 1 
TmFANWDUNSCALI = SUMAVFW/NSAMPLE 
TMEANMDUNSCALI = SUMAVEM/NSAmPLE 

*RI TE <6* 1 INRECORI » LEVEL * AVEMO*TMEAnMD< INSCAL > •NSAHPLE tAVEWO* 

. tme anwu ( inscal i ♦ ns ample 

1 FORMAT < 1 HO. 5K* INSTRUMENT CAL IP«AT 10N*5X*RECORO M£ANS*5A*RECORD*l3f 

*5X*INPUT *A i 0/1 QX* INSTRUMENT*! 0X*RECOR0 ME AN* 1 0X*CUMUL AT IVE MEAN* 

.^♦NUMBER OF RECORDS IN CUMULATIVE MEAN*/7X*MIRR0R 0!RECT10N*9XF6 
.,3*1 9XF 6.3*30X12/ 8X*Wl NU D 1 RfCT ION* 1 0XF6. 3* 19XFO. 3 . 30X 12 > 

00 400 K=1*LENARRI 

SUMS0W04 INSCAL) = SUMSQWD ( INSCAL ) ♦ tf lND IRE IK ) **2 
400 SUMSUMOC INSCAL) = SUMSQMD (INSCAL ) ♦ DIRHIPRIK) **2 
GO TO 150 

500 INSCAL = INSCAL * I 
600 ICHECK * IOECK ♦ l 
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SUBROUTINE INSTCAL TRACE 


COC 6400 FTN V3.0-P26X 0PT*0 02/10/72 13. 01,02 


60 


65 


TO 


75 


60 


65 



4r> 


95 


100 


10s 


no 


IF (ICHECK .LT. 21 00 TO 150 
TEMPWD = TEMPWD ♦ AVEWO 
TEMPMO = TEMPMD ♦ AVEMD 
IEND = ICHECK - 1 
LEVEL = 10HFULL SCALE 

WRI TE <6i 2)NREC0Rl »LEyEL •AVEMO*T€MP mO* IEND« A vtUOf TEMPwO* IEND 
2 FORMAT ( 1 HO # 5X*TEMP0 RaRY SUMS*1 OX*RECORD* 13* 1 OX*lNPuT *A10/ )0X*1NS 

,tpument*iox*record MEAN* lQX*5UM of RECORD means*! ox*numher of reco 
. ROS IN SUM*/ 7x*M 1RR0R DIRECTlON*9xF6.3,20XF6*3.25Xl2/6X *WINO OIR 
.ECT ION* 10XF6,3*20XF6*3f 25X12) 

00 700 K=l,tENARRl 

5UHSQW0UNSCAU = SUMSOwD < INSCAL ) ♦ V IN01RE IK > «*2 
700 SUMSQMDUNSCAU = SUMSQMU( INSCAL) ♦ D IRM IRR <K > **2 
800 IF <AVE*0 • GT • TME ANrfD < INSCAL- 1 > ♦ VARIINI GOTOBUF a 3HVES 
IF ( 1 CHECK .(ST. 31 GO TO 900 
IF ( GOT OBUf .EQ. 3 hYES) GO TO ISO 
SUMSQWO l INSCAL) - 0.0 
SUMS(JMOUNSCAL) = 0.0 
ICHECK = 0 
TEMPMO = 0.0 
TEMPMO = 0.0 
INSCAL a INSCAL - 1 
GO TO 150 

650 INSCAL = INSCAL * 1 
900 IEND = INSCAL - t 

LEVEL = 10HZER0 INPUT 
IF I IEND .EQ. 2) LEVEL * IOMFuLL SCALE 
STDEWDUEnO) = SQRT(SUMSOtfD(jEND)/<NSAMPLE*LENARRII - 
. TMEANWDUENO}**2> 

STOE VMO ( 1 EnD) * SORT ( SUmSQMO ( t END) / (NSAMPlE*LEnARR1 1 - 
. THE AN«D ( IENO) **21 

WRITE (6*3) I END, LEVEL, STOtVMU UENO> , STOE VWD(IENQ) 

3 FORMAT ( 1 HO • 5x *ST AND ARO DEVI AT !0N5* 1 OX*CAL I BR AT ION* I2»10X* INPUT * 
.A10/10x*lNSTRUMENT*10X*RMS*/7>r*MIRROR DIRECTION* 6XF5.3/ 8X*W1ND D 
. IRECTI0N*7xFe>.3) 

LEVEL = 10HFULL SCALE 
NSAMPLE = ICHECK - 1 
SUMAVEw=TEMP*D 
SUMAVEM = TEMPMO 

TMEANWDUNSCAL) = SUMAVEW/NS AmPLE 
TME ANHO( INSCAL! = SUMAVEM/NS AmPLE 
TEMPWD = 0.0 
TEMPMD = 0.0 
ICHECK = 0 

IF < INSCAL .LE. NIN5CAL} GO TO 150 

SLOPEWD = 1 FULSCAw-ZERO/OI / i TmEANMO (2) -TMEAnwD ( 111 

SLOPEMO * <FuLSCAM-2ER0mD1 / ITmEANMD (2)»TMEANMD (1 ) 1 

dm inter = fulSCam - slopemd*tmeanhd<2> 

WDlNTER = FULSCAV - SLOPE WD* TmEANWO «2> 

WRITE (6»4)NREC0R1*TMEANMDU> t7ER0M0*TMEANMD(2) •FULSCAMiSLOPEMO* 
.0M1NTER*TMEANWD< n fZEROWDiTMEANWD< 2 ) vFULSCAW* 5LOPEWD* POINTER 

4 FORMAT <1 HO tSX* I nstrument CAlIpRATION*SX*NUMBER of rtECOROS USED* 13/ 
. 1 0 A*M I RROR UIRECTI0N*/15A*VA| UES USEO FOR CAL JBRAI ION* 1 OX* 1NPUT*S 
.X *T APE VALUL*5X*ACTUAL VALUE*/ S3x*2ERO*6XF6 . 3t 10 XF 7*3/ A9X*FuLL 
.SCALE*4AF6. J«10XF7*J/ 15X*VALuES OBT A INE0*22X*SL0PE*SX* INTERCEPT*/ 
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«SIXF7.3»6XF5.3//10X*WlND DIRECTION*/ 15X*VALUES USED FOR CAL I BRAT 1 
»0N*10X*1NPUT*5X *TAPE VALUE*5 x*ACTuAL VALUE*/53X*2ER0*bXF6#3*10XF7 
*.3/ 49X*FULL SCALE*4XF6« 3t 1 0XF7* 3/ 15X* VALUES CUT AlNEO*22X*SLOPE*5X 
# *INTERCEPT*/SlXF7»3t6XF7*3l 
RETURN 
END 
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SUBROUTINE LASCAL 

common/8lascal/ncalrec*ns*prec* iheochk •calevel *CALV£LO»WAVLEN» 
. devfreq 

COMMON/0LOCK2/LENARR2*SYNC( 1500) » YLASER 1 15001 *NREC0R2*NFILE2t 
# ZER0TM2**RIDAT2,NTAP£2 

C0MH0N/9LASER/NFLTBAC»NPTS^P 
DIMENSION NPTR1GI20) 

next so 
NPTS*P * o 
SLASTPT = 10. 0 
CALVELO =0.0 
NSAMPLE = 0 
NSHPS = 0 

lastpt = 0 
GO = 3H NO 
leftovr - 0 
XOO CALL 6UFLAS2 

IF (NREC0R2 .GT. NCAlREO GO TO 900 
150 NTRIG = J 

00 50 1=1*20 
50 NPTRIG(I) * 0 

DO 500 K=1*LENARR2 
IF IK .GT. 1) GO TO 200 

IF I SYNC IK ) .GT . 0.0 .A. SLASTPT .LT. 0.0) GO TO 300 

GO TO 500 

200 IF ISYNCIK) .GT. 0.0 .A. SYNC(K-l) .LT. 0.0> GO TO 300 

GO TO 500 

300 IF I GO .EG. 3H NO) GO TO 400 

NPTS*P = NPTStfP ♦ K » LASTPT ♦ LEFTOVR - I 
NS*PS * NSWPS ♦ 1 

leftovr = 0 

400 NPTRIG (NTRfG) = K - 1 
NTWIG = NTRIG ♦ 1 
LASTPT = K - 1 
GO = 3HYES 
500 CONTINUE 

SLASTPT = SYNCILENARR2) 

LEFTOVR = LENARR2 - NPTR I G (NTRIG- 1 ) 

LASTPT = 0 
NEXT = 0 
NTRIG = I 

NAVEPTS = NPTSWP / NSWPS 
600 ISTART = NPTR IG INTRIG) ♦ IBEGCHK 

LAST = NPTRIG(NTRIG) ♦ NAVtPTS - NFLYBAC 
IF (LAST .GT. LENARR2) NEXT = LAST - LENARfi? 

IF (LAST .GT. LEnARR2)LAST«LEnARR 2 
DO 700 I = ISTART*LAST 
IF (YLASER(I) .GT. CALEVfcL) NFAT * 0 
If (TLASER(I) .GT. CALEVEL) Go TO 000 
700 CONTINUE 

IF (NEXT .EU. 0) GO TO 750 
CALL BUFLAS2 
DO 72S J= 1 .NEXT 

IF ( YLASER I J) .LE. CALEVEL) Gn TO 725 
I=J*LENAR»? 
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GO TO 800 
725 CONTINUE 

IF (NRECOR? .GT. NCALRECI GO TO 900 
GO TO 150 

61 750 NTRIG = NTR1G * 1 

IF ( NPTR1G<NTRIG> ,GT. 0) GO TO 600 
IF I NRECOR? .GT. NCALRECI 00 TO 900 
GO TO 100 

800 IF ( YL ASER ( 1 ♦ 1 1 .GT. YLASEK(I)I IM *1 
65 CALVELO = CALVELO ♦ I - NPTR Ifi (NTR lOl 

NSAMPLE * nsahple ♦ 1 
NTRIG = NTRIG ♦ 1 

IF (NRECOR? .GT. NCALREC) GO TO 900 
IF (NEXT .GT • 0) GO TO ISO 

70 IF <NPTRI6(NTR|GI .GT. 0) GO TO 600 

GO TO 100 

900 CALVELO = CALVELO/NSAMPLE 
npts*p = NPTSWP / NSvPS 
*RITElb.5)NRECQ&2»CALVEL0tNPTS*P 

7*; 5 FokvaMIHI .5X*0EVIAT[ON FREQUENCY CAlIBRATION*SX*NUM0ER OF RECORDS 

. USED FOR CALI6RATI0NM3/95X*wAvEL£NGTH X DEV. FREQ.*/ I OX* AVER AGE 
.NUMBER OF POINTS TO DEVIATION FREQuENCY*SX*AVERAGE NUMBER OF POINT 

.S/SWP*5X« */ 95X*2 X POINTS TO OEV. FREQ.*/ 

. 3?XF5. I * 3RX I 3 > 

80 CALVELO = (WAVLEN*DEVFREQI/(2*CALVEL0) 

VRITE(6»6)CALVEL0 
6 format (1h*.104xfs. 

RETURN 

End 
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SUBROUTINE NOISCAL 

CONHON/0NOlSCA/FLTBACK.NO!SRECtXNLEVEU275) 

COMMON/BLOCK2/LENAR«2»SYNCU5oO) tTLA$ER( 15001 «NRECOR2*NF ILE2« 
„ ZEROTM2. MRI0AT2.NTAPE2 

CONMON/BLASER/NFL YBAC.NPTSwP 
DIMENSION SUHPTS ( 300.) 

DO 100 1*1.NPTS*P 
100 SUMPTSU) - 0-0 
150 CALL BUFLAS2 
K = I 

JCOUNT = 0 

175 DO 200 I=K*LENAR«2 

IF (YLASERIU -6T. FLYBACK) GO TO 300 
200 CONTINUE 

IF tNWECOft? -GE. NOISREC I GO TO 800 
GO TO 150 

300 00 400 J-1.LENARM2 
JCOUNT = JCOUNT ♦ 1 

if < vlaseh ( ji ,lt. flyback ) Go to 450 

400 CONTINUE 

IF (NRFC0R2 ,GE, NOISREC) GO TO 800 
CALL BUFLAS2 
I = I 
GO TO 300 

450 IF (JCOUNT ,GT. IS) GO TO 500 
JCOUNT = 0 
K=J 

IF IK.LE. LENARR2) GO TO 175 
IF (NREC0H2 .GE, NOISREC) GO TO 800 
GO TO ISO 
500 M=1 

550 ISTART = J - 1 

DO 600 K = I$T ART *LENARR2 
SUMPTS(M) = SUMPTS(M) ♦ YLASER(K) 

M-M ♦ 1 

IF (M ,GT. NPTS WP - NFLYBAC) GO TO 700 
600 CONTINUE 

IF {NRECOR2 .GE. NOISREC! GO TO 000 
CALL BUFLAS2 
J= 1 

GO TO 550 

700 NSAMPLE = NSANPLE * 1 
JCOUNT " 0 

IF (K *LE. LENARR2) GO TO 175 
IF (NREC0R2 ,Lf. NOISREC) GO TO 150 
600 LAST = NPTSWP - NFLYBAC 
DO 900 l = ltLAST 

900 XNLEVELU) = SUMPTSII) / NSAMpLE 
WRITE (G»l) tXNLEVEL(I) ,1=1*LAST» 

1 FORMAT UH0,/*« NOISE LEVELS*,/. llX*lS|Ffl.3))l 
00 IOOO I=9*LAST 
1000 XNLEVEL 1 1 ) = 40.0 
RETURN 
End 
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SUBROUTINE BUFLASl 

COMMON/BBUFLAl/NTOTFIl t JCL0CK1 ♦ IEXlT tUMAOAtNREC 
COHMON/BLOCK1/LENARR1«WINDIRE(IOO) tNRECORl tNFILEl * 

. ZER0TN1*0IRM1rR< 100 > « VOLT <2*1001 »*R10AU 

C0 M M0N/8SKPE0l /LPACOAl ♦ IDEnT I tNFLSKPl *NRCSKPl 

CQMMQN/UNRK 1 / 1 T IMEI t ICOMWRU ( 20l 1 * IUATWR011000I 


CORTlMl = 3H NO 

100 BUFFER INU*1> ( IT IMEI . ICOMdRO <LPACDA1I > 

IF (UNITUH SQO * 200 *400 
200 hR1TE(6*1) nRECORI *NFILE 1 
1 FORMATdrtQ,* ThEWt AREMS* RECORDS ON FILEM3* UNIT 
,D IN BUFLA5!*) 

NREC0R1 = 0 
NFILEl - NFILE1 ♦ 1 
IF InFILEI .6T. NTOTFin 00 To 300 
If (10ENT1 «EQ* 3MYES) CALL HFADERl 
60 TO 100 
300 IEXIT= 3HYES 
RETURN 

400 NRECOW 1 = NRECOR] ♦ l 
LEN * LENGTH { 1 ) 
teRITE ( 6 ♦ 2 > NRECORi^NFILeitLEN 
2 FORMAT OHO** PARITY ERROR ON NEXT DATA* R£C0R0 # I5* 
# MBER OF COMPUTER *0RDS*l4l 
1F1LEN «NE * LPACUAl) 60 TO 100 
CALL UNPAK1 
CALL SORT l 

IF ( 4R I L>AT 1 «£Q* 3H NO) CALL DATWRIl 
60 TO 600 


l. encounter* 


F!LE*l4*5X* NU 


SOO NREC0R1 = NRECOP 1 ♦ 1 
LEN = LENGTH m 

IF a*N .NE. LPACDA1) 60 TO 7oO 
CALL UNPAK1 


CALL SORT! 

600 IF (NRECOHl .EQ. 1) ZEROTMl » ITIMU 

IF l 1TIMEI-999999 * GT * -12000 ) C0«TlHl«3HYES 

ITImFI = I T IME 1 ♦ JCL0CK1*99999R ♦ T 1HADA/ INREC-U • INRECORl-l I 

IF (C0RTIM1 # EQ« 3HYES) JCLOCkI * jCLOCKl ♦ I 

RETURN 

700 dftITE(6#3)NREC0Rl*NF[LEl*LEN 

3 FORMAT I 1 HO * * RECORD ENCOUNTERED OF IMPROPER LENGTH ON UNIT 1* 
.RD*I4« F ILE*I2* NUMBER OF COMPUTER *0R0S*I4> 

GO TO 100 

BOO RETURN 
END 


RECO 


RAGE I 


4 
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SUBROUTINE MEADER1 

COMMON/ BLOCK 1 /LENARRi INDIRE ( 100 1 tNRECORl »NF ILEl * 

« ZEROTMl ♦ D1RNIRR I 100) tVOLT <2*1001 *WRI0AT1 

C0HM0N/BHEA01/10(91 

S 50 BUFFER 1NU*OMIO<1), ID<9)1 

IF (UNIT < 1 1 ) 300*200*100 
200 PRINT 1 *NF JlE 1 

1 FORMATUHO,* EOF READ IN HEADER ON FILE*12 # UNIT 1*1 
00 TO 50 

1ft 100 PRINT 2* NflLCl 

2 FORMAT (1H0» PARITY ERROR IN HEADER ON FILE*I2* UNIT l«» 

300 LEN = LENGTHU1 

PRINT 3 * NFILEN (I0(!).T«U2)*LEN 

3 FORMAT (1 HO • * 10 ON UnIT 1* FItE*l2* IS *2A10* NUMBER OF COMPUTER N 
„OWOS*U> 

RETURN 

ENO 
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SUBROUTINE SORT! 

COHMON/BSORT 1/ IBEGSK 1*1 SKlPltF ACTOR 1 
COMMON/ BLOCK I /LENARR 1* MI NO I RE (100) .NRECORI *NF ILE1 1 
. ZEROTH I .D1RMUR 1 100 ) • VOLT (2*100) tWRIOATl 

COMMON/UNPK 1 / IT IME1 » ICOMURD <2Ql ) v IOATMRD 1 1000) 

M— IBEGSKI 

DO 100 I = ULENARRi . 

KIND IRE <1 ) = IDATKRD <M) * f ACTORl . 

V0LT<1»1> = 1 DAT KHD t M ♦ 1 1 * FACTOR1 
DIRMIRRII) = IOATkRD (N* 2> * f ACTOR 1 
VOLTIZ.I) = ID ATKRO (M*3 > ♦ FAfTORl 
100 M= ISK1P1 

IF (KRIOAT1 *EQ. 3HYES) CALL DATKR 1 1 

RETURN 

EnO 


% 
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SUBROUTINE DATWRI1 

COMMON/ BLOCK 1/LENARR1 vWINDIREclOO) .NRECOR1 «NF ILE1 » 

. 2EROTM1 tOlRHlRR(lOO) ♦ VOLT ( ?i 1 00 ) tfcRlDATl 

COMMON/UNPKl/!TIMEl*ICOmRD(2ol ) * IDATWRD 11000) 

WRITE (6*11 NREC0R1.IT1ME1 

1 FORMAT ( 1 M J » • RECORD NUMBER *t*.6X* ITIM£1«I6) 

WRITE (6*2) ( VOLT C 1 *1 1=1 ♦LENARR1) 

2 format uho*/* ux* io (Fio.s* 2 x> n- 

WRITE (6*2) (V0LT(2*I)*I=lfLENARRi) 

WRITE (6.2) (wINUlRE(I) *I=ltLFNARRl> 

WRITE (6*2) (0IRM1RRII) ,I = l,LFNARRU 

RETURN 

END 



SUBROUTINE 


S 


10 


IS 


20 


?s 


30 


3 * 



SO 
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SUBROUTINE BUFLAS 2 

COMMON/ 0 BUFLA 2 /JCLOCK 2 *IEXTlME*NREC 2 .NREC 3 fNREC 4 fT!MADJ|*TINAOJ 2 * 

T !HA 0 j 3 *T IMA 0 J 4 

COMHON/ 8 LOCK 2 /LENARR 2 *SYNCU 500 I * YLA 5 ER f 1500 1 .NREC 0 R 2 .NFRE 2 * 
ZER 0 TM 2 * WRIDAT 2 .NTAPE 2 

COMMON/ & 5 K PED 2 /L PACO A 2 • NT 0 TFI 2 * IDENT 2 *NFLSKP 2 *NRCSKP 2 * 

NTOTApE.IEX I T 2 .NT 0 TwEC* TIMADJ 


COMMON/BSPEEO/SUM VELO t IS AMPLE* I D AT AH* • SUM VOLT *T 1 MRAT I tCHANNEl » 
0 I 6 RAT 1 »T IMECHG, VQLTCHGt MULT IHE *T JMEHW.DCSUPRE 
*FRST$P 0 *WR 1 TAPE.PRINT 0 K 
COMMON /BMW I TE/T IME 2 USGQ) * VELOL AS U 500 ) t IPO |NT 

C0MM0N/UNPk2/ITIME2*LC0mmRD<O01I*L1>ATMRD(3000J 


C0RTIM2 = 3H NO 

100 BUFFER IN<2*1 > i I T IME2.LCQMWRD <LPACl)A2) 1 
IF (UNITU)I 400*200.300 
200 WRITE <6* 1 1 nRECOH2*NFILE2*NTApE2 
1 FORMAT UHO.* THERE AWE* Ip* WFCOWDS ON 
NFILE2 = NFILE2 ♦ 1 
NREC0R2 = 0 

IF INFILE2 .6T * NTOTFI2I GO Tfl 250 
225 IF < 1 DENT 2 .EG. 3HYES) CALL Hf AOEW2 
GO TO 100 

250 NTAPE2 = NTAPE2 * 1 

IF INTARE2 .GT. NT OT APE > GO To 600 
NFILE2 = 1 
NTOTFI2 = l 
CALL UNLOAD* (21 


FIIEM 3 * TARE* 12 > 


JCLOCK 2 = 0 
IEXTIME - I 5 TORTM 
260 GO TO I 22 U. 230 * 2401 *NTAP £2 
220 NTOTREC = NWEC 2 
TIMADJ = TIMA 0 J 2 
GO TO 225 

230 NTOTREC * NREC 3 
TIMAOJ * TIMA 0 J 3 
GO TO 225 

240 NTOTREC * NREC 4 
TIMADJ = TIMA 0 J 4 
GO TO 225 

300 NRECOR 2 = NWEC 0 R 2 *1 
LEN * LENGTH ( 2 1 

WRITE 16.31 NRECOR 2 .NFILE 2 tNTAPE 2 tLEN 
3 FORMAT (lHOf* PARITY ERROR ON REC 0 RD*I 6 * FIL£* 13 * TAPCM 2 * NUMBER 
.OF COMPUTER WORDS * 14 1 
IF (LEN .NE. LPAC 0 A 2 ) GO TO 100 
CALL UNPAK 2 
CALL SORT 2 

IF IWRIOAT 2 .EG. 3 H NO) CALL OATWR 12 
GO TO 500 

400 NRECOR 2 = N«tCOR 2 ♦ 1 
LtN = LENGTH< 2 > 

IF ILEN . EQ. LPACOA 2 ) GO TO 450 
WRI TE Ib* 4 > LEN*NRECOW 2 *NF ILE 2 *NT APE 2 
4 FORMAT 1 1 HO * « ENCOUNTERED WECOPO OF IMPROPER LENGTH. LENGTH WAS* 13 * 
. COMPUTER WORDS. THIS OCCURRED ON kECORO*I 5 * FIlE*I 2 * T APE* 12 * ON 
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♦UNIT 2*1 
GO TO 100 
450 CALL UNPAK2 
CALL S0RT2 

500 IF <NTOTREC .£Q. 01 RETURN 

IF ( NRECOR? .£0. 11 ZER0TM2 * ITINE2 
IF IITIHE2 - 999999 -GT. -935j C0RTIH2 * 3MYES 

ITIME2 = ITIHE2 ♦ JCL0CK2*999999*(TIHA0J/lNT0T«EC-in*lNReC0R2-ll 
* ♦ IEXTIM1 

550 IF JCORTIH? *EQ . 3HYES1 JCL0CK2 * JCL0CK2 ♦ l 
ISTORTM = IT IME? 

IF <nrecor? .le. ntotrecireturn 

rfRlT£<6*2} NREC0R2 *NT APE2 

2 FORMAT MHD*5K*REACHE0 RECORDMS* On TAPE*I2* KITMOUT E0F*1 
GO TO ?00 

800 IF (WR 1 TAPE .NE. 3HYESI CALL EXIT 
LENARR2 ' 1 
CALL LASWNIT 
700 ENOFILE 3 
ENOFILE 3 
ENOFILE 3 
ENOFILE 3 
RE* INO 3 
CALL EXIT 
RETURN 
ENO 
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SUBROUTINE HE ADER2 

COHHON/8LOCK2/LENARR2*SYNC < 1S00 > « YlASER 1 1500 » *NRECOR2vNF ILE2* 

* ZER0TH2t*RIDAT2.NTAPE2 

COMMON/BHE AD2/ ID (9) 

50 BUFFER IN ( 2 ♦ 0 ) (10(11,10(9)1 
IF (UNIT(2)> 300*1 00,200 
100 PRINT 1 *NF I L£2*NT APE2 

1 FORMATHHO,* EOF IN HEADER ON FILEM2* TAPEM2* ON UNIT 2.*l 
GO TO 50 

200 PRINT 2* NFILE2*NTAPE2 

2 FORMATUHO,* PARITY ERROR IN HEADER ON f JLE*I2* TAPE*I2* UNIT 
300 LEN - LENGTH(2> 

PRINT 3, NFILE2* NT APE2 • (ID ( I ) ♦ I»1 »2 ) f LEN 

3 FORMATUHO,* 10 ON FILEM2* TaPE*1** UNIT 2 IS *2A10* NUMBER OF CO 
.MPUTER *ORDS*UI 

RETURN 

END 
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SUBROUTINE SOPT2 

COMMON/BSORT2/I0EGSK2*1SKIP2 uorroa? M rn PJ>* 

COMMON/BLOCK2/LENARR2*SYNCa50O> * YLASER (1500 » »NPEC0«2*NF ILE2 h 
ZEROTM 2ttfPlDAT2.NTAPE2 

*C0MM0N/UNPk2/IT IHE2*lCOM*RD <O0l > »U>ATlfRO< 30001 
«=ie£GSK2 

DO 100 I-I.LENAPR2 
SYNC IIJ = loat*ro<*> 

YLASE«m = tOATy«OIH*ll*l-l.O) 

100 M*M* 1 SK I P 2 _ 

ir HR1DAT2 .EQ. iHYES ) CALL 0 ATMR 12 

RETURN 

END 



SUBROUTINE 



0ATWRI2 TRACE CDC 6*00 FIN V3.0-P26J OPT*0 02/10/72 13*01,02, 

SUBROUTINE DATWR 12 

COMMON/8LOCK2/LENARR2tSYNCU5nO) • YLASER 11500) •NREC0R2*NF IL£2« 

* ZER0TM2 *WR 10AT2*NT APt2 

C0MM0N/UNPK2/ IT IME2»LC0HWR0 < 60 1> tLDATWRO I 3000 > 

WRITE (6* i > NREC0R2 * IT IHE2 

1 FORMAT ( 1 H| ♦ • RECORD NUMBER* 14* IT1ME2=*16> 

WRITE <6. 2) (SyNC<nfI-l*LENApR2) 

2 FORMAT (lH0«/»<lAf 10(F10.S*U) M 
WRITE I6«2i <yuASERCI)*I=ltLENARR2> 

RETURN 

END 
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SUBROUTINE SKPE0F1 

COMHON/BLOCKl/LENARRltWlNDIREf 1001 #NRECORl*NFItEl* 

# Z£ROTMl»DIftMIfiR(IOO>*VOLT<2il001*WRlDATl 

COHMON/UNPK 1/ IT IME 1 * ICQMWRD <2ol 1 * IOATWRD (10001 
COMMON/BSKPEOl/LPACOAl* IOEnT l.NFLSKP 1 *NRCSKP1 

nrec = o 

IF (nFLSKPI .LE« ffl GO TO 500 
NFltSKP = 1 

100 BUFFER IN 1 1 « 1 1 (ITlMEl* lCOMWWD (LPACOAl I 1 
IF 4UNITC1I) 300 **+00. 200 
200 UEN - LENGTH(l) 

NREC = NREC ♦ 1 
NRECOR1 - NRtCORl ♦ 1 

WRITE <6 >2) NRECOR1*NFILE1pNFl5KPUnREC*LEN 
2 FORMAT UNO** PARITY ERROR in ftECGRD*14* FILE*I2* on UNIT 1, ENCOU 
.NTERED while skipping filemj/ 5X* number records skippeo*I3** 

. LENGTH OF RECORD* !<#* COMPUTER WORDS* I 
GO TO 1 00 

300 LtN = LENGTH m 

nrec = nrec ♦ l 

NRECOR 1 = NRECORl * 1 

IF <LEN .NF. LPACOAI) WR ItE (6* 3>lPACDA 1 *LEN<NRECORl tNFlLEl* 

, NREC 

3 FORMAT UHO,* A RECORD WAS ENCOUNTERED WITH LENGTH NOT EQUAL T0*14* 

. COMPUTER WORDS. LENGTH WAS* 14* .*/5X* «EC0R0*I4* FILE*I2* ON UNIT 
.1, NUMBER OF RECORDS SkIPPEDMa) 

GO TO 100 

400 WRITE (6*4lNR£C0Rl tNF ILE l * NREC * NF ILSKP*NFLSKPl 

4 FORMAT 11H0,5X*THEHE WEREMS* RECORuS ON FILE*I2* UNIT l.*/5XI3* RE 
.CORDS SKIPPED ON THIS FILE. TOTAL NUMBER OF FILES SKIPPEDM2* TOT 
.AL NUMBER TO BE SKIPPEDM2I 

NFILE 1= nfilei ♦ 1 

NFILSKP = NFILSKP ♦ 1 
NREC = 0 
NRECOR 1 - 0 

IF i IDENT 1 . EQ . 3HYESI CALL HEADER l 
IF (NFILSKP .le. nflskpii GO TO 100 
IF (NRCSKPI .GT. 0) GO TO 500 
RETURN 

500 00 900 1=1 tNRCSKPl 

BUFFER IN(1*11 (I TIME 1 * IC0HWRO (LPACUAl 1 1 
IF (UNlT*meOO*700*600 
600 LEN = LENGTH (1) 

NREC * NREC • 1 
NRECOR1 = NRECOR 1 ♦ l 

WRITE <<a*5> NREC0R1*NFILE1*NREC*N«CSKP1*LEN 

5 FORMAT (IMG* par ITT ERROR IN RFC0RD*14* flLE*I2* ON UNIT l.*/SX* NU 
.MBER RECOROS SKIPPED*I4* NUMBFR RECORDS To BE SKIPPEO* 13*, LENGTH 
. OF RECORD WAS* I** COMPUTER WORDS.*! 

GO TO 900 

700 WRITE (6 .6 > NRCSKPI * NREC* NRECOR l iNFILE 1 

6 FORMAT ( IhO * * EOF RE A(j WHILE TRYING TO SKIP*I3* RECORDS**!4* RECORDS 
.s HAVE been SKIPPED. RECORD nUMBEP»14* F1LEM2* ON UNIT 1*1 

GO TO 900 

800 NREC = NREC * 1 
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NRECOR I* NRECOP1 ♦ 1 
LEN = LENGTH! 1) 

IF <LEN .NE. LP ACDA1 > WftlTE!6»3>LPACDAltL€N*NREC0Rl.NFILEU 
. NREC 

900 CONTINUE 

WRITE (6* 7) NREC^NRCSKPJ * NRECOR 1 *NF ILEI 
7 FORMAT IlHO** COMPLETED SKIPPING* 14* RECORDS* NUMBER OF RECORDS TO 
, HAVE BEEN SKIPPED* I4/5X* RECORD NUHBERM5* FlLE*I2* ON UNIT l«l 
RETURN 

end 
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SKPF0F2 TRACE COC 6400 FTN V3.0-P26I OPT«0 62/10/72 13.01 *02. PAGE 1 

SUBROUTINE SKPE0F2 

CO M «ON/BSKPE02/LPACO*2# NT0TFI2* 1DENT2 *NFLSKP2*NRCSKP2* 

. NTOTAPE *IEXlT2.NT0TftEC*T IMAOj 

C0MM0N/B8UFLA2/ JCL0CK2. IEXT 1M£*NREC2 *NREC3»NR£C4*T 1MADJ1 .T IMAD J2» 

. TIMA0J3*TIMADJ4 

COMMON/BLOCK2/LENAR«2*SYNCU5oO> * YlASER 11500) f NREC0R2tNF ILE2. 

. ZER0TM2* WR 1DAT2.NT APE2 

COMHON/UNPK2/ITlME2.LCQMWftD(6on 'LDATWRD < 3000) 

NREC = 0 

IF (NFLSKP? .LE. 01 GO TO 600 
NFItSKP ^ 1 

100 BUFFER IN(?.1> U T IHEZ.LCOMWRD (LPAC0A2 ) ) 

IF (UNIT (2 1 ) 300*400.200 
200 NREC * NREC ♦ 1 
LEN = LENGTHI2) 

NREC0R2 = NPEC0R2 ♦ 1 

WRITE (6*2) WF lLt2*NTAPE2*N«EC0»2*NREC.L€N 

2 FORMAT (1H0.* PARITY ERROR OCCURRED WHILE SKIPPING RECOROS ON FILE 
.NUMBERS 1 2* OF T APE* 12* UNIT 2.*/5X* THE RECORD NUMBER IS*I5»2XI3* 

.RECORDS HAVE BEEN SKIPPED. ThE RECORD LENGTH WASM4* COMPUTER WOR 
.OS*) 

GO TO 100 

300 NREC 2 NREC ♦ 1 

NRECOP? = NREC0R2 ♦ 1 
LEN 2 LENGTH{2> 

IF (LEN *NE ♦ LPAC0A2) WRITE (6 .31 LPACOA2 *L£N *NREC0R2*NF lLE2»NTAP£2* 

. NRtC 

3 FORMAT OHO,* LENGTH OF A RECORD WAS NOT EQUAL 70*14* COMPUTER WORO 

.S. IT CONT A INEu* 1 4* COMPUTER WORDS. */SX* THIS OCCURREO WHEN RECOR 
-D*I5* WAS SKIPPED ON FILE* 12* TAPE*|2* UNIT 2. TOTAL NUMBER 

. OF RECORDS SKIPPED* 131 

GO TO 100 

400 WRITE (6*41 NREC0R2.NFILE2*NTAPE2*NRECiNFILSKP*NFLSKP2 

4 FORmaT(1hO.SX*ThEHE WERE* IS* RECORDS ON FILEM2* TAPE*12* UNIT 2.*/ 

./SAIJ* RECORDS SKIPPED ON THIS FILE . TOTAL NUMBER OF FILES SKIPP 

. ED* 12* TOTAL NUMBER TO BE SKl»PEO*12> 

NFILC2 = NFILE2 * 1 • < 

NF ILSKP = NFILSKP ♦ 1 

IF (NFILE2 .LE. NT0TFI2) GO TO 500 

NTAPE2 = NT APE2 ♦ l 

IF (NTAPE2 *LE. NT0TAPE> GO Tf) 4S0 

IEXIT2 2 3h YES 

RETURN 

450 NTOTFI2 = 1 

CALL UNLOADW (2) 

NFILE 2= 1 
NRECOR 2= 0 
NREC - 0 

475 GO TO (48G.46S.490>#NTApE2 
460 NTOTREC = NREC2 
TIMADJ = TIMADJ2 
GO TO 495 

485 NTOTREC 2 nrEC3 
TJMADJ 2 TIMADJ3 
GO TO 495 


% 


subroutine 


6 o 


6S 


70 


75 


00 


as 


90 


9S 


SKPEO ft TRACE COC 6400 fTN V3.0-P261 OPT»0 02/10/72 13*01.02. PAGE * 

490 NTOTREC = NREC4 
T1NADJ = TIHA0J4 

495 IF UDENT2 ,£0. 3MVES1 CALL Nf A0ER2 
I r INfllSKP «LE« NFLSKP21 GO TO 100 
IF INRCSKP? .GT « 0) GO TO 600 
RETURN 

SOO NREC0R2 = 0 
NREC = 0 

IFINFILE2 *E0 • NTOTFI21GO TO 4?5 
IF < I0ENT2 .EQ* 3MYES) CALL HFADER2 
IF INFILSKP .LE • NFLSKP?» GO TO 100 
IF (NRCSKP? .GT • 0) GO TO 600 
RETURN 

GOO 00 1000 1=1 .NRCSKP2 

BUFFER IN<2*1> < IT IHE2*LC0H WRf) (LPAC0A2I 1 
IF <UNIT<2t) 900*600*700 
700 LEN = LENGTH <2 1 
NREC = NREC * 1 
NRECOR? = NREC0R2 ♦ 1 

WRlTE(b*5) NREC0R2.NFlLE2tNTApE2*NREC*L£N 
S FORMAT MHO** PARITY ERROR OCCURREO WHILE SKIPPING RECORDS* RECORD 
. NUMBER* IS* FILEM2* TAPt*I2« ON UNIT 2.*/5Xl4* RECORDS HAVE BEEN 
.SKIPPED* LENGTH OF RECORD WAS*I4* COMPUTER WORDS**! 

GO TO 1000 

800 WRITE 5 6*61 NREC *NREC0R2 *NF ILE?*NT APE2 

6 FORMAT (1H0* AN EOF WAS ENCOUNTERED WHILE SKIPPING RECORDS** 15* RE 
.CORDS HAVE BEEN SK IPPEO .*/$X* RECOkO NUMBER*1S* OF FIl£*I2* ON TAP 
.E* 12* OF UNIT 2.*> 

GO TO 1000 
900 NREC = NREC ♦ 1 

NRECOR2 = NREC0R2 ♦ 1 
LEN = LENGf H < 2 1 

IF (LEN .HE. LPACDA2) WW HE (6 *3 >LPACUA2t LEN *NREC *NFILE2*NT APE2* 

„ NREC0R2 

1000 CONTINUE 

WRITE <0* 7) NREC *NRCSkP2 * NRECOR 2* NF ILE2 *NT APE2 

7 FORMAT HHOt* COMPLETED SKIPPING* {4* RECORDS* NUMBER OF RECORDS TO < 

. HAVE BEEN SKlPP£0*l4* RECORD NUMBER*15* FILE*I2* TAP£*I2* ON UN 

.IT 2* 1 
RETURN 
ENO 
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SUBROUTINE CONSTHR 

COHNON/6CONSTM/NAVEMIR*DlfttCMR»CH6MlR*TlMEM!R 
COMMON/ BLOCK 1/LENARR1 »W INOIRE < 1 00) .NREC0R1 tNF 1LE1 • 

. ?ER0TM1>DIRM1PR< 100) 'VOU<2')OOI'WRIOAT1 

COMMON/BCAL IBR/SLOPE (2 ) *ZER0TaP< 2) * SLOPE AN* AN INTER* SLOPEHWa 
. SLOPEWD* WD INTER *SLOPEHD#DM INTER 

AVEMIft = 0 
NREC0R1 = 0 
CALL 0UFLAS1 
100 CALL 6UFLAS1 

DO ?00 K*1,L€NARR1 
200 AVEMIR = AVEMIR * DlRMIpRlK) 

IT (NRECOR1 .LE. NAVEHIRI GO TO 100 
OIRECMR = AVEMIR/<LENARRl*NAYEM!RI 
LAST= NAVEnIR ♦ 1 
DO 300 1=1 tL AST 
BACKSPACE 1 
300 CONTINUE 

REAR {S* 1 ) CHGMlRtTlHEHlR 

1 FORMAT <A3*F6.?> 

*RITE(6*2) NRECOR1 • DIRECMRt SLOPERDtUH INTER 

2 FORMAT UMO*5A*MIRROR DIRECT I0n*5X*NUMBER OF RECORDS USEO FOR AVER* 
* GE # 1 3/ 10 A* AVER AGE VOLT AGE*SX#SlOPE«5X*TNTERCEPT*5A*OIRECT ION* OEGR 
.£E5*/UAF7*3*9XF7.3,5XF5 # 3> 

01 RE CHR = SL0PEH0*D1RECMR ♦ DmJnTER ♦ 180 
WRITE (6*3)01 RECmR 

3 FORMAT I1H*«5^AF7.3) 

NRECOR 1 = 0 
RETURN 

END 
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SUBROUTINE VOCTADJ 
COHHON/BVOLTAO/ISCALE 

COMMON/BCALIBR/SLOPE<21.ZEROTaP<2)*$LOP£An*ANINTER»SLOPEHW* 

. SLOPE WO ? WD INTER »SL0P£MD* DM INTER 

GO TO (100*200*300) » I SC ALE 
100 SLOPEAN = 47.736 
ANINTER = -0*4H 

WRITE 16* 1 ) 1 SCALE ♦ SLOPEAN* AN INTER 

1 FORMAT ( 1 H0*5x*ANEM0METER VALUfS*5X* SCALE* 12/ l 0X*SLOPE*5X*INTE«C£P 
*T*/10XF6.3*6XF6.3> 

RETURN 

200 SLOPEAN = 43*021 
ANINTER = 0.451 

wr ite (6* i > i scale • slope an* an inter 

RETURN 

300 SLOPEAN = 0.0 
ANINTER = 0.0 

WR1tE(6*1MSCALE*SL0PEAN*ANINtER 

RETURN 

ENO 
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SPEED TRACE CDC 6400 FTN V3.0-P261 OPT-Q 02/10/72 13.01*82. 

SUBROUTINE SPEED 

COHHON/BSPEED/SUHVELO*ISAHPLE« IDATAHWtSUMVOLTfT INRAT 1.CHANNE1* 

. DIGRAT i tTIHECHG* VOCTCHGfHULT 1ME*TIMEMW*DCSUPRE 

, ♦F«STSPD*WPITAPE*PR1NT0K 

COMMON/BBUhPUP/TIMEI (70?) *VELOC2<702> 

COMMON/0LOCK1/LENARR1 **1N01RE ( 1 00) *NREC0R1 *NF ILE 1 * 

, ZEROTH 1 *0IRM JRR (100) *VOLT (2*1001 «*R1DAT1 

COMMON/ 8CAL I HR/SLOPE (2) *Z£R0TAP(?> *SLOPEAN« AN INTER *SLOP£Hy » 

SLUPEwDiwD INTER* 5L0PEMD* DM INTER 
COMMON/UNPK1/ITIME1 * ICOMHRU<2ol > *lUATWRDU0001 
COMMON /B8UMP/ HR I T IME ( l DO) • VELnC 1(1001 
COMMON/BBUFLAl/NTOTFIlt JCLOCK] *IEXIItTIMAOAiNREC 
PR INT OlC = 3M NO 
KFIPST = 1 

IF (FM5TSP0 .EO. 3HYE5 .A* NRfCORI -NE. 0) GO TO 2S 
CALL BUFLAS1 

IF (IEXlT .EO. 3HYES) GO TO 2fl0 
?S DO 100 JDATAWS = KF IRST .LENARP 1 

VEL0C1 (IDATAHS) =■ SL0PE(1>* VOLT < 1 » 1QAT AHS) ♦ ZEROTAP(l) 

VEL0C1 < IDATAHS) a (SLOPEAN* Vf LOCI ( TOATANS ) ♦ ANlNTER) *0*3048 

VEL0C2 < IOATAHH 1 = SLOPE (2 ) * V0»_ T ( 2* iOAT AHS » ♦ Z£R0TAP<2I* DCSUPRE 

IF (FRSTSPO .EO. 3HYES) GO TO 50 

5UMVELO * SUMVELO * VELOC 1 I IOA T Ah'S 1 

I5AMPLE = ISAMPLE ♦ 1 

SUMVOLT = SUMVOLT * VELOC2 (1 0 aT AH* | 

50 TIME! t IOATAHHI = T IMP ATI * U IT I MU-2ERQTM1 I /1 0000 ♦ U10ATAHS 1 
.•CHANNE 1 )/OIb«ATl I 

IT <T 1ME1 UDATAHH) .GE. TIMECnG .A. VOLTCHG *EQ. 3HYES) 

. CALL voltadj 

IF < T ImEI U JATAH*) *OE. MUL T ImE*T IM£H« ) GO TO 200 
100 IOAT AH J = IDAIAHH ♦ 1 
FRSTSPO = 3H NO 

istart = ioatamh - lenapri first- i 

00 115 I =KF IRST fLENARftl 
WRITIHE(I) = TIMEKISTARTI 
115 ISTART = ISTART ♦ 1 

1?0 IF (NRECOR1 *GE. 2 -A. NRECOHl .LE. B> GO TO 125 
J = NRECOR1 - I 

IF (MOD (1*30) *EQ* 0 1 GO TO 125 
GO TO 160 

125 PRINT 3.NRECORI 

3 FORMAT Uho** ANEMOMETER VELOCITY RECORD NUMBER* 14//2X*T IMC*SEC 
# S*5X* VEtOCITT*M/SEC*10X* T IMF *SECS*SX* VELOCI T Y #H/SEC* 1 OX* TIMERS 
« ECS*5X* VELOClTY*M/SEC*/> 

PRINTOK a 3HYES 

HR I TE (6*1) (wRITIMEm *VEL0C1 ( ! * ♦ 1 = 1 *LENARR1 1 
i FORMAT (1H * lxF0.3*llXF6.3*15XF8*3*nxF6.3*15XF0*3*|lXF6.3) 

IF fuRJTAPE *NE. 3HYES1 RETURN 
160 WRITE(6*12) NRECOR1 
12 FORMAT i 1M *• NRECOR 1 -* 15) 

BUFFER OUT(3*ll (riklTlMEU) tVfLOCl (100)1 
IF (UNIT < 3 J 1 400*180*1^0 
180 WR1TE16*61 NRECOM 1 * NF ILE 1 
6 FORMAT ( 1H ,* EOF ON RECORD NUmBERMIO* FILE NUMBER*^) 

GO TO 400 
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60 


6^ 


70 


7*> 


8n 


8* 


9A 


99 


100 


109 


110 


190 WRITE (6*71 NRECOR 1 • NF ILE 1 

7 FORMAT (1H *• PARITY ERROR INPUT ON RECORD NUMBER* MO* FILE NUMBER 
.*151 
GO TO 400 

200 AVEVELO * 1 SUMVELO/ ISAMPLEl / 0.3048 
AVEVOLT * SUMVOLT/ ISAMPLE 
KFIRST “ 1DATAWS ♦ 1 
DO 250 JK = 1 * I DAT AW5 

250 WRITIMEUK) - TIMEKIOATAHW-IDATAWS ♦ JK| 

SUMVELO = 0.0 
SUMVOLT = 0,0 
ISAMPLE = 0 
MULTIME = MULTIHE ♦ 1 

HWlNTER = SORT ( AVEVELO >-SL0 PEhW*AV£VQLT**2 
DO 300 1*1*1 DAT AM* 

300 VEL0C2 < I ) * ( < SL0PEMW*VEL0C2 ( I ) **2*hW IN TER I •• 21*0.3048 
IF (NRECOR I ,GE. 2 .A. NRECORl *LE. 8) GO TO 305 
I = NRECORl - 1 

IF (M0D( I *30) *EU* 0 ) GO TO 305 
GO TO 3|0 

305 WRITE(6*2> TlMEl(l>*TIMEl(2DATAHW)»NRECOR|t<T!HEmi'VELOC2(I>» 1» 
. 1 + !OAT AMW1 

2 FORMAT MHO** HOT WIRE VELOCITY* CALCULATED FOR TIME PERIOD FROM* 
,F7 , ?* T0*F 7*2* RECORD NUMBERM3//* TIME*SECS*5X* VELOCI TVtM/SEC*10 
.X* TIME.S£CS*5X* VELOCITT*M/SEC*10X* T IME*SECS*5X* VELOOTV*M/SEC* 
.//* MX.F8.3* 12X*F6,3*l<*X*Fe.3,12X*F6.3*14X*F8.3*12X*F6,3>> 

IF (WRITAPE ,EQ. 3HYES) GO TO 310 
IF MEXlt .EQ. 3HYESIRETURN 
307 IDATAHW = i 

IF (KFIRST ,GT, LENARRUGO TO 120 
GO TO 25 
310 N * 1 

PRINT 13*NREC0R1 • IDATAHW 

13 FORMAT ( 1 H •* H, W* NRECOK 1 tt * 1 1 0* NUMBER OF WOHOS * *M0| 

LAST = 2 

IF (IDATAHW .LE» 30DLAST * 1 
00 370 !=1*LAST 
N - M * 300 

IF (I .EQ. LAST) N = IOATAHW 
BUFFER OUT (3*1) (T IHE1 (H) *T iMEl (N) 1 
IF (UNIT (31 1 340*330*320 
320 WRITE(6*S) NREC0R1»NFILE1*M*N 

6 FORMAT ( lH ** PARITY ERROR ON HW TIME* R€CORO*UO* FILE # I3* M«*|5* 
*N=* IS) 

GO TO 140 

330 WRITE16.9) NRECORl *NFILE1 •M.N 
9 FORMAT MH ** EOF ON HW TIME* »ECORU*I10* FILE«13* M**I5* N**ISl 
340 BUFFER OUTO.l) (VELOC2 <M> ♦ V£lOC2(N) ) 

IF (UNJT(3l) 370*360.350 
350 wft I TE (6* ! 0 ) NRECORl *Nf I lEI *m*n 

10 FORMAT <1H ** PARITY ERROR ON HW VELOCITY* RECOROMIO* F1LE*I3* H*« 
.15* N-*I5> 

GO TO 370 

360 WRITE (6*11) NRECORl »NF1LE1*M*N 

11 FORMAT MH *• EOF ON Hw VELOCITY* RECORD* 110* FlLE*13* M»*1S* N«*I 



ns 
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.5) 

370 M * M ♦ 301 
fOATAHw = 1 

IF (IEXIT .EO. 3HVESI RETURN 
IF (KFIRST .GT. LENARRl )G0 TO 120 
GO TO 25 
*00 RETURN 
END 
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AVEWINO TRACE COC 6*0» FTN V3.0-P261 OPT«0 02/10/72 13,01.02. PAGE t 

SUBROUTINE AVEtfINO 

COMMON/BAVEWIN/JSAMPLE •SUMWlNDtMULT INI tTlHAVnD,AVE*0 <3000 1 
. »LASTIME 

C0HM0N/8SPEED/SUMVELQ* I SAMPLE • IDAT AHW tSUMVOLT » T IHRAT 1 f CHANNE1 • 

. DIOR AT 1 *TlNECHG»VOLTCHGtHUtTlME«T JHEHWf DCSUPRE 

• ♦ FRSTSPQtWR IT APE *PR]NTOK 

COMHON/BBUMPUP/T IME 1 <702) *V£LnC2<702) 

COMMON /BLOCK 1/LENARH 1 1 W I NO I R£ f 1 00 ) ,NREC0R1 *NF iLEl « 

. ZCROTH1*DIRHJhW<100) . VOLT <2*1001* WRIDATI 

COMMON/BCAL1BR/SLOPE<2) *2EKOTaP(2> » SLOPE AN » AN INTER »slopehw« 

. SLOPE MO t fcD IN TE P *5L0pEHD*MD INTER 

J= IDATAHW - LENARR1 - 1 
DO 100 k=i,lenarri 
J - J * 1 

JSAMPLE - JSAMPLE ♦ 1 
SUM W 1 NO s SUMWINO ♦ W IND IRE (K I 
IF (K .LT* LENARRllGO TO 100 
150 A VEwO (MULT IH 1 > = $um*ind/jsample 

A VEWQ (MULT IM 1 ) = SLOPEWD*AV£tir)(MULTlMl) ♦ POINTER 
MULTIM1 = MULT1M1 ♦ 1 
5UMW1NO * 0*0 
JSAMPLE = 0 
100 CONTINUE 

IF { rfft I TAPE *EO, 3HTESI RETURN 

«RITE <6«1) T IMAVMDt MULT 1M1* ( A VEtfD ( 1 1 » I*LAST IME* MULT IMl I 
1 FORMAT <1H ,* MINO DIRECTIONS IN DECREES* MEANS OF DATA FOR*F5.1* 

.SEC INTERVALS. THIS IS INTERVAL NUMBER* I A/ <1 OF 10. 31 ) 

LASTIME = MULTIM1 

RETURN 

END 
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SUBROUTINE LASWRIT 

COMMON TIME 2(200) ♦ VELOLAS (200) * I POINT 

C0MM0N/BL0CK2/LENARR2*SYNCC1500) *YLASERU500) *NREC0R2*NFILE2» 

* ZER0TM2#WRIDAT2,NTAPE2 

WRITE (6*5) NRECOH2.IPOINT 

5 FORMAT ( 1 H .* NREC0R2-* 110* NUMBER OF WORDS **I1D» 

50 M = 1 

N = IPO J NT - 1 

BUFFER OUT (3*1) ( T JME2 ( M ) ♦ T IMF2 (N > > 

IF I UN IT ( J > 1 300*200*100 
100 WR ITE ( 6* 1 ) NRECO«2*NF!LE2tM*N 

1 FORMAT (1H .* PAR1TT ERROR ON LASER TIME RECORD NUMBER*! 1 0*F 1LE NUM 
.BER* 13* M=*IS« N = * 1 5 ) 

GO TO 300 

200 WRITE (6*2) NWECOR2*Nf ILE2*M,N 

2 FORMAT ( 1H ,* EOF ON LASER TIME RECGROMIO* FILE *13* Mc*I5* N**I5> 

300 BUFFER 01JT(3*n ( VELOL A$ IM > * VFLOLAS <N) > 

IF (UNIT (3) )60Q*SQQ*^00 
400 WR1T£(6*3> NREC0R2iNFILE2*M*N 

3 FORMATUM .* PARITY ERROR ON LASER VELOCITY* REC0RDM10* FILE*I3* 

*M=»IS* N=*I5) 

GO TO 600 

SOO WR ITE (6 *4 ) NRECOR2 *NF ILE2 *M *N 

4 FORMAT (IH ,• EOF ON LASER VELOCITY* RECORD* 11 0*F ILEM3* M**I5* N=* 

.IS) 

600 TIME2U) * T I H£2 ( I POINT ) 

VELOLASU) = VELOLAS(IPOINT) 

RETURN 

END 
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I dent unpaki 

INSERT LENGTHS OF PACKED AND UNPACKED ARRAYS 




312 

lengtha 

SET 

202 





17S0 

LENGTHS 

SET 

1000 



5 



* 









USE 

/UNPK1/ 




0 


NE 

BSS 

LENGThA' 




312 


0 

BSS 

lengths 







USE 




10 




ENTRY UNPAK1 




0 


UNPAKI 

BSS 

1 




1 

7170000001 


SX7 

IB 





7100004000 


SXO 

4000B 




2 

43214 


MX2 

14? 



15 


S110000000 c 


5A1 

nE 

AU) 

— FIRST word of arr< 


3 

6160000311 C 


SBb 

P-1 

B1J) 

BASE 



6170002261 C 


SB7 

^♦LENbTHB-i 

B(LAST> 


4 

6150000060 


SB5 

40 





6110000074 


SB | 

^0 



20 

5 

5011000001 

GET60 

$A1 

h 1*1 

GET 

Am 


6 

6166000001 

GET12 

SBb 

l 





11621 


0X6 

X2*X1 

MASK 

OUT 12 BITS 



67515 


SBB 

Ml-05 


RIGHT SHIFT 


7 

22656 


LX6 

M5.X6 


BUT 

25 


67515 

* 

SB5 

01-05 


AVOID SIGN EXTENSION 



11760 


0X7 

X6*X0 


CK FOR SIGN BIT 


10 

0307000011 ♦ 


ZR 

X7.STORB 





15660 


BX6 

-X0*X6 


MASK OUT SIGN BIT 

30 


14666 


BX6 

-X6 




1 ! 

11 21601 

56660 

0 <*67000000 

12 0550000014 ♦ 

6150000060 

13 * 321 * 

0400000005 ♦ 

14 20260 

6155777763 

15 0400000006 * 

0 

16 


DELETE 2ERQ-0IT RIGHT-FILL 
STORE IN BU> 


EQ 06.B7.DONE 
NE R5.0O. 1NMIO 
SB5 40 
MX2 12 
EQ f?E T 60 

LX2 46 
$05 H5-12 

EO GET 12 
EOU ijNPAK 1 
END 


STORAGE USED 
6400 ASSEMBLY 


44 STATEMENTS 
0.336 SECONOS 


10 SYMBOLS 
23 REFERENCES 





UNPAK2 


COMPASS - VCR 2.0 M 02/10/72 13.01.S7* 


1132 

5670 


LENGTHA 

lengths 


IOENT UNPAK2 

INSERT lengths of packed and unpacked arrays 
SET *02 
SET 3000 


0 

1132 


10 


IS 


20 


25 


30 


35 


4 0 


NE 

6 


UNPAK2 


/UNPK2/ 

LENGThA 

LENGTmB 


USE 
BSS 
B5S 
USE 

ENTRY UnPAkZ 
8SS 1 


1 

7170000001 


SX7 

IB 



7IOQ004DO0 


SAQ 

40008 


2 

43214 


HA2 

12 



5110000000 C 


SA1 

nE 

A(l) FIRST WORD OF ARR( 

3 

6160001131 C 


Sd6 

e-i 

0UI BASE 


6170007021 C 


5d7 

P*LENGTHB-1 B (LAST 1 

A 

6150000060 


SB5 

48 



61 1000007* 


$B1 

60 

GET A ( I 1 

5 

501 1000001 

GET60 

SA1 

AIM 

6 

616600000] 

GET 12 

5B6 

R6*l 



11621 


8X6 

*2*X1 

MASK OUT 12 BITS 


67515 


SB5 

B 1 *B5 

RIGHT SHIFT 

7 

22656 


1X6 

R5.X6 

BUT 


67515 


585 

RI-B5 

AVOID SIGN EXTENSION 


11760 


BX7 

X6*X0 

CK FOR SIGN SIT 

10 

0307000011 ♦ 


ZR 

X7.STORB 



15660 


8X6 

-X0*X6 

MASK OUT SIGN BIT 


14666 


BX6 

-A6 


11 


STOR0 

BSS 

0 


11 

21601 


AA6 

1 

DELETE ZERO-BIT RIGHT-FILL 


56660 


SAb 

R6 

STORE IN BfJ> 


0467000000 * 


CO 

H6.tt7.D0NE 

12 

0550000014 ♦ 


NE 

RS.SOt INM10 


6150000060 


SB5 

46 


13 

43214 


MX2 

12 



0400000005 * 


EO 

GET60 


1* 

20260 

INM10 

LX2 

48 



6155777763 


SB5 

P5-12 


IS 

0400000006 ♦ 


EQ 

GET 12 



0 ♦ 

DONE 

ECU 

IJNPAK2 


16 



END 




63417 

STORAGE 

USLO 


44 STATEMENTS 10 ' 



6400 ASSEMBLY 


0.339 SECONDS 23 1 
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UNLOAD* 






IDENT 

ENTRY 

USE 


2*6 

03111720000000000000 

CIOC 

VFD 

s 

2*7 

15230 70 00 000 00 000 000 

MSGC 

VFO 


250 

01022401000000000000 

ABTGS 

VFD 


251 

01022400000000000000 

ABTC 

WFD 


252 

1505 150000000 0 QQOOOO 

MEMC 

vfd 


253 

22031*00000000000000 

RCLC 

VFO 

10 

254 

05260400000000000000 

ENDC 

VFO 


255 

01030520000000000256 ♦ 

CNTCC 

VFO 


256 

oooooooooooooooooooo 

CNTCCB 

DATA 


257 

06111405552701235522 

MSG1 

OATA 


263 

S52516141701045517J6 

MSG2 

DATA 

IS 

266 

552516141701 0*551716 

MS63 

DATA 


271 

55202217072201 155503 

MSG4 

DATA 


21 4 

SS222516S50102172224 

MSGS 

DATA 


300 

55222516550516241 122 

MSG6 

DATA 


304 

55222516550102172224 

MSG7 

qata 

20 

307 

5527011 1241116075506 

MSG9 

DATA 


314 

552225 1 655051 6040504 

MSG8 

DATA 


321 

OOOOOOOOOOOOOOOOOOOO 

5 AWE W 

DATA 

USE 



CALLPP MACRO 
1FC 
BX7 

♦ SA5 
NZ 

♦ SA7 

♦ SA5 
NZ 

ENDM 


*S 3 CLOSER MACRO 

LOCAL 

SA3 

$A4 

1X5 

50 ZR 

MXU 

LXO 

SA3 

6X7 

55 ex* 

SX1 

6X3 


COMPASS - VER 2.0 H 02/10/72 13.0) .59 


PAGE 
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UNLOAD* 

UNLOAD* 

DATA. 

10/3LC 10 ♦ 2/ 1*40/0 

18/3LM5G*42/0 

18/3LABT»6/l*36/0 

18/3LA8T»42/0 

lfl/3LMEM»A2/0 

2 B/3LRCL *42/0 

2 6/ JLEND*42/0 

18/3LAC£*2/1* 22/0* 10/CNTCCB 

0 . IOB READ FORWARD 400 FOR BACK5P 

C*F ILE WAS REWOUND BEFORE RETURN * 

C* UNLOAD ON NON-TAPE FILE * 

C* UNLOAD ON UNDEFINED FILE • 

C* PROGRAM CONTINUED * 

L* RUN ABORTED WITH SPEC PROCESSING * 

C* RUN ENTIRtLV ABORTED - DUMP * 

c* Run aborted without dump * 

L* WAITING FOR NEXT REEL - GO 10 CONTINUE* 

C* RUN ENDEO WITH NO DUMP -NORMAL CC STREAM * 

0 


A 

NE***A** l 
X.A 
1 

X5** 

AS 

AS 

AS** 



IN EQ OUT 



UmlOAOW 


ZR 

X3tL0PE 

SKI 

*a 

ZR 

X3*L0PE 

SX1 

5008 

BX3 

Xl*XG 

NZ 

X3fLOPE 

sxo 

20 

BAG 

X0*XA 

SX5 

240 

BX6 

X6*XS 

0X6 

X7*X6 

SAG 

b2 

SAG 

02 

SAG 

CIOC 

bx 1 

A5*X6 

CALLPP 


ess 

0 

ENUM 



?S 3 


30 


35 



A 5 


REMIND 


MACRO 


CLOSER 


SA5 

02 

MXO 

ie 

LAD 

18 

BXb 

-X0*X5 

sxo 

28 

BX5 

XS*XO 

Sx<* 

SOB 

BXS 

XA*X5 

bag 

X6*X5 

SAG 

02 

SAG 

B2 

SA5 

CIOC 

0X7 

XS*X6 

CALLPP 


SAG 

MSG1 

SAS 

MSGC 

6X7 

A6*XS 

CALLPP 


£NOH 



WAITER MACRO 

local lop 

SXA M$09 

SA5 M5GC 

0X7 X4*X5 

CALLUP 


COMPASS - VER 2.0 M 02/10/72 13*01.59 


PAGE 
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* FILE NOT OPENED 

. WRITE MASK 

•NOT OPEN FOR WRITE 


. SPECIAL FUNCTION CODE 


. WRITE CODE 
. AOD PARITY BIT 
• ADD LOGICAL FILE NAME 
. RESET FIRST WORD FET 


• ADD FET ADR TO CALL WORO 
END INSTRUCTION IN MACRO 


•SAVE FILE NAME 


•rewind code 
•ado parity an 
•ADD FILE NAME 


GET WAIT DAVFILE MESSAGE 

►GET PP CALL WORD 

►AOD ADDRESS TO CALL WORO 



UNLOAD* 


sxo 

looooe 

SA5 

BO 

BX6 

xs^xo 

SA6 

A5 

SA2 

RCLC 

CALLPP 

2 

SA5 

•00 

BX5 

x5*X0 

NZ 

XS»LOP 

ENDM 




UNLOAD 


MACRO 


CLOSER 


SA5 

03 

MXO 

18 

LXO 

18 

BX6 

-A0*X5 

sxo 

20 

bxs 

X5*X0 

SX4 

boe 

BX5 

A4*XS 

BXo 

X6*X$ 

SA6 

82 

SX6 

02 

SAS 

CIOC 

8X7 

X5*X6 

CALLPP 


ENDM 



40 


45 


PMSG MACRO A 
SX6 A 

SA5 MSGC 

0X7 X6*X5 


CALLPP 

ENDM 


0 

1 503100000 1 

55 3 0303000003 ♦ 

53230 



LIST 

-R 

UNLOAD* 

BSSZ 

l 


SA2 

Al*l 

♦ 

ZR 

X3p*M 

- 

SA2 

X2 

* 

0X7 

X2 


3 10733 


COMPASS - VCR 2*0 M 03/10/72 13.01 *59. 
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•SAVE FILE NAME 


•UNLOAD CODE 
•AOO PARITY 0IT 
*AOD FILE NAME 


\ 



UNLOAD* 




5170000321 ♦ 


SA7 



63210 


SB2 


4 

67202 


SB2 



0100000000 X 


RJ 

5 

5 

0720000105 ♦ 


LT 



5152000001 


SA5 


6 

0325000106 * 


PL 



5132000002 


UNLOAq 


35 

7140000307 ♦ 


*AITEr 

10 

51 

5132000002 


REWIND 


104 

• 0400000000 + 


EG 


105 

7160000266 *■ 

NTDEF 

SX6 



0400000142 ♦ 


EQ 


106 


NTTAPE 

faSS 

15 

106 

5132000002 


RE* l NO 


141 

7160000263 * 


$X6 


U2 

5150000247 ♦ 

CONT. 

5 AS 



12765 


BX7 


U3 

5150000001 


CALLPP 

20 

146 

5150000321 ♦ 


SA5 



63250 


502 


147 

6130000005 


SH3 


150 

0632000151 ♦ 

♦ 

GE 



66200 

• 

SU2 

25 

151 

0220000153 ♦ 

♦ 

JP 


152 

0400000000 * 


EO 


153 

0400000161 ♦ 

JHP 

EO 


154 

0400000167 ♦ 

♦ 

EO 


155 

0400000201 ♦ 

♦ 

EQ 

30 

156 

0400000211 • 

* 

EO 

- 

157 

040000023* ♦ 

* 

EQ 


160 

0400000000 ♦ 

* 

EQ 


161 

7160000271 ♦ 

exiu 

PMSG 


166 

0400000000 ♦ 


EQ 

35 

167 

7160000274 ♦ 

EXIT2 

PMSG 


174 

5150000251 ♦ 


S AS 



10755 


call P p 


200 

oooooooooo 


PS 


201 

7160000300 ♦ 

EX1T3 

PMSG 

40 

206 

5150000260 ♦ 


SA5 


21? 

10755 


CALLPp 


oooooooooo 


PS 


213 

7160000304 ♦ 

EXIT4 

PMSG 


220 

5140000255 * 


SA4 

45 

221 

7160000010 

LOOP 

SX6 



5160000256 ♦ 


$A6 


222 

10744 


CALLPP 


226 

5150000070 


SA5 



0315000221 * 


N2 

50 

227 

5120000254 * 


SA2 



10722 


CALLPP 


233 

oooooooooo 


PS 


234 

7160000314 ♦ 

EXITS 

PMSG 


241 

5120000254 ♦ 


SA2 

55 


10722 


CALLPP 


?4S 

OOOOOOOOOO 


PS 


32? 



END 


COMPASS - VER 2.0 M 02/10/72 13*01.59. 


SAVE* 

Xl*0O 

BO-82 

-XGETBA 

B2*B0*NT0EF 

82*1 

XS.NTTAPE 


UNLOAD* 

M5G3 

CONI. 

0 

HSG2 

MSGC 

X6*XS 

SAVE* 

X5 

SB 

03.B2.*^l 

BO 

JMP*B2 

UNLOAD* 

EXIT1 

EXIT2 

EXJT3 

tXIT4 

EXITS 

UNLOAD* 

MSG** 

UNLOAD* 

MSG5 

A0TC 


MSG6 

ABTCS 

S 

MSG7 

CNTCC 

108 

CNTCCB 

4 

70B 

AStLOOP 

ENOC 

2 

MSGB 

tNOC 

2 
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APPENDIX A-2 


Computer Program for Determination of Velocity Profil 



PROGRAM 


5 


10 


IS 






30 


3^ 


AO 


45 


50 


55 


ANEvEL TRACE CDC 6400 FTN V3.0-P261 OPT*0 02/10/72 12*53.4?. PAGE 1 

PROGRAM ANEVEL < INPUT • OUTPUT * TaPES* INPUT * T APE6 3 0UTPUT *TAP£1 *FILMPLI 
COMMON T 1ME 1 101 > *VELOC (6* 1001 *NCHAnNC*LENARR *NF ILE*NR£C0R* 1£X| f 
COMMON /9CAIBR/ NCALVAl* ACTV0LT<5> »SL0PEI6) *ZEROTAP(fcl . 

. ST ANOEV <6*5)*vAR[ 

COMMON/BBUF ANE/ IOENT *IPARITY*lPACOAT #EOFMUL*NTOTF lL f PLOT 
COMMON/8SQRT/lBEG$KP*bRlTDAT* fSK IP, FACTOR 

COMMON/8VOL TAD/ 1 CHANGE* SCALE {M*ZEhOACT(61 *VOLTChG*T ihechg* 

. 1 SCALE (6) 

COMMON/8PLOTV£/ELEV<6) *SuMAVE( 6> • lSAMPLE *L ABEL A l 4) *LABELY<4l * 

. LTITLfc<4) ,HULTIM£*IUATAPT 

COMHON/BSK ipf.o/nf ILSkP*NRECSKP 
COMMON/UNPK/ 1 TIME* ICOMMRO<2QO 1 * I DAT WHO (10001 
DATA LABEL*/40H VELOCITY# H/ SEC / 

. *LABELT/40H flevation* M / 

. *LTITLE/40H VELOC [ f Y PROFILE / 

REWIND i 

READ! 5, 1 1 IOENT,wRITAP*EOFMUL*wRITDAT#W«lTPAP*PLOT*VOLTCMGf 

. LPACDAT t NCHANNC*NCA| VAL» ISK1P* 16EGSKP*L£NARR»NT0TF1L« 

. PL0TIME*AVETIME* UScALttn *1 = 1*6) *I£LEVII) *1 = 1 *6>*VARl* 

. OIGMAT* I ACT VOLT HI *!=l*S> *1 IMERAT ^CHANNEL *T INECHG. 

, iNSTCALtTApECAL 

1 fORMAT <3X.7A3*7I4,2FS.1,6I 12) */*6<F6.2> #2F5 .2 ,5 <F5. 1> *3F3. 1/2A3) 

WRITE 16 *31 I DENT* wRITAP*EOFM,)L*wRI T DA T tWR ITPAPtPLOT • VOLTCHG* 

# LPACOAT#NCMANNC*NCALVAL* ISK IP* IBEGS«P*LENARR*NTOTFIL* 

♦ PLOT IME*AVE TIME* ( I SCALE <1 1 *1*1*61 • <ElEV<! > * J =1 *6) • MARI t 

. DIGRAT* ( ACT VOLT U ) * f = l*S) # f IMERAT t CHANNEL * T 1MECHG* 

* INSTCAL* TAPECAL 

3 FORMATUHO,* IDENT = *A4* WRITAP =*A4* EOFmUL = *A4* WR1TDAT **A4* W 
.RITPAP =*A4* PLOT =*A4* VOLTChG =«a4/* LPACOAT = M4* NCHANNC = *I2* 

. NCALVAL =*12* 1SMP =« 13* ItKGSKP = *I3« LEN ARR =*|4* NIOTFIL =*12 
./* PLOTIME =*FS.i* A VET THE =*F5.1* ISCALEU THRU 6) =*6J2/* ELEV < 1 
. THRU 61 =*C>F6.2* V A R I =*F4.2# DIGHAT =*F5.l/» ACTVOLTU THRU 5) = 

. *5F5. 1 • TIMRAT =*F4.1* CHANNf | -*F4.1* T1MECHG =*F3*1/* INSTCAL =* 

,A4* TAPECAL 

IF IIOENT.EG. 3HfES)CALL HfcAOFR 
IE* IT = 3H NO 

CORTIME = 3M NO * 

JCLOCK = 0 
ZEROTIM * 0.0 
I CHANGE = 0 

factor = sortiz.i/iz.**? - i»oi 

MULT1ME = 1 
NRECOR - 0 
BAOOATA = 3H NO 
NEATPTS = 0 
NFILE = 1 
1 PAR I T Y = 0 
1 SAMPLE = 0 
00 100 1=) .NCHANNC 
100 SUMAVE < 1 1 = 0-0 

IF (TAPECAL .EQ. 3HYES) CALL CAL I0RA 
IF (INSTCAL .NE* 3HYES1GOTO 10* 

NFILSKP = ) 

NPFCSKP = 0 
CALL SKIPEOF 



PROGRAM ANEVEL TRACE 


CDC 6400 FTN V3.0-P261 OPT*0 02/10/72 12. 53*49, 


60 


65 


70 


TS 


6 * 


6 ^ 



102 CALL VOLTADJ 

103 CALL BUFANE 

IF { lEXlT .EG. 3HYESI T ImE UOaT APT- H * T INE C IOAT APT-1 > ♦ 0*1 
IF ClEXlT *£0. 3HYESI CALL PLOTVEL 

IF (NRECOR *LO. 1 .AND* IT1ME «NE. 0> ZEROTlH ■ ITIME 
IF (CITIME - 999999) .GT. -12oOO) C0RT1ME * 3MVES 
ITIME - ITIME ♦ JCLOCK * 999999 
IF (CGRT1ME .EG. SKIVES) JCLOCK * JCLOCK ♦ 1 
CORTIME = 3M NO 
IDATAPT = 1 

105 DO 110 I*1,NCHANNC 

VELOCU. IDATAPT)sSLOPEm*VELOCUtIOATAPTl ♦ ZEROTAPU1 
110 VELOC(ItlOATAPT) = ( V ELOC ( I • InAT AP T) *SCALE t II ♦ZERO ACT < II > *0.3046 
TIME ( 10 AT AP T I - T IMER AT* ( ( I T 1ME -ZE«0TIM)/1G0Q0.*( 1 1 OAT APT- 1 1 # 

. CHANNEL ) /D IGRAT \ 

IF (TIME (I HAT APT ) «GE • T IMECHG .ANU. VOLTCHG .EQ. 3HYES1 


. CALL VOLTADJ 

IDATAPT = IDATAPT ♦ l 

IF ( AVETIME*MULTIMt *LE. T IME ( IOATaPT-1 ) .AND. PLOT 
. .EG. 3HTE5) GO TO 120 

IF ( IDATAPT »LE* LENARR1 GO TO 105 
IF < wP I TRAP «EO • irtYES ) GO TO 134 
IF (rfRITAP ,EG • 3HYESI GO TO | 35 
GO TO 103 

120 00 130 1= l.NCHANNC 

130 SUMAVE(I) s SUMAVEUI ♦ VELOC ( 1 • IOaT APT-1> 

1SAMPLE = ISAMPLE * 1 

IF { TIME(IOATAPT-I) .&E. PLOT IH£*HULT IME) CALL PLOTVEL 
IF (IDATAPT .IE* LENARH) GO Tn 105 
IF (WRITPAP • EQ ♦ 3HVES 1 GO TO 134 
IF (VRITAP ,EQ # 3MYES) GO TO 1 35 
GO TO 103 

134 *R1TE(6«2> 

? FORMAT tlHl.AX* TIME, SECS* 1 OX* VELOC IT IE5»H/SEC*4X* LEVEL I*4X 

, • LEVEL 2*4X* LEVEL 3*4X* LEVEL 4«4X* LEVEL S*4X* LEVEL 6« 

. //I 

rfRITE (6,4) (TIME(J) * (VFLOC( I, J> * 1*1 f 6) ♦ l ,LENARR> 

4 FORMAT ( IH , l00<4AtFlQ.3t28X,6(F6.3,6X>/I> 

GO TO 103 

135 CONTINUE 
140 CONTINUE 

END 
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SUBROUTINE 


5 


I* 


15 


2 « 


25 


30 


35 


40 


45 


50 


55 


CAL ! BRA TRACE COC 6A00 FTn V3.0-P261 OPT=0 03/10/72 12. 53,4** PAGE l 

SUBROUTINE CALIBRA 

COMMON TIME UOU *VetOC (6*100) ,NCHANNC*LENARR*NFILE.NRECOR*IEXIT 
COMMON /BCAIBR/ NCALVAL* ACT VnLT (5) .SLOPE <61 « ZEROTAP (61 « 

. STANDEVI6*51VvARI 

DIMENSION SUMCAL <61 .SUMTAP {6) .SQ VAlUE 1 61 ♦ SUMAC T (6) * ACT X T AP (6) « 

. 5UMSQ (6.51 .RECMEAN<4 1 «TOTHE AN 16.51 ♦ TEPMEAN <6 ) . SUME AN <61 i 

, TEHPSUM (6) 

ICHt CK = 0 
NSAMPLE - 0 
LASTCAl - 0 
ICALVAL = l 
DO 100 1-1. NCHANNC 
SUMEANtl) = 0.0 
TEMPSUHtll = 0,0 
SU M C AL < I > = 0.0 
SUMACTU) = 0.0 
SUMTAPm * 0.0 
SQVALUEUI = 0.0 
ACT X TAP ( I ) = 0.0 
TEpmean < 1 1 = 0.0 
RECMEANI I) = 0.0 
00 100 J*l* NCALVAL 

totmeani r* j> = o.o 

100 SUMSOU.J) = 0.0 
105 CALL flUFANE 

GOTG0UF= 3H NO 

IF I ICALVAL .EO. NCALVAL) LASTCAL » LASTCAL ♦ 1 
oo no k=i.lenarr 
00 110 1 = 1 * NCHANNC 

no sumcal m= sumcal<u ♦ veloc<i*ki 
if < i check .gt. o) go to i3i 

NSAMPLE = NSAMPLE ♦ 1 
DO 125 1=1. NCHANNC 
RECMEAN<I>- SUMCAL < I l/LENAWR 
IF INRECOW .EO. 1) GO TO 120 

IF ( &ECME AN < I > « Gt . TOTMEANI I. ICALVAL) ♦ VARl .OR. RECM£AN< 1 1 .LT. 

. TOTMEANI I* ICALVAL! - VARI) GO TO 130 ' 

120 IF (I *EQ. M ^MlTEIfafllNRECOP. ICALVAL. ACTVOLTUCALVALI 
1 FORMAT l|H0»5X*RtC0WD ME ANS*4X*REC0 hD NUMBER* 14 * 7A*CAL1HRAT ION* 12 *4 
, X* I NPUT VAluL*F5.1/1 1X*CHANNEl * 1 OX»MEAN* 13X*CUMUL AT 1 VE MEAN*6K*NUM 
.BEW RECORDS FOR CUMULATIVE MEAN*! 

00 123 K=1.LENARR 

123 SUMSOH. ICALVAL! = SUMSO < 1 * I C aLV AL I ♦ VELOC < I »K) **2 
SOME AN ( 1 1 = SUMEANII) ♦ RECmEaNII) 

SUMCAL < 1 ) = 0.0 

TOTMFANU* ICALVAL) = 5UMEAN { I \ /NSAMPLE 
125 WRITE 1 6.2) I .WFCmEAN U 1. TOTmEAm < I . ICAL VAL ) .NSAMPLE 

2 FORMAT < 1H .12*. 12.10 X.Fd. 4. l<*t.Fd. 4*25X13) 

GO TO 105 

110 NSAMPLE= NSAMPLE - 1 
ICALVAL = ICALVAL ♦ 1 

131 IF < ICALVAL .GT. NCAlVAl .ANU. LASTCAL .GT. 3) GO TO 160 
ICMfCK * 1CHECK *■ 1 

MR I TE < b. 3) NWECOR * ICALVAL. ACT VOLT < ICALVAL I 

3 FORMAT (1H0.5X*TEMPORAWY ME AN $*8X •RECORD NUMBER* J4. 1 OX •CAL I BRAT ION* 



SUBROUTINE CAL I BRA TRACE 


CPC 6400 riN V3.0-P261 OPT«0 02/10/72 12.53*49 


PAGE 
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6ft 


69 


7ft 


7S 



105 


tin 


. J2.4X* INPUT VALUE*FS.1/!1X*CHaNNEL* 10X*M£AN*I 
00 140 1 = 1* NCrtANNC 
PECHEANU) = SUMCALUl/LENARR 
WRITE<6«4) I.RECHEAN(I> 

4 FORMAT < 1 H , 12X12. 10XF8.4* 

suMCAtm = o.o 

IF (1CHECK .£Q, 1) GO TO 135 
DO ! 37 K=1.LENARR 

137 SUwS0U*ICM.VAU=SUMSQ<IfICALVAL> ♦ VELOC 1 1 » K1 **2 
TEMPSUM ( 1 I = TEMPSUM ( 1 1 * RECmEAN<I> 

135 If < REC M f AN ( [ 1 .GT. TOTMEAN 1 1 * ICAL VAL-11 ♦ VARI .OR. RECMEANUO 
* .LT. T0TM£AN(I,ICALVAL-1> - VARI) G0TQBUF=3NYES 

140 CONTINUE 

IF IICHECK.GT. 3) GO TO 160 
IF (GOTOBUF .EQ. JHYESI GO 10 105 
DO 150 1=1. NCHANNC 
TEMPSUM! I ) = 0.0 
150 sumsou.icalvad-o.o 
ICHECK = 0 

ICALVAL = ICALVAL - l 
GO TO 105 

160 I END * ICALVAL - l 

WMlTE!6.S> IEND.ACTVOLT(IENO) 

5 FOR M AT (1H0*5x*ST ANDARD OE V I AT IOnS* 1 OX*CAL IBRAT ION* I2*5X*|NPUT VALU 
.t»f 5. 1/1 1X*CHANNEL*10X*RMS«) 

00 170 1=1. NCHANNC 

STAnOEV i I . 1 C AL V AL-1 ) = SORT <SUm$Q I I • ICALVAlM ) / INSAMPLE«LENARR 1 - 
. TOTMEAN 1 1 . ICALVAL-l I **2> 

170 WRITER. 6) I.STANDEV(I.IEND) 

6 F ORMAT ( 1 H , 12A12.7XF9.3) 

ns ample = i check - l 

00 ITS 1=1. NCHANNC 
SUM E AN | I ) s TEMPSUM { I 1 

TOTMEANU. ICALVAL) = TEMPSUM < \ ) /NSAHPLE 
ITS TEMPSUM It) ^ 0.0 
ICHECK = 0 

If MCALVAl .LE. NCALVAL) GO TO 10* 

180 MR 1 TE <6* 7 > NWECOR 

7 format < ) ho * 5 x * actual vs tape voltage*iox*le*st square hetkoo*sx*nu 

.MBER RECORDS USEO FOR CALCUL AT IONS* 131 
DO 200 1=1. NCHANNC 
DO 190 J=l. NCALVAL 
SUMTAp(I)=SUMTAPm ♦ totmean<i*ji 
SO VALUE 1 1 > =SOVALUE 1 1 1 ♦TOTMEAN 1 1 • J> **2 
ACT X TAPII1 * ACT X TAP<I> ♦ TOTMEAN < I • JI*ACT VOLT tJ > 

190 SUMACT(I) = SUMACT(I) ♦ ACTVOlTU) 

SLOPfc(I) = <SUMACTU)*SUMTAP<|) -NCALVAL*ACT X TAP< 1 1 1/ 

. ISUMTAPtI >*«2-NCALVAL* SUVALUEim 

/ERn T AP < I ) = <SUMTAP! 1 1 • ACT X TaPU) - SUMACT ( I) *SOV ALUE ( 1) 1/ 

. { SUMT AP < | ) **2“ NCALVAL*SQVALUE(I>) 

*R I TE < 6. B> I. (ACIVOLT(J) ♦ TOTMF. AN I 1 * J ) *J=l *NCALVAL) 

8 FORMAT UHQ.10X*C H A N N L L* 13/1 5X*V ALUES USEO FOR LEAST SQUARE 
.CALCULATIONS*! OX* INPUT VALUE*5X*UPE VALUE*/ <69XF4*i*llXF6*3)l 

200 *RITE<6*9) SLOPt ID.^EROTAPd) 

9 FOHHATUH . 1 5X* VALUES OBTAINED FROM LEAST SQUARE CALCUL AT ! ONS*7X*S 


t 



SUBROUTINE CAUBRA TRACE CDC 6400 FIN V3.0-PZ6I OPT»0 02/10/72 12.S3.49. 

.L0Pe*8K»!NTCRC£PT*/68XF5*3*lHF5.3) 

RETURN 

ENO 
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% 



SUBROUTINE BUFaNE TRACE CDC 6400 FTN V3.0-P261 OPT*0 02/10/72 12.53.4R, 

SUBROUTINE 8UFANE 

COMMON T !ME aOH * VCLOC I 6* 100 > »NCHAnNC*LENARR*NF I LE ♦ NRECOR* I EX IT 
COMMON/BBUF ANE/ IOENT «IPaRITV*lPACDAT ♦EOFMUL #NTOTF IL »PLOT 
common/unpk/ i time » i C omoro < 200 > * idatvrd< ioooi 

S 100 00 105 I=1«LPAC0AT 

ios icowwpom = o 

LUNPDAT = LPACOAT * 5 

oo no i - i*lunpoat 

J10 IDATWRO(l) = 0 

10 115 NRECOR * NRECOR ♦ l 

120 BUFFER IN U * l ) ( i T 1ME ♦ iCOMwRO (LPACDATM 
125 IF lUNITUM 140*130*135 
130 NRECOR = NRECOR - 1 

WRITE (6*1) NRECOR* NFIlE 
1^ NRECOR = 0 

NFILE = Nf ILE *1 

IF (NFILE .GT* nTOTFRI GO TO 136 
IF (IOENT *EG • 3HYES)CALL HEAqER 
GO TO 100 

2A 135 IPARITY » J P AP I T Y ♦ 1 

WRITE (6*2) NRECOR* NFILE 
NRECOR = NRECOR - 1 
WRITE (6*3) (PARITY 

GO TO ns 

2* 136 JF (PLOT «FQ • 3HYES) IEXIT = 3HYES 

IF (PLOT *EO. 3HYES) GO TO 150 
CALL EXIT 
140 CALL UNPAK 
CALL SOPTANE 

30 1 FORMAT ( 1 MO * * THERE ARE *14* RECORDS ON FILE NUHBERM31 

? FORMAT UNO.* PARITY ERROR OCCURRED ON RECORD NUHBER*|4* FILE NUMB 
. ER # 1 3 ) 

3 FORMAT ( 1 HO * * THERE HAVE BEEN»13* PARITY ERRORS* 1 
150 RETURN 

3^ END 



Q> 
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SUBROUTINE 






HEAOEft TRACE CDC 6400 FTN V3.0-P261 0PT=0 02/10/72 12*53.49. PAGE 1 

SUBROUTINE HEADER 

COMMON T!HEUOl>*VELOCl6tlOO>.NCHANNC»LENARR»NFlL£fNR£COR*iexlT 
COMMON/BH£ADER/IO(2) 

BUFFER IN < 1 • 0 ) U0m*I0(2) > 

IF (UNIT ID) lG0*llO»UO 
110 WRITE <6«2> NF1LE 

2 FORMAT (1H0.* PARITT ERROR OR EOF OCCURRED IN HEAOER OF FILE NO* 

♦ n> 

100 WRITE!*.)) IDfNFILE 

1 FORMAT ( iHl « * HEADER IN BINARY *2AiQ* ON FILE NUMBER *12) 

120 RETURN 
END 





SUBROUTINE SORTANE TRACE 


CDC 6400 FTN V3.0-P261 OPT*Q 02/10/72 12.53.49 


S 


14 


IS 





SUBROUTINE 50RTANE 

COMMON TIME non • VELOC I6t 1001 .NCHANNC»LENARR»NFILE*NftECoft'IEXIT 
COMMON/PSOPT/lBEGSKP.WRJTI)AT*!SKlPfFACTOR 
COMHON/UNPK/1TIME«ICOMWROI200I tlOATWROUOOO) 

M= I BE6SKP 

00 100 Isl.LENAPR 

VELOCUtll = IDATWRO(M) * FACTOR 

VEL0CI2*I> - I0AT*R0<M*1> • FACTOR 

VEL0CI3*n = JDAT*RD|M*2 )*FACtOR 

VELOC 1 A * I ) = 1QATWRD|M*3> * FACTOR 

vfLoc(Sti) - idat*pdim*4) * factor 
VELOC 16 ♦ 1 1 = IDaT#RD<M*5) * FaCTOR 
100 M = M ♦ ISK1P 

IF ( VRITOAT .EO. 3HTES> CA|_L 1>ATA*RI 

RETURN 

ENO 



subroutine oatamri trace 


COC 6400 FTN V3.0-P26I OPT-O 02/10/72 12,53,4R 


PAGE 


1 




lft 


SUBROUTINE OATAMR! 

COMMON TIME 1 101 1 « VELoC (6* 1001 *NCHANNC*LENARB*NF ILE*NRECoR* iexit 
COMMON/UNP*/ [TIME •ICOMWRD <2001 *I0ATWRO|1000) 

00 10 I=I*NCHANNC 

10 WRITE <6*11 1 « (VELOC (!«J)*J*1*LENAMR> 

1 rORMATUH0»10A,* ANEMOMETER VELOCITY DATA* LEVEL NUMBER* 12/1 01 

. ItOFll.S/n 

PRINT 2 » I T IME 

2 FORMAT <1H0«10X»MTIME AT BEGINNING OF RECORDMIOI 
RETURN 

END 



SUBROUTINE SKIpEOF trace 


COC 6400 FTN V3.0-P261 0PT*0 02/10/72 12.S3.49. • 


PAGE 


l 


5 


1* 


15 


?r> 


25 


30 


35 



50 


55 


SUBROUTINE SKlPEOF 

COMMON TIME (lOU *VELOCl6*lOOJ *nchannc.lenarr ?nf ile. nrecor* I ex it 
COMMON/BBUFANE/IOENT. lPAPITr.LPACOAT.EOFMUL.NTOTFlL*PLOT 
COMmON/UNPK/ I TINE* 1C0MWRD (200 > ♦ IDAT WRD (1000) 

COMMON/dSK I peo/nf ilskp .nrecskp 
NREC = 0 

IF INFILE *GT. NFILSkP* GO TO 125 
100 BUFFER 1N<1.]HITImE,IC0M¥RD<lPACDaT>> 

IF <UNlTnnU5.120iU0 
110 NREC = NREC ♦ 1 
LEN = LENGTH ! M 
NRECOR = NPECOR * I 
PRINT 2tNFlLEiNPEC0R,NRFC*LEN 

2 FORMAT UH0.5X*PARITV ERROR OCCURRED WHILE SKIPPING FILEM2* RECORO 
,* I4/7X 1 3* RECORDS HAVE Ht£N SKIPPED. LENGTH WASM4) 

GO TO 100 

115 LEN = LENGTH! 1 > 

NREC = NREC ♦ 1 
NRECOR = NRECOR ♦ 1 

IF (LEN *NE. LPACOAT ♦ 1) WRUEC6.J) LEN .NRECGR.NF ILE .NREC 

3 FORMAT (1H0, 5X*REC0RD SKIPPED OF IMPROPER LENGTH. LENGTH WASM4/7X 
♦•RECORD*!** F ILE* 12*2X12* RtCOROS SKIPP£0*l 

GO TO 100 

120 WRITE 16.41 NREC * NF ILE .NRECOR 

4 FORMAT ( 1 HO ♦ 5X*Sk IP PE D - * 13* RECORDS ON F ILE*I2* THERE WERE*I4* RECOR 
.DS ON THIS FILE*! 

NFILE = NF1LE * 1 

IF IIPENT .EQ. 3HYES1 CALL HEADER 
NREC = 0 
NRECOR = 0 

IF INFILE .LE. NFILSKPIGD TO 100 
125 IF (NRECSKP .EQ. 01 RETURN 
DO 160 1 - 1 • NPECSKP 

BUFFER IN ! 1 « I > I IT IME » ICOMWRD (LPAcOAT ) I 
IF < UN IT (1 ) ) 130*160.140 
130 NRECOR - NRECOR ♦ 1 
NREC = NREC ♦ 1 
LEN = LtNGTH(l) 

IF ILCN .NE. LPACOAT * \i WRirE(6ii> LENtNRECOHtNF ILE.NPEC 

GO TO 160 

140 NREC = NREC ♦ 1 

NRECOR = NRECOR • 1 
LEN = LENGTH (11 

WRITE 16.51 NRECSKP.NF ILE .NRECOR *LEN* NREC 
5 FORMAT UH0*5A*PAR ITT ERROR OCCURRED WHILE SKIPP1NGM2* RECORDS ON 
.F1LFM2/7A* RE CORO *14* LCNGTHM4* RECORDS SKIPPEDU3) 

GO TO 160 

150 wR I TE (6.6) NF ILE .NRECOR .NREC 

6 FOWMATUH0.5X*EOF OCCURRED WHILE SKIPPING RECORDS ON F ILE* I2/7X*LA 
.ST HECORO* 14*2X13* RECORDS HAvE BEEN SKIPPED*) 

160 CONTINUE 

WRITE<6»7) NPEC. NFILE 

7 FORMAT (1H0.5X.I3* RECORDS HAVE BEEN SKIPPED ON FILE*I21 
RE Turn 

END 




SUBROUTINE VOLT AO J TRACE COC 6400 FTN V3.0-P261 OPT*0 02/10/72 12.S3.49. 

SUBROUTINE VOLT AO J 

COMMON TIME 4 1 01 1 . VELOC 16. 100 > .NCHANNC.LENARR.NF ILE.NRECOR. I£xlT 
C0MM0N/UV0LTAD/ICHAN6E. SCALE 16) .ZEKOACT 16) .V0LTCHG.TIMLCH6. 

. ISCALE (61 

S IF ( ! CHANGE .OT. 0) READ (5. I H ISC ALE ( I) * 1 = 1 .6) . VOLICHG. T IMECHG 



20 

DO 90 1*6 

GO TO 130 *40*50* 60 *70*601 * 


30 

GO TO 131*32*33)* ISCALE ( I ) 


31 

SCALEU) - 40.166 

1ft 


ZEROACTU) - 2.799 


32 

GO TO 90 

5CAL E ( I ) = 70.067 

IS 

33 

zEPOAcnn = 2,4i3 
GO TO 90 
SCALE!!) * 0.0 


40 

zeroactui - o.o 

GO TO 90 

GO TO (41,42*43) *I5CALE(1) 


41 

SCALE < 1 > = 42*161 

?« 


ZERO ACT ( [ ) s 2.163 


42 

GO TO 90 

SCAlEtn = 01.437 

2S 

43 

ZEROACTU) = 2.261 
GO TO 90 
SCALE!!) « 0.0 


50 

ZEROACTU) * 0.0 
GO TO 90 

GO TO (51 *52*531 * ISCALE (I) 


51 

SCALE!!) = -2.961 

3ft 


ZEROACTU) = 2,057 


52 

GO TO 90 

SCALE! 11 « 03.606 

3S 

53 

ZEROACT ! 1 ) - 1.083 
GO TO 90 
SCALEU) » 0.0 



ZERO AC Ml) = 0,0 
GO TO 90 


60 

GO TO (61 *62*63) ISCALE*!) 


61 

SCALEU) = 4*2.069 

4ft 


ZEROACT i !.) * 3.674 


62 

GO TO 90 

SCALE! 1) a 83.224 

4S 

63 

ZEROACT(f) =3.065 
GO TO 90 
SCALEU) * 0.0 


70 

ZEROACTU) = 0.0 
GO TO 90 

GO TO (71*72*73) » ISCALE ( 1 1 


71 

SCALE U ) * 47.070 

SO 


ZEROACTU) = 0.075 


72 

GO TO 90 

SCALE!!) = 93.300 

ss 

73 

ZEROACTU) a 0.330 
GO TO 90 
scale < \ ) a 0.0 


PAGE 1 



SUBROUTINE VOLT AD J TRACE 


CDC 6400 fTN V3.0-P261 OPT*0 


02/10/72 12.53.49* 


PAGE 


60 


6*i 


70 


7S 


ZEROACTH) * 0.0 
GO TO 90 

60 GO TO (81*62*83) ♦ I SCALE ( 1 1 

81 SCALEUI = 40.217 
ZEROACTUI * 3.764 
GO TO 90 

82 SCALEUI “ 77.313 
ZEROACTUI » 3.639 
GO TO 90 

63 SCALE U I = 0.0 
ZEROACTUI = 0.0 
90 CONTINUE 

ICHANGE = ICHANGE ♦ 1 

1 FORMAT 16(12)* A3 *F5. 31 
*RIT£(0.2I 

2 FORMAT ( l HQ .5* ° ACTUAL VOLTAGE V$ WElOCI T V*5X»REGRESSI0N VALUES*/10< 
.*L£vEL*5X*SLQPE*SAMNTERC£PT*) 

DO 100 1 = 1* NCHANNC 
100 *«!TE<&.3> IfSCALEUltZEROACTU) 

3 FORMATUM * 11 X II * 7 AF6 .3*6XFS. T) 

RETURN 

tNO 



subroutine; plotvel trace 


CDC 6400 FTN V3.0-P261 0PT*0 02/10/72 12*53*49 


PAGE 


I 


s 


1ft 


is 


?0 


SUBROUTINE PLOTVEL 

COMMON TIME <1011 * VELOC <6* ) 00) .NCHANNC»LENAHR *NF ILEtNRECOR. IEXI T 
COMMON/BPLOTVE/ELEV (6) *SUMAVE f 6 i t IsAMPLE iLAflELA (4) 'LABEL V (4) . 

. LTITLE14) .MULTlHEilDATAPT 

DIMENSION AVEVELIM 
MULTIME = MULTlMt ♦ 1 
*«ITE tb'M TIME(IDATAPT-1)'N«EC0R 

1 FORMATI1mQ,SX«YELOCITY PROFILE PLOITEO AT T INE«F9* 3*5A*RECO«D NUNfl 
*EP*I4/10X*VALUES USED FOR PL0T*5X*LEVEL*5X*ELEYATI0N* H,*5X*VEL0CI 
.TY, M/SEC*) 

00 100 I = I.NCHANNC 
100 AVEVELUI - 0*0 

DO MO I = 1'NCHANNC 
AVEVELUI = SUMAVE ( 1 )/ 1SAMPLE 
MHITE<6'2> I'ELEVUMAVEVELM) 

2 FORMAT MN ♦ 3bX [ 1 * 1 QXF6. 3. 1 3XF* , J* 

110 SUM A YE 1 I } r 0.0 

CALL IDIOT ( AVE VEL' ElE Y *6' 2*Ot)M*LAbELX*LA0ELY *LT I TLE« ♦ 1 1 

isample = n 

IF f It A I T .FO. 3HYES) CALL EXIT 

RETURN 

FNO 



UMPAK 


COMPASS - VER 2.0 M 02/10/72 12»S4.I6. 


IDEnT UNPAK 

INSERT LENGTHS OF PACKED AND UNPACKEO ARRAYS 


10 


15 


20 


25 


30 


35 


AO 


IS 



311 


lengtha 

SET 

201 




1 750 


LENGTHS 

£ 

SET 

1000 







USE 

/UNPK/ 



0 



NE 

BSS 

LENGTHA 



311 



B 

faSS 

LEnGTHB 







USE 








ENTRY UNPAK 



0 



UNPAK 

BSS 

1 



l 

717000000] 



SX7 

ia 




7100004000 



sxo 

40006 



2 

43214 



MX2 

lZ 




5110000000 C 



SA1 

nE 

Ad) 

FIRST WORD OF ARRA 

3 

6160000310 C 



586 

fi-1 

B4JI 

BASE 


6170002260 

c 


Sd7 

r*lengthb-i 

B (LAST I 

4 

6150000060 



SBS 

48 




61 10000074 



SB1 

60 



s 

5011000001 


GE 160 

SAi 

Al*l 

GET 

All) 

6 

6166000001 


GET 12 

SB6 

R64| 




11621 



BXb 

*2*X1 

mask 

OUT 12 BITS 


67515 



SBS 

P 1 -B5 


RIGHT SHIFT 

7 

22656 



LA6 

RS*X6 


BUT 


67515 



SBS 

Rl-d5 


AVOID SIGN EXTENSION 


11760 



BX7 

*6*X0 


CK FOR SIGN BIT 

10 

0307000011 ♦ 



ZR 

< 7 * STqRB 


MASK OUT SIGN BIT 

15660 



BXb 

-X0«X6 



14666 



BX6 

- A6 



11 



STORB 

BSS 

0 



11 

21601 



AX6 

\ 

DELETE ZERO-BIT RIGHT-FILL 


56660 



SA6 

P6 

STORE IN BIJ) 


0467000000 

♦ 


€Q 

P6 *87 « DONE ♦ 


12 

0550000014 ♦ 



NE 

rS»B0« INH10 



6150000060 



SB5 

48 



13 

43214 



HX2 

U 




0400000005 ♦ 



EQ 

T-ET60 



14 

20260 


INMID 

LX2 

4d 




6155777763 



SBS 

HS-12 



15 

0400000006 * 



EQ 

GET 12 




0 

♦ 

DONE 

EGU 

UNPAK 




16 


46302 


END 

STORAGE USED 
6400 ASSEMBLY 


44 STATEMENTS 
0*341 SECONDS 


10 SYMBOLS 
23 REFERENCES 
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APPENDIX A-3 


Computer Program for Determination of Temperature 
and Humidity Profiles 



PROGRAM 


TEMPHUM TRACE 


PAGE 




1 A 




?n 


2* 


3" 


3S 


4n 


4* 


Sn 


S«5 - 


CDC 6400 FTN V3.0-P261 OPT*0 


02/10/7? 08. 46. 42. 


PROGRAM TEMPHUM ( INPUT tOUTPUT p T APE5= INPUT » T AP£6=0UTPUT ,T APE 1 # 

. filmpu 

CS DEBUG 

Cf ARRAYS 

COMMON NF IlE »NR£ COR, LENA RR.NChANN, TEMP < )(>* 1001 « WRITDAT 
C0MM0N/8S0RT/ I0EGSKP .FACTOR *JSK IP 

COMMON /0 INSTCA/NINSC AL • VARl IN.EXC1TVO.RESISC1 0,2) tCALRLS 1 10*21 v 

• GAIN ( 10) 

COHMON/BCAL I BR/ SLOPE f 10 1 .ZEROTAPC 1 0 1 

COMMON/ B0UFT EM/ I DENT ,HULE0F.1FXIT»NT0TF IL,PL0T » 3 PAR I T Yf ZERO! IH, 

• JCLOCK 

COMMON/UNPK/ITIME, ICOMmpD( 200I • IOATMRO (1000) 
COMMON/HSKPEOF/aAOATA,LPACUAT,NRECSKP.NFlLSKP 
COMMON/BHUM I l)/S I GM AtBARPPES « MF ATLAT • CP.CPV.HUM JOI CIO) 
C0HM0N/8TAPECA/NCALVALtVARnP«ACTV0LH5> 

COMMON/BPLOTEM/AVETEMPC 10) , I 1m£ »EL£V (6 ) ,LA0£LX (4) ,LABELY (4) .LTITLE 

• (4) 

DIMENSION RESISR(IQ) «SUmTEMP<10> 

DATA LA0ELX/4OH TEMPERATURE ,C / 

. .LABELY/40H FLEVAUON. M / 

• tLT ITLE/**0H TEMPfRATUkE PROFILE / 

read is. 1) C4L TAPE. CAL INST. *RI TO AT. IDENT.mJLEQF, PLOT. NCHANNiLENARR. 

• NAVERFC, IBEGSKP* 1 5K [P .N INSCAL .LPACOAT .NTOTF IL.EXCIT VO, 

• PESISRU) . A.B'C.O.E. T IMRAT.VAR! IN .SIGMA. BARPRES.CP. 

. HEATLAT.CPV.VARITP. <ACTVoLTU).I = lt5).CRESISU.n# 1*1,10 

• ) . « RESISC1 ,2),] = 1, in) .(CALRESCI ,11,1 = 1 ,10) • CCALRES (1,2 1 , 

. 1 = 1.10) , (ElEVCU .1 = 1*6 ).NCALVAL 

I FORMAT l6A3,8!4,3F&.3/4FR.3*2FW,5F7.3/iFS.2,5<F5.2),lQ(F5,2)/ 

. 10(ES»2) *6lF5.2)/4CF5.2),l()<F5.2)/6|F6.3l . 13) 

REMIND l 
PRINT 3 

3 FOPMAT<1HO.SX*NOTE.*. CHANNEL I IS LEVEL 2. AMBIENT TEMPERATURE*/ 12 
.X*CHANnEL 2 IS LEVEL 3. ORY*/l2X*ChANN£L 3 IS LEVEL 1. 0RY*/12X*CM 
* ANNEL 4 IS LEVEL 1, WE T*/ i£X*CHANNtL S IS LEVEL 4, DRY*/ 12X*CHANNE 
.L 6 IS LEVEL 4* MET*/ 12 X*CH ANnEL 7 IS LEVEL S. ORY*/ 12X*CHANNEL 8 
.IS LEVEL 5. *ET*/12X*CHANNEL R IS LEVEL 6. ORY*/ 12X*CHAnnEL 10 IS 
.LEVEL 6, WfT*1 

PRINT 4, CALTAPE,CAlINST.VRItOAT,IDENT,MULEOF,PLOT.NCHANN,LENARR. 

. N AVEREC, IREGSkP. ISK IP .N InSCAL .LPACOAT iNTOTF1L*EXC1TVO, 

. RESISRU ) t A .b. C, D. E, T IMR AT. VARIIN. SIGMA, BARPRES.CP, 

. HEAtLAT .CPV*VAWITp.(ACTVl>LT (I). 1 = 1, 5) , (RES IS 1 1 . 1), 1 = 1, 10 

. ) , (RES IS (1 .2) ,1 = 1 • in) i (CALRESUil > » 1*1. 10 1 • (CAlRESC 1.2) , 

. 1 = 1,10), CElEV (!)♦,! =1 ,6) iNCALVAL 

4 FORMAT I IHO , * CALTAPE =*A4* CAt INST =*A4* WRITOAT =#A4* JDENT = *A4 
.* MULEOF =*A<+* PLOT =*A4/* NC^ANN =*I3* LENARR =*I4* NAvERLC =*I5 
.* JPEGSKP =* 12* ISKIP =* 13* NlNSCAL =*I2* LPACOAT = *U* NTOTFlL =* 
.J2/* EXCITVO =*FS. 3* RES I SR ( 1 ) =*fW* A =*F6.2* B =*F7.3* C =*F7 
..3* 0=*F5.3* f = *F6.3/* T 1MRA T =*F4,1* VARUN =*FS.2* SIGMA =*F?.5 
.* BAHPRES =*FB.2* CP **FS.J* hEATLAT -*F6«1* CPV **F5.3* VARITP =• 
.F^.2/*ACTV0LTU THRU S>**5FS.|/* RLSISd THRU 10,11 =»10F5.2/* RES 
.IStJ THRU 10,2) =• 1 0F5. 2/* CAi.RESU THRU 10,11 =*10F5*2/* CALRES U 
«1 THRU 10,?) -*10F5.2/* ELEVCJ THRU t I =*6F7.3* NCALVAL =*121 

I EX I T = 3h NO 
JCLOCK = 0 
BAUATA = 3H NO 



PROGRAM 


TtHpHUM TRACE 


COC 6400 FfN V3.0-P261 OPT*0 02/10/72 08.46.4 2 


PAGE 


fen 


6* 


7tx 


1* 


Hn 




90 


9S 


mo 


10« 


im 


FACTOR = S0PT(2.O) / <2.0*«*9-l .0) 

NFILE # 1 
NRECOR = 0 
2ER0TIM » 0.0 
IPARITT * 0 
MULRLC * 1 
DO 100 1=1«NCHANN 

100 SUHTEMPin ■ 0*0 

IF ICALTAPE .EG. 3HYES) CALL TaPECAL 

if (cal inst .eq. 3hyesi call instcal 

NFILSKP = l 
NRFCSKP a 0 

IF (NFILE .le. 1) call skipeof 
JCLOCK = 0 
2EP0TIM - 0.0 
NRECOR * 0 
RFA|)(5,5)FPStftEC 
S F0PHAT(A3I 

IF (FRSTREC .NE. 3HBADI GO TO 200 
NFUSKP = 0 
NkFCSKP = I 

CALL skipeof 
200 CALL ttUFTEPP 

IF ( IEX I T .EQ. 3HYES1 GO TO 400 
00 300 1=1 .NCMANN 
00 300 K=1*LENARR 

300 SUMTEHP(I) = SUMTEHP ( f 1 ♦ TEMPlItKl 

IF (NRECOR .LE. NAVEREC*NULREC> GO TO 200 

GO TO 500 

400 NAVEREC - NRECOR- (NAVEREC*IMUlREC-U ) 

500 DO 700 1= l »NCH ANN 

AVf TEMP ( 1 1 = 5URTEMP ( I ) / (NAV£rCC*LENARRI 
SUHTEMP(I) = 0.0 

avetempui = (Slope ( 1 1 *avetehp ( 1 1 ♦ zepotapum/gainuj 

FACTORl = AVETEMPUI/EXCITVO 

FACT0R2= RFS 1 S < I » 1 )/ IRES I SR I I ) ♦ RESISHilH 

AVETEMPd) = RE5IS(1«2I/<FACT OR 2*f ACTOR 1 1 - RESI5U.2) 

IF (I .GT. 1 ) GO TO 700 
DO 600 J=2 *NCHANN 
600 RESISR(J) = A VETEMP ( 1 ) 

700 AVETEMP ( I i = A *9* AVt'TEMP ( l) ♦ C« AVETEMP ( I ) «»2 ♦D*AVETEMP(1 1**3 ♦ 

. E* AVETEMP f I ) **4 

TIME = (TlPRATMITIME-ZEROYlMj/lOOOO.I 
4R I TE < 6 * 2) T IMF » NRECOR* (AVETEMP ( II « 1 = 1 vNCHANN > 

2 FOP«AT(lM0.» TEMPERATURES AvEpAGED OVER 10 MINUTE INTERVALS. TINE 
* *F9.3» SECS. RECORD NUMBER* I 4//2*«CHANNEL i *4 X*CHANnLL 2*4X*CHAN 
.NEL 3*4X *Cm ANNEL 4*4X»CHANNEL S*4X*CHANNEL 6*4X*CHANNEL 7*4X*CHANN 
.EL fc»*4JM>CHANNEL 9*4X*CHANNEL lO*/3x*LCVEL 2*6X*LEVEL 3 # 6X*LEVEL 1* 
•6X*I EVEL 1*6X«LEVEL 4*6X*LLvEt 4*6X*LEVEL 5*6X*LEVEL 5*6X*LEVEL 6* 
♦6X«LE VFL 6*/3X* AMUIENT # BX*0RY*1 0X*ORY* 1 0X*tfET* 10X*DRT*1 0X*WET # lOX* 
# DRY*10X*WET*IOX*ORY*1CX*WET*/3XF6. J* C*5XF6.3* C*5XFt>.3* C*5XF6.3* * 

. C*5xF6 • 3* C*5 xF6. 3 * C*SXfb.J* C*5 aF 6.3* C*5XFb.3* C«* 5 XF 6 . 3 # C*> 

CALL HUMID 

IF (PLOT .EG. 3HYES) CALL PLOtEMP 
MULRLC = MULREC ♦ l 
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IF ( 1EXIT .NE. 3HTESI GO TO 2ft0 
END 
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SUBROUTINE TAPECAL 

COMMON/BTAPECA/NCALVAL* VARITP* ACTVOLT I 51 
COMHON/BCAL IBR/SLOPE UO > ♦ZEROTAP ( Ko> 

COMMON NF ILE «NRECOR*LENARRtNCHANN« TEMP 1 10. 100) » VRITDAT 

dimension sumcal no) *$uhtap( in) »sq value t io) » Sum act ( 10 ) §act x tap (I 

• 01 * SUMSQ I 10 tS) tRECMpAN I 10 ) »TOTM£AN I 10«5) » TEPHEAN (10k *SU 

* MEAN! JO) .TEMP SUN 110 1 iSTAnDEV 110*5) 

1CHECK “ 0 

nsample = o 
lastcal = 0 

ICALVAL = ! 

DO 100 I=a*NCHANN 
SUMEANCII = 0.0 
TEMPSUM ( I ) s 0.0 
SUMCAL(I) = 0.0 
SUM ACT ( I ) a 0.0 
SUM 7 AP 1 1 ) a 0.0 
SQV AtUE ( 1 1 s 0.0 
ACT X TAPI!) - 0.0 
TEPMfANII) = o.o 
RECMEANII) * 0.0 
DO 1 GO J= | tNCAL V AL 
STAnOE^(J#j> = 0.0 
TOT m£an ( 1 • J I a 0.0 
100 SUMSQ ( 1 * J ) * 0.0 
200 CALL BUFTEMP 

GOTOBUF = 3H NO 

IF UCALVAL .eo. ncalvau lastcal - lastcal ♦ 1 

DO 300 K=1 .LENARR 
DO 300 I=1«NCHANN 

300 SUMCAL I I ) * SUMCAL 1 1 ) * TEMPI I »K) 

IF IICHECK .GT. 0 ) Co To BOO 

nsample = nsample ♦ I 

MR I T£ <6* I I NRECORg IC ALVAL « ACT VOLT I1CALVAL) 

FORMAT I 1H0 *5X*RECORD ME ANS**X*R£CORO NUMBER* 14 *7X*CAL I8RATI0NM2* 
.AX - INPUT VALUE*F5.)/1 lX*CHANNFL* 1 0X*N£AN* I 3X*CUMUL AT I VE MEAN* 6X 
.•NUMBER RECORDS FOR cumulative MEAN 4 ) 

DO 600 I-I *NCHANN 

RECMEANII) = SUMCAL I I) / LENARR 

WRirEI6*2U.KECMEANm 

2 FORMAT ( 1H .12X.!2.10X*F«.4) 

IF INRECOR ,LQ. 1) GO TO 400 

IF (RECMEANII) .GT . TOTMEANII.IC At V AL ) ♦ VARITP .0. RECMEANII) 

• .LT. TOTMEANU tICALVAL) - VARITP) GO TO 700 
400 00 SOO * = 1 .LENARR V 

SOO SUMS 0 1 1 1 1CALV Al ) = SUMSOII* ICAU/AL) * TEMP 1 1 »K )**2 
SUDANI I) = SUM E AN ( I ) ♦ Rt CMLaN ( I ) 

SUMCAL 1 1 > = 0.0 

TOT ME AN 1 1 » ICALVAL ) = SUmE An ( 1 1 /NSAMPLE 
600 4RITE(6i4> TOTMEANU . ICALVAL) .NSAMPLE 

4 format UH»* 4GX.Fb.4*?sx*i3) 

GO TO ?00 

700 NSAMPLE a NSAMPLE - I 
ICALVAL = ICALVAL ♦ I 

800 ir UCALVAL .GT. NCALVAL .A. l*STCAL .GT. 3) GO TO 1300 


I 
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7* 


7* 


80 


«« 


9n 


9* 


100 


10* 


no 


I CHECK = ICMECK ♦ I 

WRI TE <6*51 NRECOR* ICAL VALt ACTVOLT < 1CALVAL) 

S FORMAT <lH0*5X*TEMP0RARr MEANS*6X*R£CQ«D NUMBER* 14* 1 0X*CAL IBRAT ION 
,*I2*4X*INPUT VALUE *fS.l/llX*CHANNEL*10X*MEA!|*! 

00 1100 1=1. NCHfiNN 

RECME AN < 1 1 * SUMCALUl/LENAPR 

*RITE<6*6)IfRECMEAN(I> 

6 FORMAT < 1H *12XI2*10XFB.4> 

SUMCAUU = 0-0 
'IF (ICHECK .EG. 11 GO To 1000 
00 900 K= l «LEN ARR 

900 SU«SO< niCALVALl = SUHSO ( I • ICaLVAL) ♦ T£HPUtK)**2 
TEMPSUMUI * TEMPSUM ( I ) ♦ RECuEANUI 
1000 IF CRECMEAMII .GT* TOTMEANII.ICALVAL-I) ♦ VARITP .0. RECME AN < 1 i 
. *L T • T0TmEANU*ICALVAL- 1> -VARlTPl GOTOBUF » 3HVES 

1100 CONTINUE 

IF ( ICrtECK .GT. J) GO TO UoO 
IF (GOTOBUF .EQ. 3HYES) GO TO 200 
DO 1200 1-1 *NCHANN 
TEMPSUMd) = 0.0 
1200 SUM5Q(I*ICALVAL> - 0.0 
ICHFCK = U 

ICALVAL = ICALVAL - l 
GO TO 200 

1300 IENO = ICALVAL - l 

HM1TE(6*BMFN0* ACT VOLT ( I END I 

8 FORmaT(1hO./*5x*STANOARO DEVI *T I0N5* 10X*CAL IBRAT ION* I2*5X* INPUT VA 

♦ LUE*FS.|/HX*CMANNEL*10X*RHS*) 

DO 1400 !M«NCHANN 

STANDEV<ItrCALVAL-l> * SORT (SuMSQ I I • lCALVAL-1 I / <NSAMPLE*LENARR) - 
. ToTMEAN < I * ICALVAL- 11 **2) 

1400 «VRlTE(6*91I*STAN0EV(i*lEN0> 

9 FORMAT ( 1 H *12XI2*7XF9,3> 

NSAMPLE * ICHECK - 1 

00 1500 I*1*NCHANN 
SUME AN < ! 1 - TEMPSUM ( 1 1 

TOTNEAN ( { * ICALVAL I * TEMPSUM < 1 1 /NSAMPLE 
1500 TEMPSUM < I > * 0.0 
ICHECK =0 

IF (ICALVAL .LE. NCALVAL) GO TO 200 
1550 *R! TE (6*10) NHECOR 

10 FORMAT () HO, 5X* ACTUAL VS TAPE VOLTAGE*! 0X*LE AST SQUARE M£TH0D*5X*NU 
*MBER RECORDS USED FOR CALCULAT IONS* 13) 

DO 1700 I -1 .NCHANN 
00 1600 J=1*NCALVAL 

SUM TAP (I ) = SUMTAPUI ♦ T0Tm£aN(1*J) 

SOVALUtm = SOVALUCin ♦ T0 Tm£AN(1*J1**2 
ACT X T AP I T ) * ACT X TAP < 1 1 ♦ TO TM£AN< I* J) * ACTVOLTIJl 
1600 SUM4CT ( 1 I * SUMACTUl ♦ ACTVOlT(J) 

SLOPEm -(SUMACTU) * SUMTAR(I) - NCALVAL*ACT X TAPUM/ 

. (SUMTAP 1 1 ) **2-Nr AL V AL*SOVALUE (I) I 

ZEROTAPU) = ( 5UMT API 1 ) *ACT X TAPUI - SUMACT 1 1) *$QVALUE< 1 1 1/ 

* (SUMTAP( n**2-NCALVAL*SQVALUE<m 

WRJ TE <6* 1 1 ) I * (ACT VOLT ( J) * TOTmEAN U » Jl * J=l «NCA lVAL> 

11 FORMAT UH *10X*C N A N N E L* 13/ 15x*VALUES USED FOR LEAST SQUARE C 


|o 


SU9R0UT1NE 
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.ALCULATIONS*10X*INPUT VALUE*5**T APE VALUE*/ «69XF4. I ♦ U XF6.3M 
1700 WRITEI6* 12) SLOPE < U • ZEROTAP U > 

12 FORMAT < 1H0* 15X* VALUES OBTAINED FROM LEAST SQUARE CALCULATIONS* 7X* 
.SL0PE*8X*lNTERCEPT*/6BXFS.3#llXF5.31 
RETURN 
END 
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SUBROUTINE INSTCAL 

COMMON/0 INSTCA/NINSCAL tVARI IN *EXC ! T VO .RES I S O 0 .2 > » CALRES U0*2> * 

. GAIN (101 

COHHON/0CAL IBR/SLOPE ( 10) #2EROTAP llO 1 

COMMON Nf ILE * NRECOR »LENARR»NCHANN» TEMPI 10*1001 tNRITOAT 

COMMON/BSKPEor/BAOATAiLPACOAT,NftECSKP*NFIL5KP 

DIMENSION SUMAVE (10) #TSUMAVE (101 . SgVALUE <1 0) *ST ANOEV < 10 ♦ 2) * 

. TOTAVE (10*2) .SUMTEMPUO) •TSUMSQ CIO > t ACTUAL! lOl » 

. AVERECIlO) 

nrecor = 0 

BAOATA = 3HTES 
DO 100 IM.NCMANN 
SUM A VE ( I ) = 0.0 
TSUMAVE(I) = 0.0 
Tsuwsom = 0.0 
SQVALUEU) * 0.0 
SUMTEMP ( I ) 3 0.0 
DO 100 K s 1 . N I NSCAL 
STANDEVU.K) = 0.0 
100 TOTAVEU.K) = 0.0 
1NSCAL = 1 
ICHECK = 0 
NSAMPLE « 1 
LEVEL = IOh/ERO INPUT 
200 CALL BUFTEMP 

IF (Nf ILE ,GT. 1) GO TO 1050 
GOTOBUF = 3H NO 
DO 400 I=1*NCHANN 
DO 300 K=1 .LENARR 

300 SUMTEMP ( I ) = SUMTEMP (II ♦ TEMPI I *K> 

AVEPECII) = SUMTEMP! D/LENARR 
400 SUMTEMP (I) = 0.0 

IF (ICHECK .GT. 0) GO TO 800 
IF (NRECOR .EG. I) GO TO 500 

IF ( ABS ( AVEREC ( 1 >I .GT. ABS ( ToT AVC ( l * I nSCAL )> ♦ VARllNl GO TO 700 
500 KHITEI 6 . 3 ) NRECOR* INSCAL .LEVEL 

3 FORMAT ( 1 HO .5* •RECORD M£ANS*<.A*RECORO NUMB€R*I4f7X*CALIBRATlON*12i4 
.X*INPUT *A10/H X«CHANNELM0 X*mEANM 3X*CUMULAT I VE ME AN*6X*NUM8£R Of 
. RECORDS FOR CUMULATIVE MEAN** 

DO 600 I = 1 • NCM ANN 
SUMAVEUl = SUM AVE I I ) ♦ AVLRECUI 
TOTAVEII. INSCAL) * SUMAVE(I) / NSAMPLE 
*R HE (6 *2) 1 .AVEREC M > * TOTAVE ( T * INSCAL ) .NSAMPLE 
2 FORmaTHh *12X*U*10x*F8,4.l4J«*F8.4f25Xf 13) 

00 600 K~ l .LENARR 

600 sovalufi!) = sqvaluE(I) ♦ tempu*k)** 2 
NSAMPLE a NS AMPLE ♦ 1 
GO TO 200 

700 INSCAL a INSCAL ♦ I 
LEVEL = lOhFULL SCALE 
800 ICHECK a ICHECK * 1 

IF (ICHECK .EQ. 1) GO TO 200 
I END = ICHECK - 1 

WRITE 16.41 NRECOR. INSCAL. LtVEL 

4 FORMAT (|H0, SX^TFMPORART SUM OF MEANS*4X*RECQ«D NUMBER* 14 * 7X*CAL 16R 
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SUBROUTINE INSTCAL TRACE COC 6400 FTn V3.Q-P261 0PT=0 02/10/72 08*46.42. 

• AT ION* I2*4X* INPUT *Al 0 /UX#CHaNNEL*10X*$UM OF MEANS* 13**NUMBER OF 
.RECORDS IN SUM*! 

DO 900 1“ 1 *NCNANN 

TSUMAVE I ! ) * TSUMAVE ClI ♦ AVEPECU) 

WRI TE (6*5) I*TSU«AVECn, IEND 

5 FORMAT (IN •13AfI2»16xF6.3t27XT2) 

DO 900 K*WLENARR 

900 TSUMSO(I) = TSUMSOU! ♦ TEMP<ItK)**2 

IF (AB$(AVFRECU >) .GT* A0S < TOTAVE ( | • INSCAL-II I ♦ VAflllN) 

GOTOBUF = 3HYES 
IF ( 1 CHECK «GT • 3! GO TO 1100 
JF f&OTOBUF *EQ, JhYeS) GO TO 200 
DO 1 000 1 - 1 * NCHANN 
TSUMAVEU) s 0*0 
1000 TSUMSQm r o.o 

LEVEL * 10N7ERO INPUT 
ICHECK * 0 
1NSCAL = INSCAL - l 
GO TO 200 

1050 INSCAL = INSCAL ♦ 1 
NSAMPLE = nsample - 1 
1)00 If NO = INSCAL - I 

LEVEL = 10H2ERO INPUT 

IF IIEND .CO, 21 LEVEL a IOhFuLL SCALE 

WRITE <6*61 IENOtLEVEL 

6 FORMAT ( 1H0 ♦ 5X*ST ANDARD 0£V I AT ION$*10X*CAL 1BRAT ION* 12 * 5X* INPUT *A10 

• / 1 1 X*CHANNEL *1 0X*RMS* 1 • 

DO 1200 1=1. NCHANN 

STANDEVd . INSCAL-n * SORT ( SOyALUE III/ (NSAHPLE*LENARRI • XOTAVElIt 
. INScAL-1 ) **2> 

1200 *RITE<6*7M*5TANDEV<I*JEND! 

7 FOPMATllH ,12X12*7XF9.3) 

NSAMPLE = ICHECK - 1 

DO 1300 1=1 ♦NCHANN 
SUMAVE 1 1 ) = TSUMAVEU) 

TOTAVE t I • INSCAL) = SUMAVE < 1 1 /NSAMPLE 
SOVALUE < 1 1 = TSUMSOdl 
TSUMAVEU) ^ 0.0 
1300 TSUMSQU) =0.0 
IChECK = 0 

NSAMPLE = NSAMPLE ♦ ) 

LEVEL = lOHFULL scale 
IF (INSCAL .LE. NInSCAL) GO To 200 
1400 DO 1600 I=I*NCHANN 
DO 1500 K-l fN INSCAL 

1SOO TOTAVEU.M = SLOPE ( I > * TOTAvEUtK) ♦ 2EROTAP ( 1 ) 

ACTUALUI = EXcnvo*<(RESlSU*n/(HE51S<Nl)*CALRE$CItlm - (RES1 
. SU*2)/(hESISU*2| ♦ CALRES< I »2) 1 1 ) 

GAINU) = (TOTAVE < 1 * 2 ) - TOT AvE (I » 1 1 > /ACTUAL (I I 
1600 WRITE (6*8) l»T0TAVt(I*l)*TQTAvE(I*2) • ACTUAL U ) •GAINU I 

8 format tiHO* iox*c h a n n e l*t3/15x*values used for gain calculati 

• ONSMOX* INPUT *SX® TAPE VALUE*5k*ACTUAL VALUE*/59X*ZER0*7XF6«3*12X*0 

• • 0 */56X*FUL L SCALE *4xF6.3*l 1 *F6 . 3// 1 5X *GA I N COmPUTED*F 1 0 .3) 

R£ TURN 

END 
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BUFTEMP TRACE COC 6400 FTN V3.0-P261 OPT«Q 02/IO/T2 00*46.42, PACE I 

SUBROUTINE BUFTEHP 

COMHON/9BUFTEH/IDENTtHUtEOF.lEAITfNTOTFIL«PLOTf IPARITVtZEROTINt 

* JCtOCK 
COMMON/BSKPEOF/BADAT AtLPACDATf NRECSKP*NFILSKP 

COMMON NF IL E ♦ NRECOR *l£NARR*NCh ANN • TEMP 4)0*1001* UR IT DAT 
COMHON/UNPK/ITIME*1COMWRD<200) * IDATWRD U 000 ) 

IF IIDENT • EQ • ' 3HYE5 *A. BAOATA ,EU* 3H NO ,A. NRECOR .EG. 0) 

* CALL HEADER 
CORTlME = 3H NO 
HADAT A s 3H NO 

100 NRECOR * NRECOR ♦ 1 

BUFFER IN ( 1 * l J <1TIME*IC0MW«D<LPAC0AT>> 

IF(UNITM) 1500*200*400 
200 WftITE<6*l> NRECORtNFILE 

1 FORMATdHO** there ARE*IS* RECORDS ON F1LE*I3) 

NRECOR = 0 

IF (MULEOF .NE. 3HVES) GO TO 300 
NFILE = NFfLE ♦ 1 
IF (NFILE .GT- NTOTFILI GO TO 300 
GO TO 100 

300 IF < PLOT *E0 • 3HYES) IEXIT = 3HYES 
IF (PLOT . EQ* 3HTESI RETURN 
CALL EXIT 

400 WRITE <6*21 NRECORtNFlLE 

2 FORMAT ( 1H0* * PARITY ERROR IN DATA* RECORD*IS« FlL£*I3> 

[PARITY * IPARlTY ♦ 1 

WRITE (6*3) 1PARITY 

3 FORMAT ( l HO * * THERE HAVE 6EtN*f3» PARITY ERRORS* > 

CALL UNPAK 

CALL SORT 

IF (WRITOAT ,EQ, 3H NO) CALL DATWRlT ’ 

GO TO 600 
500 CALL UNPAK 
CALL SORT 

600 IF (NRECOR .EO* 1 *A. ITIME *nE. 0) ZEROTIH * ITIHE 
IF ( ITIME-999999 *GT, -12000) CORTlME - 3MYES 
ITIME = ITIME ♦ JCLOCK * 999990 
IF (CORTlME *EG* 3HYES) JCLOCk * JCLOCK * I 
RETURN 
END 
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/ 

SUBROUTINE HEADER 

COMMON NFIlE*NRECOR*LENARR*NChANNiTEMP < | 0> 100) tMRlTOAT 
DIMENSION 1 0 12 > 

BUFFER IN { 1 *0 ) (I0<l>»l0<2> I 
Jf (UNITUM 200* 100*100 
100 *RIT£<fc*l> NFlLE 

\ FORMAT (1H0** PARITY ERROR OR FOF IN HEADER OF FILE NUMBER* 131 
200 W<MTE<A*2> NT lLEt(IDIIl*I»i.2l 
Z FORMATOHOo* HEADER ON F ILE*I3* IS *21101 
RETURN 
END 
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SORT TRACE CDC 6*00 FTN V3.0-P26I OPI*0 02/10/72 00. *6. *2. 

SUBROUTINE SORT 

COMMON NFlLE.NRECGR.LENARRtNCHANN.TENPtlOtlOO) •WRITDAT 
COMMON/BSORT/I0EGSKP»FACTOR*ISKIP 
COHHON/UNPK/ IT IME* I COMORO (200 > i IDATWRO UOOOI 
M = 1BEGSKP 
L = 0 

00 200 1= WLENARR 
DO 100 K- 1 1 NCMANN 

TEMP (Kt 1 1 = IDATWROIM ♦LI * FACTOR 
100 L - L ♦ l 
L * 0 

200 M=M ♦ 1SKJP 

IF ( WR I TDAT .EQ. JMTES) CALL OATHfilT 

RETURN 

ENO 
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SUBROUTINE DATWRIT 

C0«H0N/UNPK/ITlME»ICO*WRD<200t »10ATWRO<1000> 

COMMON NFlLE.NRECOR,LENARR*NCHANNtTEMPUO»IOOI»WRITDAT 
kWRITE1S.ll NRECORt ITIME 

1 FORMAT <]H0t* RECORD NUMBERM5* ITImE IS*110) 

DO 100 1-1 »NCHANN 

100 wRITE16*2> l*<TEMP(I»K>«K-LtLFNARR) 

2 FORMAT {IHO.* TAPE CHANNEL NUMBER*! 3t/» U 6FB. 3) 1 
RETURN 

END 
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SUBROUTINE SKIPEOF 

COHMON/eSKPEOF/aAOATAtLPACOATtNRECSKPiHFlLSKP 

COMMON/BBUFTEM/IDENTfMULEOf*lEXlT.NTOTFIL*PtOT*lPARlTV. 2 EROTIH* 

• JCLOCK 

COMHON/UNPK/ 1 T I HE* ICOHWRO ( 200 ) • 1 DATKRD f 1 0 Q 0 1 

COMMON NFILE'NPECOP*lENAPR'NChAnn»I£HPUOiIOO) •WftlTDAT 

NREC = 0 

IF (NFILE «GT. NFILSKP) go TO 500 
100 BUFFER INU, 1 ) (ITIME, ICOMWRD (lPACDAT t ) 

IF (UNIT dn 300 * 400*200 
200 LEN = LENGTH ( 1 ) 

NREC = NREC ♦ 1 

NRECOR = NRECOR ♦ I 

WR 1 TE( 6 # 1 > NFIlE*NRECOR,NPEC*lEN 

S FORMAT UHQ • 5 **PAR 1 TT ERROR OCCURRED WHILE SKIPPING FILE*I 2 * record 
♦*I 4 *.* 13 * RECORDS HAVE BEEN SKIPPED. LENGTH OF REC 0 R 0 M 4 ) 

GO TO 100 

300 LEN * LENGTH(l) 

NREC = NREC ♦ 1 
NWECOR = NRECOR * 1 

IF ILEN *NE » LPACOAT * \) Wftl TE <6#2INREC0fi #NF I lE f LEN#NREC 

2 FORMAT ( 1 HO, 5 X*RECORO ENCOUNTEpEO OF IMPROPER LENGTH. RECORD* 14 * Tl 

• ILEM 2 * LENGTH* 14 # 2 X 1 3 * RECORDS HAyE BEEN SKIPPED*) * 

GO TO 100 

400 WRITE 16*3) NREC«NF ILE .NRECOR 

3 FORMAT ( 1 HO * 5 X 1 3 * RECORDS HAVE BEEN SKIPPEO ON FILEM 2 *. THERE WER 

• E* 1 4 * RECORDS ON THIS FILE*) 

NFILE = NF IlE ♦ 1 

NREC - 0 
NRECOft “ 0 

IF (IOENT *EQ • 3 HVES) CALL HEAOER 
IF (NFILE .LE. NFILSKP) GO TO 100 
500 IF (NREC 5 KP *EQ* U> RETURN 
DO 900 IsItNRECSKP 

bUFFER IN Cl # 11 UTIME# ICOMWRD (lPACDAT ) ) 

IF (UNIT ( 1 ) > 800 * 700*600 
600 NREC * NREC ♦ l 

NRECOft = NRECOft ♦ 1 
LEN = LENGTH(l) 

WRITE < 6 * 4 ) NRECOR# NFILE #NfttC# LEN 

4 FORMAT | ) HO *SX*PAR|TY ERROR OCCURRED WHILE SKIPPING REC 0 R 0 M 4 * ON F 

• ILE« 12 *. * 2 X 13 * RECORDS HAVE 8 f£N SkIPPEO* LENGTH*| 4 I 
GO TO 400 

TOO WRITE( 6 # 5 !N*FC 0 R#NFlL£«NRfcC 

5 FORMAT < 1 m 0 . 5 X*£ 0 F OCCURREO WHILE SKIPPING RECORDS* LAST RECORD WA 
,S*I 4 * ON F IlE* 12 * 2 X 13 * RECOROS HAD BEEN SKIPPED.*) 

GO TO 900 

BOO NRECOft * NPtCOft ♦ l 
NREC = NREC ♦ 1 
LEN = LENuTH C 1 ) 

IF (LEN.NE. LPACDAT ♦ 1 ) WRITE ( 6 • 2 ) NRECOR #NF ILL* LEN# NREC 
900 CONTINUE 

WRITE ( 6 * 61 NREC *NF ILE* NRECOR 

6 FORMAT <|hO**aI 3 * RECORD (S> HAVE BEEN SKIPPED ON FItE*I 2 *. RECORD 
.NUMBER 1 S*I <0 
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RETURN 

END 



